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Final  report  Grant  No  DAAL03-86-K-0145 
Principal  Investigator:  R.  Lhermitte,  University  of  Miami 

Project  title: 

REMOTE  SENSING  USING  A  GROUND  BASED  94  GHZ 

RADAR 


1. INTRODUCTION 

Tht  original  purpose  of  this  research  was  to  explore  ike  usefulness  and 
performance  of  a  94  GHz  pulse  Doppler  radar  for  the  remote  sensing  of 
clouds  and  precipiiaiion.  The  approach  was  first  to  design  and  assemble  a 
94  GHz  Doppler  radar.  Then  observations  of  clouds  and  precipitation  with 
that  radar  were  conducted  along  with  the  interpretation  of  the  observational 
data  in  an  effort  to  make  contributions  toward  a  better  understanding  of  the 
physics  and  dynamics  of  clouds  and  precipitation  systems. 

This  project  was  originally  based  on  the  design  and  fabrication  of  the 
94  GHz  Doppler  radar  as  well  as  the  implementation  of  hardware  and 
software  methods  for  Doppler  radar  signal  processing.  This  was  initiated 
under  a  previous  contract  but,  during  the  contract  period  covered  by  this 
report,  the  radar  was  substantially  modified  and  upgraded. 

Also,  in  order  to  place  the  94  GHz  radar  in  proper  perspective,  theo¬ 
retical  calculations  of  radar  reflectivity  and  attenuation  coefficient  under 
various  cloud  and  precipitation  conditions  were  performed  for  radar  fre¬ 
quencies  ranging  from  35  GHz  to  240  GHz, 

The  field  experiments  were  all  conducted  with  the  ground-based  9iGHz 
radar  operated  in  a  vertically  pointing  beam  mode.  The  meteorological 
events  observed  with  the  radar  were:  cirrus  clouds,  fair  weather  cumuli 
and  stratiform  and  convective  rain.  The  data  acquired  with  the  radar 
were:  vertical  profiles  of  radar  reflectivity,  vertical  profiles  of  mean  Doppler 
velocity  and  occasionally  Doppler  spectra  observed  at  vertical  incidence. 

This  research  was  concerned  solely  with  ground  based  atmospheric  mea¬ 
surements.  However,  the  experience  acquired  in  gathering  and  analyzing 
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millimeter  wave  radar  observations  presented  the  opportunity  to  explore 
the  capabilities  and  feasibility  of  a  millimeter  wave  radar  installed  on  an 
orbiting  satellite  platform  for  the  purpose  of  remote  sensing  of  earth  cloud 
and  precipitation. 

The  ARO  supported  research  yielded  more  that  ten  published  articles 
which  are  listed  at  the  end  of  the  report  and  numbered  so  that  they  can 
be  used  as  references  in  this  writing  which  essentially  summarizes  various 
aspects  of  the  results.  A  copy  of  some  of  these  articles  is  attached  to  this 
report.  Reprints  of  any  of  the  remaining  papers  is  available  on  request. 

2.RADAR  AND  SIGNAL  PROCESSING 

The  radar,  which  was  developed  in-house  using  commercially  available 
parts,  is  described  in  details  in  [1],  [2],  [3],  and  [5].  The  first  version  was 
based  on  the  use  of  two  antennae,  one  for  transmitting  and  the  other  for 
receiving.  It  was  more  recently  redesigned  to  accomodate  one  antenna, 
using  a  circulator  and  two  crystal  protector  switches  along  with  pulse  timing 
circuits  designed  and  implemented  for  that  purpose. 

The  main  characteristics  of  the  radar  are: 

•  Peak  power:  1.2  kW 

•  Pulse  width  is  adjustable,  but  400  ns  and  100  ns  were  selected  for 
vertical  profiles  and  high  resolution  Doppler  spectra,  respectively. 

•  Pulse  repetition  rate:  10  kHz  for  all  measurements  reported  here. 

•  Receiver  noise:  6.5  db  DSB  noise  figure 

•  Antenna:  90  cm  cassegrain 

•  Doppler  method:  Coherent-on- receive 

•  Doppler  signal  processing:  Pulse-Pair  and  FFT 

•  Minimum  detectable  reflectivity  at  2  km  with  3  s  signal  dwell  time: 
-50  dBZ 

The  radar  receiver  circuits  include  an  orthomode  transducer  yielding 
two  separate  orthogonal  polarization  channels.  The  transmitted  pulse  sig¬ 
nal  polarization  is  horizontal. 
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The  radar  sensitivity  is  expressed  by  the  following  equation: 


IQlogr]  =  lOlogPr  —  IQlogPt  +  lOlogiirR^  —  lOlogh  —  lOlogAe  (1) 

where  g  is  the  radar  reflectivity,  Pr  is  the  received  power,  Pt  is  the 
transmitted  peak  power,  h  is  the  scattering  volume  radial  dimension,  Ae 
is  the  effective  area  of  the  receiving  antenna,  and  R  is  the  target  distance. 
With  the  radar  characteristics  above,  h  —  60m  (400  ns  pulse),  and  a  3.5  dB 
antenna  efficiency  correction,  we  have: 


lOlogT]  =  P,(d6m)  -  23.5  +  20logR 

Or,  in  dBZ: 


(2) 


dBZ  =  Pr{dBm)  +  53  +  20logR  (3) 


At  a  3km  range,  a  cloud  reflectivity  of  -33  dBz  yields  a  signal  ••qual  to 
receiver  noise.  With  a  signal  integration  time  of  3s,  a  radar  reflectivity  less 
than  -50  dBZ  can  be  measured.  Presentation  and  discussion  of  the  radar 
sensitivity  parameters  and  the  radar’s  ability  to  detect  very  small  clouds 
can  be  found  in  more  detail  in  [5],  [6],  [7],  [8j,  [9]. 

The  signal  procfssing  system  includes  a  hardwired  signal  amplitude  in¬ 
tegrator  and  a  hardwired  Pulse  Pair  [8]  processor  delivering  mean  Doppler 
estimates.  The  system  was  designed  to  process  and  record  mean  signal 
intensity  and  mean  Doppler  at  250  gates  spaced  by  400  ns,  thus  covering  a 
15fcm  range  (altitude  for  vertically  pointing).  The  mean  Doppler  and  mean 
signal  intensity  data  were  processed  and  displayed  in  a  format  suitable  for 
analysis  (i.e.  vertical  profiles)  using  a  PC/AT  assigned  to  the  project  and 
equipped  with  a  plotter. 

A  unit  digitally  recording  the  I  and  Q  signals  at  16  selected  gates  simul¬ 
taneously  was  later  added  to  the  radar  data  acquisition  system.  This  unit 
allows  continuous  digitization  and  recording  on  half-inch  magnetic  tape 
(1600  bpi),  of  8192  samples  simultaneously  at  each  of  16  selected  range 
gates.  The  minimum  gate  spacing  is  0.1  fis,  which  is  equivalent  to  a  range 
interval  of  15  m.  At  the  10  kHz  pulse  repetition  rate,  the  recording  time 
(including  data  acquisition  overhead)  is  approximately  1  s  and  a  new  set 
of  16  gates  data  (possibly  automatically  shifted  to  cover  a  greater  range 
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interval)  is  acquired  at  approximately  2  s  time  intervals.  Doppler  spectra 
were  calculated  off  line  from  the  recorded  SEimples  with  a  PC/AT  in  which 
software  developed  FFT  algorithms  were  implemented.  An  example  of  the 
excellent  range  resolution  allowed  by  the  method  can  be  seen  in  (12,  Fig.7]. 

3.  ATTENUATION  AND  BACKSCATTERING  OF  MIL¬ 
LIMETER  WAVE  RADIATION  BY  CLOUD  AND  PRECIPI¬ 
TATION 

The  absorption  and  scattering  of  millimeter  wave  radiation  by  hydrome¬ 
teors  was  investigated  in  an  effort  to  provide  the  general  background  needed 
for  meteorological  interpretation  of  the  94  GHz  radar  observations.  These 
investigations  were  also  used  for  a  compeirison  between  the  performance  of 
the  9iGHz  radar  and  that  of  other  radars  using  different  wavelengths  in  the 
millimeter  wave  part  of  the  spectrum.  This  work  is  presented  and  discussed 
in  [15]  and  the  following  is  a  summary  of  the  approach  and  results. 

Only  spherical  droplets  or  raindrops  (also  spherical  hailstones)  were 
considered  so  that  the  absorption,  scattering,  and  backseat tering  cross  sec¬ 
tions  of  these  targets  could  be  computed  using  the  classical  Mie  functions. 
This  was  done  for  water  and  ice  and  for  a  particle  diameter  ranging  from 
10/im  to  a  large  raindrop  si2e  (6  mm  diameter)  and  at  selected  frequen¬ 
cies  between  35  to  240  GHz,  The  water  and  ice  index  of  refraction  data 
which  were  required  for  these  computations,  were  determined  from  previ¬ 
ously  published  results  (see  [15]).  An  exiimple  of  the  Mie  backscattering 
calculations  at  94  GHz  is  shown  in  Fig.l. 

In  order  to  apply  these  calculations  to  rainfall,  a  Marshall-Palmer  rain¬ 
drop  size  distribution  was  assumed.  This  together  with  the  theoretically  ex¬ 
pressed  Mie  functions  and  a  raindrop  terminal  velocity  vs  diameter  function 
determined  in  an  analytical  form,  allowed  the  determination  of  relationships 
between  precipitation  and  millimeter  wave  radar  observed  parameters. 

More  precisely,  absorption  and  scattering  coefficients  as  well  as  radar 
reflectivity  calculations  were  performed  for  precipitation  of  various  intensity 
and  related  to  rainfall  rate  or  ra,in  liquid  water  content.  This  included 
theoretical  evaluation  of  Doppler  spectra  observed  in  various  precipitation 
conditions  at  vertical  incidence.  The  results  are  reported  in  detail  in  [15]. 
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4.  OBSERVATION  OF  PAIR  WEATHER  CUMULI 

Observing  fair  weather  cumuli  with  the  radar  revealed  its  unique  capa¬ 
bility  to  map  both  the  cloud’s  radar  reflectivity  (down  to  a  -50  dBz  detec¬ 
tion  threshold)  and  the  mean  Doppler  velocity  (beam  vertical  incidence) 
which,  in  that  case  of  very  low  droplets’  terminal  velocity,  can  be  identified 
with  Eur  vertical  velocity  (updrafts  or  downdrafts).  These  observations  are 
presented  and  discussed  in  [9]  and  [11]. 

Observing  these  clouds  with  a  vertically  pointing  beam  yields  vertical 
profiles  of  radar  reflectivity  and  vertical  velocity  continuously  observed  dur¬ 
ing  the  cloud  passage  overhead.  If  there  is  no  time  evolution  of  the  fields 
when  the  cloud  is  passing  overhead,  the  variation  of  mean  Doppler  or  radar 
reflectivity  from  one  profile  to  the  next  only  relates  to  the  fixed  radw  beam 
intersecting  different  parts  of  the  cloud  due  to  the  cloud  motion.  If  the  cloud 
translation  speed,  u,  is  known,  the  time  coordinate  can  thus  be  replaced 
by  a  space  coordinate  x  ~  ut.  If  the  cloud  is  passing  overhead  with  the 
fixed  radar  beam  going  through  the  cloud  center,  a  representative  vertical 
cross  section  of  intra-cloud  radar  reflectivity  and  vertical  velocity  fields  can 
be  obtained.  When  conducting  the  cloud  observations,  we  were  fortunate 
to  find  that  these  conditions  were  met  on  certain  occasions.  Observations 
acquired  in  this  manner  are  presented  and  discussed  in  [5]  and  [11]. 

Examples  of  the  results  are  illustrated  in  Figs. 2a  and  2b.  The  radeu: 
reflectivity  field  is  shown  in  Fig.  la.  The  maximum  reflectivity  is  -  31 
dBZ  and  the  minimum  detectable  reflectivity  -52  dBZ.  There  is  a  sharply 
defined  cloud  top  due  to  a  temperature  inversion  clearly  indicated  by  the 
radiosonde  data  acquired  that  day.  The  cloud  base,  however,  is  not  as  well 
defined  and  exhibits  a  noticeable  variability. 

Inspection  of  the  vertical  velocity  field  cross  section  observed  with  the 
radar,  which  is  shown  in  Fig.2b,  reveals  the  presence  of  a  strong  updraft  re¬ 
gion  (maximum  updraft  velocity  3  ms"').  Comparing  the  velocity  and  the 
radar  reflectivity  fields  indicates  that  the  bulk  of  the  updraft  is  just  below 
(and  slightly  lagging  in  time)  a  region  of  high  radar  reflectivity,  although 
the  downdrafts  are  always  in  a  region  of  very  low  radar  reflectivity.  The 
main  updraft  location  coincides  with  a  local  rising  of  the  cloud  top  (working 
against  the  temperature  inversion  lid),  although  the  downdrafts  are  associ¬ 
ated  with  a  local  lowering  of  the  cloud  bottom  that  can  considerably  distort 
the  cloud  base. 

Maximum  liquid  water  concentration  in  the  cloud  can  be  estimated 
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from  radar  reflectivity  assuming  a  droplet  size  distribution  function,  or  its 
general  shape  and  a  droplet  size  median  diameter  Do.  Simplifying  further 
by  assuming  that  all  the  droplets  have  the  same  diameter  Do  provides  a 
fairly  representative  estimate  of  the  relationship  between  radar  reflectivity 
and  cloud  liquid  water,  M  (in  gmT^).  In  these  conditions,  we  have:  M  = 
0.52i?pZ.  Assuming  Do  ^  10 fim  and  applying  the  above  relationship  to  the 
maximum  radar  reflectivity  observed  in  the  cloud  —SOdB-?),  indicates 
that  the  cloud  liquid  water  in  the  vicinity  of  the  cloud  center  in  a  region  of 
maximum  updraft  is  approximately  0.5  gm~^. 

The  fair  weather  cloud  observations  presented  above,  which  are  more 
extensively  discussed  in  [5]  and  [11],  have  not  been  yet  duplicated  in  any 
other  radar  experiment.  Our  results  indicate  that  such  new  data  can  be 
very  useful  in  the  modelling  of  boundary  layer  “fair  weather”  cloud  circu¬ 
lation  and  possibly  their  dynamics  and  physics  if  radiosonde  and  airplane 
measurements  yielding  observations  of  the  cloud’s  environment  are  included 
in  the  project. 

5,  OBSERVATIONS  OF  HIGH  ALTITUDE  CIRRUS 

Cirrus  clouds  were  systematically  observed  with  the  millimeter  wave 
radar.  Nearly  all  high  altitude  clouds  seen  by  the  naked  eye  are  detected  and 
their  radar  reflectivity  was  found  to  be  between  approximately  —ZbdBZ 
(minimum  detectable  signal  at  a  10  km  range)  and  s:  0  dBZ.  The  mean 
Doppler  velocity  was  always  observed  and  some  of  these  observations  re¬ 
vealed  the  presence  of  strong  upward  and  downward  velocities  (possibly 
reaching  1  ms~^).  Incidentally,  observations  of  stratocumulus  layers  also 
revealed  the  presence  of  well  organized  updraft-downdraft  cells. 

6.0BSERVATI0N  OF  STRATIFORM  PRECIPITATION 

The  general  purpose  of  this  part  of  our  research  was  to  observe  and  study 
the  physics  and  dynamics  of  stratiform  precipitation.  Such  a  precipitation 
system  is  essentially  characterized  by  very  low  air  vertical  velocity,  as  com¬ 
pared  to  convective  rain  in  which  updrafts  and  downdrafts  are  predominant 
features.  Our  approach  was  to  acquire  vertical  profiles  of  mean  Doppler  ve¬ 
locity  and  backscattered  signal  intensity  using  the  94  GHz  Doppler  radar, 
again  operated  in  a  vertically  pointing  beam  mode.  This  work  is  presented 
and  discussed  primarily  in  [10]  and  [12]. 

Examples  of  vertical  profiles  of  signal  intensity  and  mean  Doppler  are 
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shown  in  Figs. 3a  and  3b,  respectively.  Fig.3a  shows  that,  from  the  ground 
up  to  approximately  2.3  km,  the  received  signal  profile  exhibits  a  systemat¬ 
ically  decreasing  intensity  with  a  quasi-linear  slope.  Since  the  mean  rainfall 
intensity  did  not  vary  appreciably  from  the  ground  up  to  the  2  km  level 
(as  indicated  by  the  mean  Doppler  profile  shown  in  Fig.3b),  the  signal 
intensity  decrease  with  height  is  attributed  to  signal  attenuation  by  pre¬ 
cipitation.  The  slope  of  the  function  can  then  be  interpreted  in  terms  of 
the  attenuation  coefficient  alone  which,  in  this  case,  was  found  to  be  8 
dB  km~^  (two-way  path).  This  attenuation  coefficent  measurement  was 
compared  to  measurements  by  other  researchers  and  found  in  good  agree¬ 
ment  with  their  results.  It  was  also  compared  to  the  theoretical  attenuation 
coefficients  based  on  Mie  function  that  we  calculated  (see  above).  These 
investigations  are  discussed  in  more  detail  in  [15]. 

For  the  data  shown  in  Fig.3  the  air  temperature  is  below  freezing  above 
the  2.5  km  altitude  level  and,  beside  the  likely  presence  of  very  small  size 
supercooled  droplets  “unseen  ”  by  the  radar,  the  larger  hydrometeors  found 
there  must  be  of  the  ice  crystals  and  snowflakes  type  (graupel  is  more  likely 
found  in  convective  storms  in  which  updraft  velocity  is  significant).  The 
experience  acquired  in  observing  stratiform  precipitation  with  centimeter 
wave  radars  indicates  that  in  such  conditions,  moving  down  from  an  al¬ 
titude  well  above  the  freezing  level,  radar  reflectivity  increases  suddenly 
due  to  the  partial  melting  of  snowflakes  when  falling  in  the  transition  zone 
where  air  temperature  becomes  above  freezing.  Below  this  altitude,  the 
backscattering  signal  intensity  decreases  due  to  a  lower  particles’  concen¬ 
tration  associated  with  an  increase  of  their  fall  speed  occurring  when  they 
evolve  from  a  wetted  snowflake  structure  into  a  lower  drag  raindrop  shape. 
The  combination  of  these  processes  thus  yields  a  reflectivity  profile  max¬ 
imum  at  an  altitude  slightly  below  the  freezing  level.  This  characteristic 
feature  of  radar  signal  intensity  vertical  profile  in  stratiform  precipitation 
was  named  “bright  band  ”  because  of  the  association  between  signal  inten¬ 
sity  and  radar  display  brightness.  However,  the  bright  band  is  not  observed 
at  94  GHz.  The  drastic  change  in  signal  reflectivity  due  to  snowflake  melt¬ 
ing  appears  clearly  but  there  is  no  vertical  profile  maximum.  Furthermore, 
a  minimum  of  signal  intensity  (“dark  band”)  just  above  the  freezing  level 
is  often  found  in  the  vertical  profiles.  It  was  speculated  that  the  difference 
between  signal  intensity  profiles  observed  with  centimeter  wave  opposed 
to  that  observed  with  millimeter  wave  radars  may  be  attributed  to  the 
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rig.  3.  Presentation  of  10  profiles,  derived  from  data  acquired  on  30 
November  at  5s  time  intervals,  of  signal  intensity  and  mean  Doppler, 
Signal  intensity  is  corrected  for  r’  but  not  for  attenuation  and 
indicated  in  relative  units.  Note  the  well  defined  signal  intensity 
vs  height  slope  equivalent  to  0  dB  km'^  two-way  attenuation.  Mean 
Doppler  data  are  not  corrected  for  air  density. 


peciUarities  of  Mie  backscattering,  i.e.  ein  increase  of  particle  size  may  be 
associated  with  a  decrease  of  signal  intensity.  In  the  case  of  centimeter 
wave  radars,  Rayleigh  backscattering  prevails  which  thus  implies  that  an 
increase  of  particle  size  is  unquestionably  associated  with  an  increase  of 
backscattering  cross  section.  These  considerations  are  discussed  in  more 
detail  in  [11]  and  to  a  lesser  degree  in  [15], 

Fig.3b  shows  the  vertical  profile  of  mean  Doppler.  The  observed  val¬ 
ues  appear  to  be  much  smaller  than  those  observed  with  centimeter  wave 
radars  for  the  same  precipitation  intensity.  This  is  due  to  the  fact  that  the 
Rayleigh  backscattering-raindrop  diameter  relationship,  which  has  a  D® 
form  strongly  enhancing  the  contribution  from  large  raindrops,  does  not 
apply  and  must  be  replaced  by  Mie  scattering  functions.  Indeed,  the  Mie 
functions  emphasize  contribution  from  raindrops  having  a  diameter  near 
the  Mie  maximum  (1  mm)  at  the  94  GHZ  frequency.  The  observed  mean 
Doppler  velocity  data  are  in  agreement  with  the  theoretical  mean  Doppler 
values  determined  in  [15]. 

The  systematic  increase  of  mean  Doppler  with  greater  altitudes  clearly 
seen  in  Fig.3b  is  attributed  to  the  decrease  of  air  density  experienced  when 
moving  up  in  the  atmosphere.  The  slope  of  the  mean  Doppler  function  is  in 
agreement  with  a  theoretical  assessment  of  air  density  effects  on  raindrops’ 
terminal  velocity  [10]. 

Comparing  Fig.3a  and  Fig.3b  indicates  that  the  mean  Doppler  pro¬ 
files  variability  (in  both  altitude  and  time)  is  not  reproduced  in  the  signal 
intensity  profiles  which  are  amazingly  steady.  This  is  attributed  to  the 
fact  that  most  of  the  mean  Doppler  variability  is  due  to  the  contribution 
of  small-scale  air  velocity  (updraft-  downdraft)  fluctuations  which  appear 
surprisingly  high  for  data  acquired  in  such  steady  rain  conditions  [10],  [12]. 

In  addition  to  the  processing  of  spectral  moments  (signal  intensity  and 
mean  Doppler)  discussed  above,  complete  Doppler  spectra  were  calculated 
and  displayed.  This  was  accomplished  by  performing  a  FFT  operation  on 
the  ’’coherent  video”  I  and  Q  signals  delivered  by  the  phase  coherent  radar 
receiver,  simultaneously  sampled  at  16  selected  range  gates.  A  high  speed 
(15  MHz  sampling  rate)  signal  digitizer  was  used  and,  since  there  was  no 
real-time  processing  of  the  radeir  signals,  the  spacing  between  range  gates 
was  reduced  100  ns.  For  these  measurements,  a  100  ns  radar  transmitter 
pulse  was  selected  so  that  a  15m  vertical  resolution  was  effectively  obtained. 

The  Doppler  spectra  calculated  from  the  FFT  procedure  applied  to 
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data  collected  in  10-  20  mmhr~^  stratiform  rain  at  vertical  incidence  using 
this  method  revealed  the  systematic  presence  of  a  well-defined  notch  in  the 
spectrum.  An  example  of  such  observations  is  shown  in  Fig.  4.  The  spec¬ 
tral  notch  was  found  to  be  related  to  the  first  backseat tering  cross  section 
minimum  in  the  Mie  function  at  94  GHZ  which  is  shown  in  Fig.l.  This 
minimum  occurs  at  a  1.6  ram  raindrop  diameter,  and  thus  its  position  in  the 
spectrum  yields  direct  measurements  of  the  vertical  velocity  of  raindrops 
of  that  size.  This  also  allows  an  objective  determination  of  air  velocity 
by  comparing  the  observed  spectrum  to  a  theoretically  predicted  spectrum 
(shown  in  Fig.  4,  see  captions)  based  on  Mie  functions,  a  terminal  velocity 
vs  size  relationship  (corrected  for  the  air  density  at  the  altitude  at  which 
the  observations  are  made),  and  a  Marshall-Palmer  dropsize  distribution. 

This  concept  was  implemented  using  a  more  refined  correlation  method 
[10], [12].  To  avoid  problems  related  to  departure  of  an  actual  dropsize  dis¬ 
tribution  from  the  M-P  model  (especially  for  very  large  or  very  small  rain¬ 
drop  size),  the  correlation  calculations  were  limited  to  a  relatively  small 
diameter  interval  around  the  Mie  minimum  occurring  for  a  1.6  mm  rain¬ 
drop  diameter  in  the  predicted  spectrum.  The  procedure  applied  to  our 
stratiform  rain  data  showed  that  up-  downdrafts  with  a  magnitude  some¬ 
times  reaching  ±  40  cms“*  were  often  observed.  The  correlation  method 
yields  an  objective  measurement  of  vertical  air  velocity  which  can  be  used 
to  correct  the  Doppler  velocity  at  vertical  incidence  so  that  a  dropsize  dis¬ 
tribution,  whose  expression  is  neccessarily  based  on  a  terminal  velocity  vs 
diameter  relationship,  be  derived  from  the  Doppler  spectrum.  An  example 
of  a  dropsize  distribution  obtained  using  this  procedure  is  shown  in  Fig.5. 

The  above  considerations  as  well  as  more  general  results  related  to  anal¬ 
ysis  and  interpretation  of  the  Doppler  spectra  including  some  application 
to  downdraft  measurements  in  convective  storms  are  more  extensively  dis¬ 
cussed  in  [10]  and  [12]. 

6.  94  GHZ  RADAR  ON  A  SATELLITE  ORBITING  PLAT¬ 
FORM 

The  experience  acquired  with  the  94  GHz  Doppler  radar  we  developed 
and  used  extensively  for  observations  of  clouds  and  precipitation  from  the 
ground  provided  us  the  opportunity  to  explore  the  capability  of  such  a  radar 
for  remote  sensing  of  cloud  and  precipitation  from  space.  The  proposal  of  a 
^AGHz  radar  installed  aboard  a  low  orbit  satellite  was  thus  considered  and 
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Fig.  4.  Radial  velocity  (Doj^ler)  spectrum  observed  in  a  15  intvhr  rain  at 
ve'^-ical  incidence  together  with  a  predicted  spectrum  (dashed  line) 
based  on  Hie  back scattering  cross  section  and  an  exponential  dropsize 
distribution  with  the  slope,  A,  adjusted  for  best  fit.  The  shaded 
area  shows  the  difference  between  the  observed  and  predicted 
spectrum.  The  most  noticeable  feature  is  the  spectral  dip  at  6.2  nv/s. 


an  article  presenting  and  discussing  the  radar  concept  and  design  as  well  as 
an  assessment  of  the  radar’s  capabilties  for  probing  hurricane  circulation 
from  space  was  published  [13]. 

7.  CONCLUSION  AND  RECOMMENDATION  FOR  FU¬ 
TURE  WORK 

This  project  has  demonstrated  that  the  ground-based  94iGHz  Doppler 
radar  we  developed  and  used  for  precipitation  and  cloud  remote  sensing  is 
a  viable  tool  for  meteorological  research  and  may  also  be  successfully  used 
for  remote  sensing  from  an  airborne  or  satelliteborne  platform. 

Our  prototype  was  developed  without  any  background  reference  to  any 
otner  similar  system  so  that  we  experienced  growing  pains  with  some  of 
the  components.  A  significant  amount  of  time  was  spent  to  maintain  the 
eqtiipment  in  an  operational  condition,  especially  the  pulse  modulator  part 
of  the  radar  which  failed  often,  sometimes  at  critical  times.  We  believe  that 
this  sort  of  experience  helped  the  manufacturers  of  some  of  the  components, 
such  as  Varian  Associates  and  Pulse  Technology  Inc.,  to  develop  more  reli¬ 
able,  commercially  available,  products.  We  also  believe  that  the  problems 
experienced  in  the  early  stages  of  the  radar  development  have  now  been 
solved.  Therefore,  a  94  GHz  Doppler  radar  of  the  kind  we  developed  can 
be  operated  with  the  reliability  of  a  well-designed  X-band  radar. 

The  coherent-on-receive  method  for  the  extraction  of  the  Doppler  fre¬ 
quency  shift,  which  saves  the  much  greater  cost  of  a  four  cavities  EIA 
(Extended  Interaction  Amplifier),  appeared  very  reliable  and  yielded  sur¬ 
prisingly  low  noise  Doppler  spectra  considering  the  very  high  (94  GHz) 
operating  frequency.  A  large  part  of  this  success  is  due  to  the  excellent 
characteristics  of  the  EIO  (Extended  Interaction  Oscillator)  tube. 

Considering  a  new  equipment,  sensitivity  can  be  increased  by  use  of  a 
larger  antenna,  i.e.  a  180  cm  dish  in  a  Cassegrain  configuration.  Such  an¬ 
tenna  can  be  manufactured  at  a  reasonable  cost  at  the  94  GHz  frequency 
and  would  increase  the  receiver  sensitivity  by  6dB  compared  to  that  avail¬ 
able  with  the  90cm  dish  used  in  our  radar.  In  addition,  the  EIO  peak  power 
can  be  increaised  to  2  it W  and  possibly  slightly  more.  Also  the  receiver  noise 
figure  can  be  improved  by  possibly  2dB.  All  of  these  improvements  would 
yield  a  radar  sensitivity  increase  of  approximately  10  dB  compared  to  that 
of  the  radar  we  developed  originally.  This  would  represent  a  substantial 
improvement  of  small  cloud  detectability.  Our  experience  indicates  that 
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although  fair  weather  cumuli  of  maritime  origin  are  always  detected  and 
observed,  fair  weather  cumuli  of  continental  origin  are  sometimes  not  ’seen’ 
by  the  radar.  This  may  be  due  to  the  fact  that  the  median  diameter  of  the 
population  of  cloud  droplets  in  the  maritime  cloud  is  9-12  pm  compared  to 
7-9  pim  in  continental  cloud.  The  10  dB  increase  in  sensitivity  would  allow 
that  all  fair  weather  clouds  be  observed. 

Observations  of  fair  weather  cumuli  by  the  radar  can  be  significantly 
improved,  not  only  by  the  increase  of  sensitivity  quoted  above,  but  also  by 
a  (approximately)  ±10®  beam  angle  scanning  from  vertical  in  a  direction 
perpendicular  to  cloud  motion.  This  beam  scanning  procedure  would  min¬ 
imize  contribution  to  Doppler  from  the  cloud  horizontal  velocity.  In  addi¬ 
tion,  since  the  scanning  is  limited  to  a  small  angle  from  vertical,  practically 
air  vertical  velocity  alone  is  observed  by  the  radar.  The  beam  scanning 
combined  with  the  cloud  motion  would  then  yield  observation  of  the  fields 
of  reflectivity  and  vertical  velocity  in  three  dimensional  space.  This  is  a 
drastic  improvement  over  the  fixed  beam  method  producing  only  a  vertical 
cross  section  of  the  fields. 

Using  the  equation  of  mass  continuity,  the  dwfdz  term  inferred  field 
from  the  vertical  velocity  field  can  be  interpreted  in  terms  of  horizontal 
divergence  of  the  motion  fields  which,  in  turn  can  be  used  to  specify  the 
cloud  entrainment.  This  procedure,  assisted  by  the  three  dimensionally 
measured  radar  reflectivity  could  serve  as  an  input  for  the  development  of 
cloud  circulation  and  dynamics  models. 

One  of  the  advantages  of  the  radar  is  its  ability  (confirmed  by  experi¬ 
ence)  to  observe  small  nearby  meteorological  targets  without  ground  clutter 
problems,  thus  allowing  extremely  short  range  operation.  This  is  due  to  the 

Rayleigh  backscattering  enhancement  (60  dB  with  respect  to  S-band) 
which  yields  a  high  meteorological  target-ground  clutter  contrast  and  al¬ 
lows  low  power  detection  of  even  weak  weather  targets.  Incidentally  the 
very  short  range  operation  yields  observation  of  extremely  low  reflectivity 
targets  (primarily  above  water)  with  some  of  them  not  explained  by  the 
presence  of  hydrometeors. 

An  application  which  was  not  considered  in  our  research  is  to  observe 
a  heavy  pollution  environment  containing  a  few  pim  size  particulates.  Us¬ 
ing  the  Doppler  capability,  the  air  motion  fields  can  be  determined  leading 
to  an  estimate  of  factors  governing  pollution  occurrence.  We  believe  that, 
especially  with  the  improved  sensitivity,  such  94GHz  Doppler  radar  obser- 
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vations  can  be  done  from  a  minimum  remge  of  50  meters  to  a  few  kilometers, 
if  the  particulate  size  reaches  several  fim's.  These  considerations  can  also 
be  applied  to  the  observation  of  fog. 

We  also  believe  that  the  results  reported  here  can  serve  as  incentives 
for  airborne  (and  possibly  satelliteborne)  use  of  the  94  GHz  radar. 

A  last  word.  Other  millimeter  wave  frequencies  in  atmospheric  ab¬ 
sorption  windows  ranging  from  35  GHz  to  240  GHz  could  be  considered. 
However,  the  94  GHz  frequency  seems  to  offer  a  very  good  compromise  be¬ 
tween  frequencies  which  are  too  low  to  take  a  real  advantage  of  the  Rayleigh 
scattering  enhancement  and  frequencies  which  are  too  high  so  that  signal 
attenuation  problems,  low  transmitter  power,  and  hardware  development 
costs  become  overwhelming.  A  94GHz  radar  for  which  commercially  avail¬ 
able  hardware  exists,  represents  a  good  compromize  and  that  frequency  may 
be  considered  as  an  adequate  choice  for  a  millimeter  wave  meteorological 
radar. 
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ABSTRA(rr 

This  article  is  concerned  with  invest^tions  of  scattering  and  absorption  of  35, 94, 140,  and  240  GHz  radiation 
by  clouds  and  precipitation.  The  computations  of  radar  relkctivi  -y,  Doppler  spectra,  and  absorption  coefficients 
arc  performed  using  Mie  functions  and  ezponenlial  drop>sizr  distributions.  The  results  can  be  used  for  an 
assessment  of  the  applicability  of  short-wavdength  Dopier  rrdars  for  observation  and  study  of  clouds  and 
precipitation  either  from  the  ground  or  from  a  space  platfonr . 


1.  Background 

The  backscattering  cross  section  of  small  scatterers 
(size  much  smaller  than  the  incident  radiation  wave¬ 
length,  \)  is  proportional  to  l/^^  so  that,  compared 
to  that  observed  at  a  3  GHz  frequency  (X  =  10  cm), 
the  radar  reflectivity  of  a  cloud  of  small  droplets  (e.g., 
diameter  less  than  100  /im)  is  greater  by  42  dB  at  35 
GHz,  60  dB  at  94  GHz,  and  74  dB  at  240  GHz.  There¬ 
fore,  radars  operating  at  a  frequency  within  the  milli¬ 
meter  wave  radio  spectrum  (wavelength  interval  I  to 
8.5  mm )  appear  as  more  sensitive  cloud  observing  in¬ 
struments  than  centimeter  wave  radars. 

The  choice  of  a  wavelength  suitable  for  radar  op¬ 
eration  in  the  earth  atmosphere  is  restricted  to  spectral 
regions  in  which  absorption  of  the  radar  signal  by  at¬ 
mospheric  gases,  namely  oxygen  and  water  vapor,  is 
minimal.  The  presence  ofal^rption  lines  at  23  GHz 
(H2O),60GHz(O2 absorption  band),  ll8GHz(02), 
and  183  GHz  (H2O)  restricts  the  wavelength  choice  to 
windows  centered  on  35,  94,  1 40  GHz,  and  between 
200  and  3(X)  GHz.  The  residual  atmospheric  atten¬ 
uation  in  the  windows  is  primarily  due  to  the  skirts  of 
water  vapor  absorption  lines  and  increases  systemati¬ 
cally  with  shorter  wavelengths  ( for  a  complete  view  of 
atmospheric  gas  absorption  in  the  millimeter  wave 
spectrum  see  Liebe  ( 1985)  and  Lhermitte  ( 1987)]. 

The  first  millimeter  wave  radar  designed  for  mete¬ 
orological  use  was  a  35  GHz  system  (see  Paulsen  et  al. 
1970).  The  radar  was  originally  assigned  by  the  US 
Air  Force  to  cloud  deck  monitoring,  essentially  as  a 
replacement  for  the  cmtical  ceilometer.  More  recently. 
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35  GHz  radars  have  been  equipped  with  Doppler  ca¬ 
pabilities  and  used  again  for  the  observation  of  clouds 
and  precipitation  (Hobbs  et  al.  1985;  Pasqualucci  et 
al.  1983;  Omar  and  Sauvageot  1987.  It  was  also  dem¬ 
onstrated  thai  an  even  shorter  wavelength  (\  ~  3.2 
mm  or  94  GHz)  radar  (Lhermitte  1987),  operating 
with  1  kW  transmitted  peak  power  and  3-foot  size  an¬ 
tenna  was  a  viable  meteorolo^cal  tool  providing  ef¬ 
fective  observations  of  liquid  water  content  and  internal 
circulation  in  small  clouds.  Observing  rainfall  at  ver¬ 
tical  incidence  with  the  94  GHz  Doppler  radar  also 
revealed  that  some  of  the  raindrop  size-related  Mie 
backscattering  oscillations  can  be  clearly  seen  and  thus 
can  be  used  for  raindrop’s  diameter  identification 
(Lhermitte  1988b). 

Present  research  and  development  of  millimeter 
wave  tubes  (such  as  the  “gyrotron”  for  example)  is 
expected  to  contribute  to  the  availability  of  more  pow¬ 
erful  transmitters,  possibly  reaching  50  kW  peak  power 
and  50  watts  average  power.  A  94  GHz  radar  equipped 
with  such  a  transmitter  and  a  2  m  antenna  would  have 
a  sensitivity  16  dB  higher  than  that  of  the  94  GHz 
radar  mentioned  above  (Lhermitte  1987, 1988a).  Such 
sensitivity  allows  detection  of  clouds  having  a  radar 
reflectivity  as  small  as  -60  dBZ  (e.g.,  1  mg  m”^  of 
liquid  water  divided  into  8  Mm  droplets )  at  2  km  range, 
and  probably  would  allow  short  range  observation  of 
heavy  pollution  layers  in  which  particulates  may  reach 
a  few  micrometer  size.  Radars  operating  at  240  GHz 
(Narayanan  et  al.  1988)  may  in  the  future  represent  a 
viable  solution  for  the  observation  of  small  clouds,  es¬ 
pecially  ice  clouds.  Therefore  the  characteristics  of 
scattering  and  attenuation  of  240  GHz  radiation  by 
hydrometeors  is  also  considered  in  this  article. 

Millimeter-wave  radars  have  two  distinct  advantages 
over  centimeter-wave  radars: 
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•  Very  narrow  beams  are  produced  with  small  an¬ 
tenna  size. 

•  Rayleigh  scattering  gain  reduces  the  need  for  high 
power  transmitters. 

A  3.2  mm  wavelength  radar,  equipped  with  a  2  m  an¬ 
tenna  and  a  relatively  low  power  ( I  kW  peak)  trans¬ 
mitter,  has  the  same  0. 1  °  antenna  beamwidth  and  sen¬ 
sitivity  on  weak  distributed  targets  (such  as  a  small 
cloud),  as  would  have  a  10  cm  wavelength  radar  with 
a  60  m  antenna  and  a  1  MW  peak  power  transmitter. 
The  difference  in  dimension,  weight,  and  power  supply 
requirements  between  these  two  radars  is  enormous. 
However,  even  if  the  millimeter  wave  radiation  fre¬ 
quency  lies  in  a  spectrum  window,  it  is  significantly 
more  attenuated  than  centimeter  wave  radiation  when 
propa^ting  through  clouds  and  precipitation. 

Observations  of  small  (e.g.,  boundary  layer)  cumuli 
can  be  successfully  conducted  with  millimeter  wave 
radars  without  much  concern  about  signal  attenuation 
by  cloud  droplets  (Lhermitte  I988a)<  However,  deep 
stratus  and  cumulus  clouds  possibly  reaching  the  pre¬ 
cipitation  stage  bring  significant  attenuation  problems 
and  heavy  precipritation  may  attenuate  millimeter  wave 
radiation  so  much  that  the  effectiveness  of  the  radar 
as  a  weather  probing  tool  may  be  significantly  reduced. 
For  instance,  the  successful  operation  of  the  satellite- 
borne  Doppler  radar  mentioned  in  the  conclusion — 
even  if  its  main  capability  is  the  probing  of  motion 
fields  in  precipitation  systems — may  be  severely  limited 
if  the  radar  radiation  does  not  reach  sufficient  low  al¬ 
titude  levels  in  the  storm.  Therefore,  the  evaluation  of 
a  millimeter  wave  radar’s  performance  as  a  weather 
remote  sensing  tool  must  rely  on  assessment  of  the 
physical  conditions  governing  scattering  and  absorp¬ 
tion  of  millimeter  waves  by  clouds  and  precipitation. 

The  scattering  and  absorption  of  both  millimeter 
wave  and  microwave  radiation  by  cloud  droplets  (or 
any  particulate  or  hydrometeor  smaller  than  approx¬ 
imately  100  fim)  can  be  treated  using  the  Rayleigh 
approximation  (scatterers  much  smaller  than  the  radar 
wavelength).  Evaluating  performance  is  easy  since,  in 
these  conditions,  the  familiar  backscatteiing  a,  scat¬ 
tering  Q„  and  absorption  Qa  cross  sections  are  given 
by 

<T  =  ir’/)V(AlA:|^)  (1) 

Q,  =  2ir^D^mX*\K\^)  (2) 

Qa=  (3) 

where  |/^|^  and  lm(- /C)  are  derived  from  the  complex 

index  of  refraction  m  =  n’  +  m"  {n'  is  the  real  part  and 
n'is  the  imaginary  part).  Here  | /Cl ^  decreases  slightly 
with  shorter  wavelengths  and  Im(-A^)  varies  with 
temperature  and  wavelength  (see  Table  I ), 

If  the  particles  are  small  compared  to  the  radiation 
wavelength  (R  tyleigh  limit),  radar  signal  attenuation 
by  scattering  is  always  small  compared  to  the  effect  of 


absorption.  Absorption  of  radiation  by  ice  particles  is 
also  negligible  even  at  millimeter  waves,  but  scattering 
attenuation  by  ice  particles,  while  being  negligible  for 
Rayleigh  scatterers,  becomes  significant  if  larger  par¬ 
ticles  are  considered. 

Although  cloud  particles  (less  than  100  /rm  diame¬ 
ter)  can  be  always  considered  as  Rayleigh  scatterers 
regardless  of  radar  wavelength,  this  does  not  hold  for 
precipitation  particles  and,  at  millimeter  wavelengths, 
relationships  between  precipitation  intensity  and  radar 
reflectivity  or  signal  attenuation  must  be  investigated 
using  Mie  functions.  Even  at  35  GHz,  Rayleigh  back- 
scattering  is  not  rigorously  valid  for  a  particle  diameter 
greater  than  1  or  2  mm.  At  94  GHz  and  above,  scat¬ 
tering  by  all  raindrops  falls  in  the  Mie  region. 

This  article  investigates  the  characteristics  of  ab¬ 
sorption  and  scattering  of  millimeter  wave  radiation 
by  hydrometeors  in  an  effort  to  provide  the  basic  in¬ 
formation  needed  for  an  assessment  of  the  perfor¬ 
mance,  as  a  cloud  or  precipitation  observing  instru¬ 
ment,  of  a  ground-  or  satellite-based  millimeter  wave 
radar.  There  are  scattering  and  absorption  theoretical 
solutions  for  nonspherical  particles  but,  since  the  ge¬ 
ometry  of  these  nonspherical  hydrometeors  (except  for 
cloud  droplets  and  some  raindrops)  is  difficult  to  spec¬ 
ify,  only  spherical  raindrops  or  hailstones  identified 
with  ice  spheres  are  considered  here.  Because  of  the 
dominant  effect  of  surface  tension  for  small  water 
spheres  the  assumption  of  spherically  shaped  raindrops 
is  adequate  for  diameters  up  to  approximately  2  mm. 

In  this  article,  the  absorption,  scattering,  and  back- 
scattering  cross  sections  of  spherical  particles  are  com¬ 
puted  theoretically  using  the  classical  Mie  functions 
(Mie  1908).  The  first  task  was  to  implement  Mie  func¬ 
tion  programs  for  tue  computation  of  scattering  and 
absorption  cross  sections  of  raindrops  (temperature 
ranging  from  0®  to  20®C)  and  ice  spheres  as  a  function 
of  particle’s  diameter,  which  was  done  at  the  35,  94, 
140,  and  240  GHz  frequencies  mentioned  above.  These 
computations  required  specification  of  the  complex 
index  of  refraction  of  water  and  ice  for  all  frequencies 
and  temperatures  considered.  These  data  were  com¬ 
piled  from  several  sources  found  in  the  literature  on 
the  subject  (Ray  1972)^  The  lm(-/^)  and  pa¬ 
rameters  computed  from  the  index  of  refraction  values 
have  also  been  added  to  this  list;  the  results  are  pre¬ 
sented  in  Table  1< 

The  Mie  functions  were  calculated  for  the  frequen¬ 
cies  and  temperatures  listed  above  and  then  stored  for 
use  in  programs  assembled  to  compute  several  radar 
related  precipitation  parameters. 

2.  Drop-size  distribution  and  precipitation  related  pa¬ 
rameters 

In  addition  to  the  Mic  functions  the  programs  de¬ 
signed  to  calculate  radar  related  precipitation  param¬ 
eters  require  specification  of  a  drop-size  distribution. 
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Table  1 .  Index  of  refraction  of  water  and  ice  at  selected  frequencies  and  temperatures. 


Frequency 

(GHz) 

Waier/icc 

Temperature 

(°C) 

n* 

n’ 

Im(-A) 

35 

water 

0 

4.18 

2.58 

0.892 

0.1027 

35 

water 

10 

4.67 

2.78 

0.906 

0.0835 

35 

water 

20 

5.2 

2.87 

0.912 

0.0676 

35 

tcc 

0 

1.8 

0.0001 

0.182 

0.0001 

94 

water 

0 

2.846 

1.48 

0.711 

0.1938 

94 

water 

10 

3.128 

1.75 

0.77 

0.17 

94 

water 

20 

3.41 

2.02 

0.828 

0.147 

94 

ice 

0 

1.78 

0.0004 

0.176 

0.0001 

140 

water 

0 

2.50 

1.05 

0.569 

0.2006 

140 

waier 

10 

2.68 

1.25 

0.643 

0.1952 

140 

waler 

20 

2.86 

t.53 

0.724 

0.190 

140 

tee 

0 

1.76 

0.0004 

0.169 

0.0001 

240 

waier 

0 

2.3 

0.6 

0.414 

0.148 

240 

waler 

10 

2.45 

0.8 

0.494 

0.169 

240 

water 

30 

2.45 

l.t 

0.577 

0,215 

240 

ice 

0 

1.75 

0.0003 

0.169 

0,0001 

In  the  past  considerable  effort  was  devoted  to  selection 
of  an  appropriate  drop-size  distribution  and  the  rain 
parameters  attached  to  it  (Atlas  1 964;  Atlas  and  Wexler 
1963;  Joss  1968;  Sekhon  and  Srivastava  1971;’ Willis 
1984;  Ulbrich  1983;  1988.  However,  the  exponential 
distribution  N(D)  =  Noe\p{- \D),  derived  from  the 
original  Marshall  and  Palmer  (M-P)  (1948)  distri¬ 
bution  is  a  very  convenient  model  and  will  be  used 
throughout  this  work. 

The  M-P  distribution  was  originally  conceived  to 
be  a  unique  function  of  the  rain  intensity,  R,  expressed 
by 


/?  =  (7r/6)3.6  10^ 


[N(D)D^V{D)](iD, 


(4) 


where  R  is  in  mm  h"',  N{D)  is  in  cm■^  D  is  in  cm, 
and  the  raindrop  vertical  velocity,  V(  D)  is  in  cm  s“^ 
In  the  original  M-P  distribution,  N„  is  fixed  and  A  is 
a  function  of  R  only.  However,  an  analytical  solution 
for  A  as  a  function  of  R  and  No  can  found  if  an 
integral  such  as  (4)  has  the  following  form; 

f*Iexp(-AD)f)VZ)]  =  r(rt+  1)/A"^'.  (5) 

Jo 


with  the  right  member  of  the  equation  showing  the 
solution. 

Equation  (4)  can  thus  be  solved  in  close  form  using 
the  V{D)  relationship  proposed  by  Spilnaus  (1948): 
V{D)=  1420D“^,  yielding  the  analytical  solution;  A 
=  82.2^0^^’/?'^^*,  The  Spilhaus  relationship  exhibits 
significant  departures  from  the  measurements  of  ter¬ 
minal  velocity  of  raindrops  by  Gunn  and  Kin/er  (G- 
K)(I949).  The  following  I^(£))  equation  (Lhermitte 
1988b),  is  more  appropriate  since  it  matches  the  G- 
K  experimental  measurements  within  5  cm  s'*  in  the 
50  fim  to  6  mm  raindrop  diameter  range:, 

y{D)  =  t;(l  -  exp(-6.8Z)^  -  4.88D)],  (6) 


where  V(D)  is  in  cm  s"',  D  is  in  cm,  Vg  =  923(po/ 
Pi)°  ’  with  fio  and  p-  being  the  air  density  at  the  ground 
and  at  the  altitude  z  at  which  V{D)  is  expressed,  re¬ 
spectively.  Using  the  K(£>)function  (6)  in  Eq.  (4)(pj 
=  Pa)  and  solving  numerically  leads  to  a  log(  A )  versus 
log(J?)  relationship  slightly  nonlinear  but  in  close 
agreement  with  a  mean  slope  (0.03  to  400  mm  h"‘ 
domain )  yielding  the  solution: 

A  =  69.7(/?/Ao)-®^'..  (7) 

With  No  =  0.08,  (7)  reduces  to  the  familiar  expres¬ 
sion  A  =  41/}'“^'  proposed  by  M-P.  If  the  relationship 
between  log(  A)  and  the  logarithm  of  a  given  parameter 
can  be  represented  by  a  linear  function,  Ng  is  a  constant 
term  independent  of  D.  For  instance,  if  liquid  water 
content,  M  in  g  m'^  is  considered,  A  =  42.1  (A// 
No)~''*  (also  proposed  by  Willis  1984).:  Using  this 
expression  and  Eq.  (7)  we  have: 

=  0.1.33(A„)°“/?°”.  (8) 

With  No  -  0.08  this  reduces  to  A/  =  0.0887/?®*^, 
which  is  in  close  agreement  with  the  expression  pro¬ 
posed  by  Sekhon  and  Srivastava  (1971 ).  In  the  same 
way  the  relationship  between  the  familiar  radar  reflec¬ 
tivity  factor  Z  =  J  [N(D)D^dD]  (in  mm‘  m”^)  and 
/?  is  Z  =  89.7 1  ( which,  with  No  =  0.08 
reduces  to  Z  =  293 

A  relationship  such  as  (8)  for  instance,  or  the  Z-R 
relationship  above,  allows  the  use  of  exponential  dis¬ 
tributions  in  which  No  may  have  a  different  value  from 
that  in  the  M-P  expression,  such  as  the  “thunderstorm 
distribution”  (No  =  0.014  cm'^)  and  the  “drizzle  dis¬ 
tribution”  (Ao  =  0.3  cm“^  proposed  by  Joss  et  al. 
(1968)].  If  the  No  term  is  constant,  the  slope  of  the 
log  Z  versus  log/?  relationship  is  independent  of  Z  or 
/?..  A  large  number  of  empirical  logZ-log/?  relation¬ 
ships  have  been  proposed  in  the  last  40  years.  Within 
the  context  of  Rayleigh  scattering  (also  no  polarization 
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effects)  and  a  linear  exponential  drop*size  distribution, 
a  logZ  versus  log/?  linear  regression  slope  different  from 
the  1.47  value  mentioned  above  may  be  ex  {Gained  by 
a  possible  statistical  dependence  between  No  and  R  in 
a  dataset  including  different  types  of  rain. 

Actual  raindrop  size  distributions  may  not  have  a 
linear  exponential  slope  in  the  full  raindrop  size  domain 
but  rather  exhibit  a  systematic  increase  of  that  slope 
with  raindrop  size,  so  that  the  function  becomes  as- 
symptotic  to  a  value  which  represents  the  largest 
possible  drop  diameter  found  in  that  particular  type  of 
rainfall.  This  effect  is  more  noticeable  for  high  rain 
intensity  and  thus  may  reduce  the  slope  of  the  logZ 
versus  log/?  relationship  on  the  high  end  of  that  func¬ 
tion.  Nevertheless,  as  a  consequence  of  Mie  scattering, 
the  log*)  (i?  is  the  radar  reflectivity)  versus  log/?  (or 
logA/)  fiinction  is  not  linear  at  millimeter  wavelengths. 

Other  drop-size  distributions  have  been  proposed: 
the  gamma  (Ulbrich  1983;  Willis  1984)  and  the  log¬ 
normal  (Fcingold  and  Levin  1986)  functions,  which 
essentially  remove  the  M-P  linear  exponential  slope 
restriction.  However,  at  millimeter  wavelengths  radar 
reflectivity  is  not  overly  dominated  by  contributions 
from  large  size  raindrops  as  it  is  for  Rayleigh  scattering. 
Furthermore,  except  for  very  light  rain  intensity  for 
which,  anyhow,  the  linear  slope  of  the  size  distribution 
seems  to  extend  to  very  small  drop  size,  the  contri¬ 
bution  to  radar  reflectivity  from  drops  having  a  di¬ 
ameter  smaller  than  O.S  mm  is  not  significant.  Indeed, 
as  seen  below,  the  bulk  of  the  contribution  to  radar 
reflectivity  and  attenuation  (both  scattering  and  ab¬ 
sorption)  comes  from  raindrops  with  a  diameter  be¬ 
tween  O.S  and  3  mm,  so  that  a  linear  slope  exponential 
distribution  is  an  acceptable  choice,  requiring  only 
specification  of  No  or  an  assessment  of  the  influence 
of  No  on  the  results.  Furthermore,  the  simple  form  of 
the  exponential  distribution  is  easy  to  manipulate  and 
leads  to  easily  interpretable  results.  In  the  course  of 
this  research  we  concluded  that  a  linear  exponential 
distribution,  in  which  No  is  allowed  to  vary,  is  a  con¬ 
venient  model  yielding  representative  results. 

The  rainfall  parameters  selected  for  computation  are: 
radar  reflectivity,  mean  Doppler  velocity  and  Doppler 
spectrum  width  observed  with  a  vertically  pointing  ra¬ 
dar  beam,  and  attenuation  coefficients. 

These  parameters  are  presented  as  a  function  of  rain 
intensity..  However,  liquid  water  content,  M,  and  pre¬ 
cipitation  parameters  such  as  radar  reflectivity,  v,  and 
radar  signal  attenuation  coefficients  are  related  to  vol¬ 
ume  concentration  of  hydrometeors.  The  concept  of 
rain  intensity  (which  is  really  only  valid  at  the  g;o\ind) 
involves  a  vertical  flux  of  the  liquid  water  at  a  given 
mean  vertical  speed,  vv.  Here  w  is  a  "raindrop  mass 
weighted”  mean  velocity  expressed  by  w  =  /?/A/,  and 
can  be  calculated  assuming  a  given  drop-size  distri¬ 
bution  and  no  air  vertical  velocity..  Using  an  expo¬ 
nential  drop-size  distribution  Ng  exp(— A/)),  w 
=  5250(A)'®’*  is  a  function  of  A  only.  The  vv  varies 


from  2l6  cm  s"*  for  A  =  66  cm~'  (0.1  mm  h”’  M-P 
rain)  to  794  cm  s”'  for  A  =  12  cm"'  (2{X)  mm  h"' 
M-P  rain).  In  the  case  of  significant  up-downdrafts 
(downward  motion  is  usually  taken  as  negative),  Wa> 
/?  =  A/(  w  ±  Wfl).  In  this  perspective,  the  choice  of  rain 
intensity  to  represent  what  is  really  a  volume  concen¬ 
tration  of  liquid  or  solid  water  may  appear  questionable 
when  applied  to  hydrometeors  above  the  ground. 
However,  the  use  of  rain  intensity  as  a  valid  parameter 
was  discussed  in  more  detail  by  Kessler  (1987)  and 
Lee  (1988)  and  their  conclusions  are  used  as  a  justi¬ 
fication  for  the  selection  of  this  quantity  as  a  bonafide 
variable  for  the  expression  of  the  radar  related  param¬ 
eters  in  this  work. 

All  parameters  (including  rain  intensity)  are  nu¬ 
merically  computed  using  an  exponential  drop-size 
distribution  and  when  needed,  the  raindrop  size  versus 
terminal  velocity  relationship  (6)  and  the  Mie  func¬ 
tions. 

3.  Backscattering  cross  sections 

The  Mie  backscattering  cross  section,  (Tmic.  of  water 
spheres  calculated  at  the  35,  94,  140,  and  240  GHz 
frequencies  is  shown  in  Figs.  1  and  2  for  water  at  aO°C 
temperature.  Rayleigh  scattering  does  not  apply  to  the 
full  range  of  raindrops  and  even  the  clearly  defined 
first  Mie  function  minimum  is  within  that  range  for 
all  the  selected  frequencies.  That  minimum  occurs  at 
approximately  D  =  4.5  mm  at  35  GHz,  1.67  mm  at 
94  GHz,  1.10  mm  at  1 40  GHz,  and  0.68  mm  at  240 
GHz.  The  140  and  240  GHz  functions  clearly  show 
that,  for  large  diameters,  the  Mie  oscillations  are  re¬ 
duce  and  that  (Tmic  tends  to  be  equal  to  the  raindrop’s 
geometric  cross  section,  (Tmic  increases  slightly  with 


Fia .  1 .  Radar  cross  section  of  sphencal  raindrops  as  a  function  of 
their  diameter  at  35  and  94  GHz.  The  dashed  lines  indicate  the  Ray' 
leigh  scattering  approximation. 
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temperature,  approximately  0. 1  dB  at  35  GHz,  0.7  dB 
at  94  GHz,  1 .0  dB  at  1 40  GHz,  and  1 .3  dB  at  240  GHz, 
over  the  0®C  to  20‘’C  temperature  range. 

Figure  3,  which  shows  the  ffMie/<^R«yi(>»h  backscatter- 
ing  cross  section  ratio,  indicates  the  maximum  diam* 
eter  value  for  which  Rayleigh  approximation  is  ac¬ 
ceptable  (Rayleigh  function  is  also  indicated  by  the 
dashed  lines  in  Figs.  1  and  2);  i.e.,  approximately,  2.5 
mm  at  35  GHz,  0.9  mm  at  94  GHz,  0.5  mm  at  140 
GHz,  and  0.3  mm  at  240  GHz.  Above  these  diameter 
values  the  Mie  function  oscillates  but,  on  the  average, 
departure  from  Rayleigh  scattering  increases  with 
raindrop  size.  For  a  5  mm  diameter  drop,  fCH'  instance. 


Fig  3.  Ratio  ( in  decibels)  of  Mic  to  Rayleigh  backscatlenng  cross 
sections  at  the  indicated  radar  frequency  as  a  function  of  raindrop 
diameter. 


Flo.  4.  Backacattering  cross  section  of  water  (0°  and  23*0  tem¬ 
perature)  and  ice  spheres  as  a  function  of  diameter.  Frequency  is  94 
GHz. 


the  actual  Mie  scattering  is  approximately  20  dB  below 
Rayleigh  at  35  GHz,  35  dB  at  94  GHz,  45  dB  at  140 
GHz,  and  more  than  55  dB  at  240  GHz. 

Figure  4  shows  the  94  GHz  radar  cross  section  of 
both  water  spheres  at  two  temperatures  0  and  23  ®C, 
tTwaier,  Rnd  ice  spheres,  oic,.  Although,  in  the  Rayleigh 
region,  water  raindrops  are  better  scatterers  than  ice 
spheres  of  the  same  size,  the  effect  is  reversed  for  larger 
diameters  with  ice  particles  now  being  much  more  ef¬ 
fective  scatterers  than  water  particles.  Also  note  the 
relatively  small  variation  of  the  backscattering  cross 
section  of  raindrops  as  a  function  of  their  temperature 
(approximately  I  dB  for  a  0°  to  23°C  temperature 
change).; 

Figure  5  shows  the  ratio,  between  water 

and  ice  spheres  backscattering  cross  sections  at  the  35, 
94,  140,  and  240  GHz  frequencies.  In  the  Rayleigh 
region,  is  cQual  to  the  ratio  of  the  liCp  term 

for  water  and  ice;  i.e.,  -6.95  dB  at  35  GHz,  -6.4  dB 
at  94  GHz,  -5.8  dB  at  140  GHz,  and  -4.7  dB  at  240 
GHz.  This  holds  within  1  dB  up  to  Z)  4  mm  at  35 
GHz,  D  «  1.5  mm  at  94  GHz,  D  1  mm  at  140 
GHz,  and  D  as  0.6  mm  at  240  GHz.  Aljove  these  values 
the  effect  is  reversed  and  ice  become}  a  more  effective 
scattering  agent.  Figure  5  shows  that,  above  this  “crit¬ 
ical”  diameter,  ice  spheres  are  consistently  stronger 
scatterers  than  water  spheres  with  (rjee/ (Twaier  stabilizing 
at  an  average  of  1 5  dB  for  a  particle  diameter  greater 
than  approximately  4  mm  and  for  frequencies  above 
90  GHz.  While  the  general  wisdom  in  radar  meteo¬ 
rology  is  to  assume  that  the  same  amount  of  rainwater 
would  produce  less  scattering  if  that  water  is  converted 
to  ice,  ^is  is  not  systematically  true  at  millimeter  waves 
(frequency  above  90  GHz);  even  relatively  small  ice 
particles  (graupel)  such  as  commonly  found  in  con- 
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as  a  function  of  diameter  at  the  indicated  radar  frequencies. 


vective  storms  may  exhibit  a  radar  reflectivity  greatly 
exceeding  that  of  intense  rainfall. 

4.  Radar  reflectivity  of  rain  at  millimeter  wavelengths 

The  contribution,  dr},  to  the  radar  reflectivity,  r), 
arising  from  dtops  with  a  diameter  between  D  and  D 
+  dD  is  given  by 

dr}  =  N{D)a{D)dD  (9) 

where  <?(£))  is  the  backscattering  cross  section  as  a 
function  of  D.  Assuming  a  M-P  drop-size  distribution 
for  N{D),  dr}ldD  was  calculated  with  the  results  pre¬ 
sented  in  Fig.  6  for  the  35,  94,  140,  and  240  GHz  fre¬ 
quencies,  fora  rain  intensity  of 40  mm  h'‘  ( A,,  =  0.08 
cm  ■*,  A  =  1 9  cm  ‘ ' ).  These  results  reflect  the  behavior 
of  the  Mie  function  and  clearly  show  the  predominant 
contribution  to  radar  reflectivity  of  a  relatively  narrow 


range  of  raindrop  diameters.  The  maximum  contri¬ 
bution  occurs  at  approximately  D  =  2.5  mm  at  35 
GHz,  At  94  GHz,  the  contribution  to  reflectivity  ex¬ 
hibits  a  narrow  peak  around  D  =  I  mm  and  another 
weaker  maximum  dX  D  ^  l.\  mm.  At  140  and  240 
GHz,  several  closely  spaced  maxima  occui  within  the 
drop-size  range  for  which  significant  drop  concentra¬ 
tion  exists,  but  the  average  contribution  extends  from 
approximately  0.3  to  3  mm  diameter.  This  is  partic¬ 
ularly  noticeable  at  240  GHz. 

The  radar  reflectivity  computed  as 


»?  = 


N{D)a{D)dD 


(10) 


using  a  M-P  drop-size  distribution  and  a{D)  values 
derived  from  Mie  functions,  is  presented  as  a  function 
of  rain  intensity  for  the  35,  140,  and  240  GHz  fre¬ 
quencies  in  Fig,  7.  Also  indicated  in  Fig.  7  is  the  cm- 
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FKj.  6.  Contribution  (/■)/<//)  to  radar  reflectivity,  as  a  function  of  raindrop  diameter  in  a  M-P  rainfall  with  A  -  I9(R  =  40  mm  h''). 

pineal  relationship  between  tj  and  R  at  35  GHz  pro-  For  a  0. 1  mm  h"‘  rain  intensity,  the  radar  reflectivity 

posed  by  Wallace  ( 1987)  as  a  least  squares  fit  of  ex-  at  35  GHz  (still  approximately  40  dB  above  that  at  3 
perimental  data.  Wallace's  35  GHz  data  agree  well  with  GHz)  is  approximately  14,  15,  and  16  dB  lower  than 
the  theoretical  computations,  which  supports  the  se-  that  at  94,  140,  and  240  GHz,  respectively.  The  radar 
lection  of  a  M-P  distribution  here.  reflectivity  at  94  GHz,  140  GHz,  and  240  GHz  is  nearly 

The  influence  of  iV„on  the  »?  versus  R  relationship  the  same  for  a  0.3  mm  h~' rainfall  but  still  lOdB  more 
is  shown  in  Fig.  8,  which  presents  the  ij  versus  R  func-  than  that  at  35  GHz.  For  a  10  mm  h"'  rainfall,  the 
lion  at  94  GHz  for  three  A'u  values;  0.0 14  (Joss  “thun-  radar  reflectivity  at  all  the  selected  frequencies  appears 
derstorm  distribution”),  0.08  (original  M-P  distribu-  to  be  almost  identical  (within  2  dB)s  Here  ij  can  be 
tion),  and  0.3  (Joss  ‘‘drizzle  distribution” ).  Below  R  converted  into  the  familiar  equivalent  radar  reflectivity 
=  0.01  mm  h"',  Rayleigh  scattering  dominates  so  that  factor,  Zf,  which  is  defined  by  the  Rayleigh  scattering 
7j  is  proportional  to  ( A’„) Around  R  =  0. 1  value  that  would  produce  the  same  reflectivity, 
mm  h'‘,  n  becomes  fairly  independent  of Nu  and,  above 

I  mm  h“',  the  effect  ofA'^is  reversed  since  an  increase  Z,.  =  10'*?jA^/(|Rl^x*)  (11) 

of  that  parameter  is  now  associated  with  an  increase 

of  1).  A  similar  behavior  is  observed  at  35  GHz,  but  where  Z,.  is  the  mm*  m*^  and  ?;  in  cm“‘.  For  a  100 
the  independence  from  is  reached  for  a  20-30  mm  mm  h"  ‘  rainfall,  Z,  values  are  approximately  40  dRZg 
h’‘  rain  intensity  instead  of  0.1  mm  h“‘  at  94  GH/..  at  35  GHz,  30  dB/,. at  94  GHz,  20 dBZ,, at  140  GHz, 
The  least  sciuares  fit  to  94  GHz  experimental  data  pre-  and  10  dBZ,.  at  240  GHz.  compared  tc  more  than  50 
sented  by  Wallace  ( I9R7 )  is  also  shown  and  is  seen  to  dB2  at  centimeter  waves. 

fall  between  the  =  0.08  and  the  A'„  -  0.3  curves.  Indeed,  departure  from  Rayleigh  seatterii^  at  mil- 
This  suggests  that,  for  the  set  of  data  used  by  Wallace,  limetcr  waves  results  in  compre-.  sion  of  the  radar  re- 
the  mean  A„  is  probably  between  these  two  values.  flectivity  dynamic  range.  For  a  rain  intensity  variation 
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Fio.  7.  Radar  reflectivity  of  a  M-P  rainfall  aa  a  function  of  rain 
intensity  at  the  indicated  radar  frequencies.  The  liquid  water  derived 
from  the  drop-size  distribution  is  also  indicated  in  the  upper  scale. 
The  dashed  line  shows  the  empirical  function  proposed  by  Wallace 
( 198S)  u  a  fit  to  3$  GHz  experimental  data. 


ofO.  1  to  100  mm  h"‘  (associated  with  a 45  dB  variation 
of )» at  centimeter  waves),  tj  varies  by  37  dBat  3 5  GHz, 
22  dB  at  94  GHz,  19  dB  at  140  GHz,  and  16  dB  at 
240  GHz.  These  results  should  not  be  interpreted  as  a 
failure  of  the  short  millimeter  waves  to  probe  a  wide 
range  of  rainfall  intensity  but  rather  as  requiring  better 
accuracy  in  the  reflectivity  measurement.  The  mean 
rf(  log )?)/</(  log  ij)  slope,  which  changes  by  a  factor  of 
two  from  35  GHz  to  94  GHz,  indicates  that  a  1  dB 
uncertainty  in  radar  reflectivity  measurements  pro- 
duces  a  mean  uncertainty  in  the  estimate  of  rain  in¬ 
tensity  of  15  percent  at  35  GHz  and  30  percent  at  94 
GHz.  The  log/{  versus  logti  function  is  not  linear  so 
that,  for  the  same  dB  uncertainty,  the  variance  of  the 
rain  intensity  estimate  is  greater  (smaller)  than  the 
mean  value  for  high  ( light )  rain  intensity.  Incidentally, 
the  smaller  range  of  ij  values  may  remove  the  need  for 
a  logarithmic  compression  receiver,  which  may  be 
considered  as  an  advantage.  Nevertheless,  the  primary 
concern  with  the  reliability  of  t}  quantitative  measure¬ 
ments  at  millimeter  waves  is  signal  attenuatio''. 

5.  Radar  reflectivity  of  Ice  at  millimeter  wavelengths 

The  structure  of  ice  crystals  and  snowflakes  is  com¬ 
plex  and  certainly  far  from  being  spherical  so  that  their 
radar  cross  section  cannot  be  evaluated  theoretically 
using  the  simple  Mie  solution.  The  main  purpose  of 
this  paper  is  to  calculate  radar  reflectivity  on  the  basis 
of  Mie  scattering  functions  and  specified  drop-size  dis¬ 
tributions  and  not  to  consider  empirical  fils  of  obser¬ 
vational  data  except  if  compared  with  theorctical  com¬ 
putations.  Dry,  solid  ice  hailstones  that  have  a  near- 
spherical  shape  can  be  treated  as  icc  spheres,  at  least 
as  a  first  approximation  of  their  scattering  character¬ 


istics.  Assuming  a  known  hailstone  size  distribution, 
radar  reflectivity  can  thus  be  calculated  theoretically. 
To  introduce  the  basic  diflerence  between  hail  and  rain 
backscattcring  at  millimeter  wavelengths.  Fig.  9  shows 
the  density  distribution  of  radar  reflectivity  as  a  func¬ 
tion  of  particle  diameter,  dr\ldD,  for  both  icc  spheres 
(solid  ice,  applicable  to  dry  hailstones)  and  raindrops, 
assuming  an  exponential  size  distribution  for  both  ice 
and  water  with  No  *  0.08  cm”*  and  A  =  15  cm”‘,  At 
35  GHz,  most  of  the  contribution  to  radar  reflectivity 
comes  from  I  to  4  mm  diameter  raindrops  with  water 
drops  having  consistently  more  reflectivity  than  ice 
spheres.  At  94  GHz,  water  raindrops  having  a  diameter 
smaller  than  1 .5  mm  are  still  better  scattcrers  compared 
to  ice  particles  of  the  same  size  but,  for  larger  diameters, 
the  contribution  to  r\  by  ice  significantly  exceeds  that 
from  water,  This  effect  is  more  pronounced  at  140  and 
240  GHz  with  icc  now  showing  an  overwhelming  con¬ 
tribution. 

In  the  data  shown  in  Fig.  9,  the  A  =  15  cm”'  ex¬ 
ponential  slope,  although  smaller  than  the  A  =  19  cm”' 
selected  in  Fig.  6,  still  virtually  prv  ents  particles  larger 
than  D  -  7-8  mm  to  contribute  significantly  to  radar 
reflectivity.  Hailstones  can  grow  to  a  much  larger  size 
than  raindrops.  For  instance,  Douglas  (1964)  reported 
that  observations  of  a  typical  hailstone  size  distribution 
can  be  represented  by  the  following  expression: 

N(D)IM  =  32.2  cxp(-3.09D)  ( 12) 

where  N(D)  is  the  hailstone  concentration  in  number 
per  cubic  meter  and  M  is  the  hailstone’s  water  content 
in  grams  per  cubic  me*'  r.  In  order  to  understand  the 
advantage  of  using  a  size  distribution  normalized  to 


Fig.  8.  Radar  reflecliviiy.  q,  as  a  Tunclion  of  rain  miensily.  R,  for 
exponential  drop-size  distributions  with  diflereni  A’„.  Frequency  is 
94  GHz.  The  dashed  line  shows  the  ij  versus  R  empincal  relationship 
proposed  by  Wallace  on  the  basis  of  94  GHz  expen  mental  data. 
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Fig.  9.  Contribution  to  radar  reflcclivily  (labeled  “ice”)  as  a  function  of  hailstone  diameter  with  a  A  =  IS  M-P  drop-size  distribution. 
The  curves  labeled  “water"  show  the  same  calculations  assuming  raindrops  of  the  same  size. 


water  content,  let  us  assume  N(D)  =  No  exp(-AD) 
and  express  Mas  a  function  of  Nq  and  A  in  the  manner 
discussed  above  and  inserting  the  ice  density,  6.  We 
have; 

A^(£>)/M  =  (exp(-A£))A'')/(5tr).  (13) 

Therefore,  the  primary  advantage  of  the  Douglas’  dis¬ 
tribution  is  that  it  is  independent  of  Nq.  Incidentally, 
comparing  Eqs.  (12)  and  (13)  indicates  that  Douglas 
selected  a  5  =  0.9  value  that  is  acceptable  for  dry  ice. 
The  normalization  to  hail  water  content  is  a  very  useful 
concept  that  can  be  extended  to  calculating  the  ratio, 
tj/M,  represented  by 

)j/M  =  (AV5;r)  )  (exp(- AD)<r(£))</£)]  (14) 

Jo 

which  is  only  function  of  A.  The  tj/M  values  for  ice 
calculated  using  ( 14)  are  presented  in  Fig.  lOas  a  func¬ 
tion  of  A  forthe  35, 94,  140,  and  240GHz  frequencies. 
The  A  domain  extends  from  2  to  10(X).  The  large  A 


values  are  not  realistic  for  hail  (or  even  for  rain),  but 
they  indicate  the  tendency  of  the  functions  to  reach  a 
A"^  slope  arising  from  tj/M  =  0.07|A'|^A'’X"^  an 
expression  that  is  applicable  to  Rayleigh  scattering.  In 
Fig.  10,  tj/M  is  also  shown  for  a  10-cm  wavelength 
assuming  that  all  hailstones  are  Rayleigh  scatterers, 
which  is  really  applicable  only  to  diameters  smaller 
than  approximately  3  cm.  All  the  »j/M  versus  A  func¬ 
tions  at  millimeter  waves  exhibit  a  maximum  for  A 
as  15  cm" ‘  but  this  is  not  the  case  forthe  X  =  10  cm 
function  that  increases  monotonically  for  decreasing 
A  values  so  that,  at  X  =  10  cm,  the  laigest  hailstone 
sizes  produce  the  highest  radar  reflectivity  per  water 
content.  Figure  10  indicates  that,  at  millimeter  wave¬ 
lengths,  ice  spheres  having  a  relatively  small  size  (e.g., 
size  distributions  with  A  =  15  cm"’}  yield  a  greater 
radar  reflectivity  per  liquid  water  content  than  that 
produced  by  the  very  large  hailstones  with  the  size  dis¬ 
tribution  (A  =  3.09  cm"’ )  quoted  by  Douglas  ( 1964). 
At  94  GHz  and  A  =  15  cm"',  tj  =  1.5  10"*  cm"’ 
(approximately  40  and  1 0  dB  more  than  that  observed 
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FH}.  10,  Radar  reflectiviiy  al  Ihe  indicated  radar  fiequencics  of  a 
medium  composed  of  hailstones  with  an  exponential  size  distribution 
N’’  Nc  expf- AZ>),  as  a  function  A.  The  ndar  reflectivity  based  on 
Rayleigh  scattering  at  3  GHz  is  also  shown. 


at  3  and  3S  GHz,  respectively),  a  value  still  10  dB 
greater  than  that  observed  for  rain  with  the  same  A 
and  yVo(  ~200  mm  h‘‘  rainfall  intensity).  In  this  con* 
text,  medium  size  hail  and  graupel  remain  the  strongest 
backscattering  medium  at  millimeter  waves.  Assuming 
a  very  high  value  of  10  g  m"^  for  M,  a  maximum  hail 
radar  reflectivity  of  approximately  2  X  10  cm  " '  may 
be  observed  at  94  GHz  for  a  size  distribution  ap¬ 
proaching  that  of  the  graupel  typically  found  a  few 
kilometers  above  the  freezing  level  in  a  convective 
storm.  However,  for  very  large  size  hail  such  as  quoted 
by  Douglas  (1964),  the  maximum  radar  reflectivity 
will  be  reduced  to  approximately  3  X  10"’  cm"'. 

6.  Doppler  spectrum  at  vertical  incidence 

The  data  shown  in  Fig.  6  indicate  that  the  itmk  versus 
D  oscillations  should  be  observable  using  a  millimeter 
wave  Doppler  radar  probing  rainfall  in  a  vertically 
pointing  radar  mode.  The  Doppler  spectrum  obtained 
in  these  conditions  is  the  density  distribution  of  the 
radar  reflectivity  as  a  function  of  raindrop  velocity.  If 
there  is  no  air  velocity  contribution  the  raindrop’s  ver¬ 
tical  velocity  producing  the  Doppler  shift  is  a  unique 
and  known  function  of  its  diameter,  so  that  the  radar 
reflectivity  density  distribution  can  be  easily  changed 
from  diameter  to  velocity  intervals  and  expressed  as  a 
function  of  velocity  ( Lhermitte  1988b). 

The  Doppler  spectrum  at  vertical  incidence  is  the 
density  of  radar  reflectivity  as  a  function  of  the  particle 
vertical  velocity  F,  A(K)  =  drtfdVand  can  be  expressed 
byA'(K)  =  [N{D)a{D)]dD/dV.  Using  the  Mie  func¬ 


tions,  the  raindrop’s  vertical  velocity  versus  diameter 
function  in  Eq.  (6)  for  V(D),  and  a  M-P  drop-size 
distribution  for  JV(Z))  a  theoretical  representation  of 
the  Doppler  plectrum  observed  at  vertical  incidence 
in  a  rainfall  can  be  achieved.  The  results  for  A  =  19 
cm"'  (40  mm  h"‘  rain  intensity  according  to  a  M-P 
drop-size  distribution)  are  presented  in  Fig.  1 1  for  35 
and  94  GHz,  and  in  Fig.  1 2  for  1 40  and  240  GHz.  The 
Doppler  spectrum  at  35  GHz  fails  to  show  even  the 
first  Mie  scattering  minimum.  The  ^icctrum  at  94  GHz 
clearly  shows  the  first  minimum  at  5.9  m  s"'  and  pos¬ 
sibly  two  others  at  7,9  and  8.8  m  s"‘.  At  140  and  240 
GHz,  two  or  three  minima  can  be  identified  even  iri  a 
noisy  spectrum. 

It  must  be  remenftieted  that,  even  if  there  is  no  mean 
velocity  bias  arising  from  the  presence  of  vertical  air 
velocity,  the  vertical  velocity  of  a  raindrop  may  include 
a  random  component  due  to  air  turbulence  and  oscil¬ 
lation  of  raindrop  shape.  This  should  smear  an  ob¬ 
served  spectrum,  thus  essentially  reducing  the  relative 
amplitude  of  closely  spaced  spectral  minima  and  max¬ 
ima  (Lhermitte  1988b).  Shape  oscillations  are  more 
severe  for  raindrops  having  a  diameter  greater  than 
approximately  two  millimeters.  Hence,  the  spectrum 
smearing  caused  by  shape  oscillations  will  primarily 
affect  the  high  velocity  part  of  the  spectrum. 

Experience  acquired  from  stratiform  rain  observa¬ 
tions  with  a  94  GHz  Doppler  radar  at  vertical  incidence 
indicates  that  the  first  Mie  minimum  can  always  be 
identified  if  the  rainfall  intensity  is  greater  than  ap¬ 
proximately  5  mm  h"'.  Figure  13  shows  a  Dopple' 
spectrum  observed  in  a  8  mm  h"'  rainfall  at  1.2  km 
altitude.  A  spectral  dip  at  6.2  m  s~ ',  in  good  agreemer, 
with  the  terminal  velocity  of  a  1.57  mm  raindrop  at 
the  altitude  considered  ( 1 .2  km ) ,  is  clearly  seen.  These 
results  suggest  that  Mie  oscillations  in  a  Doppler  spec¬ 
trum  can  be  used  to  effectively  test  the  relationship 
between  raindrop  vertical  velocity  and  size,  and  as  a 
method  for  the  measurement  of  both  air  vertical  ve¬ 
locity  and  drop-size  distribution  (Lhermitte  1988h).: 

In  the  case  of  ice  particles  the  first  minimum  of  the 
Mie  backscattering  function  is  not  as  deep  but  should 
still  be  easily  identified  in  a  vertical  incidence  Doppler 
spectrum.  It  appears  at  D  -  4  mm  at  35  Ghz  and  D 
=  1.5  mm  at  94  GHz.  Therefore  a  35  GHz  Doppler 
radar  used  in  a  vertically  pointing  mode  may  be  useful 
in  the  measurement  of  vertical  velocity  of  small  hail 
or  graupel  in  a  convective  storm.  Using  the  method 
for  the  identification  of  large  size  hail  may  be  done 
more  effectively  at  9  GHz,  for  which  the  first  Mie  min¬ 
imum  appears  at  approximately  1.5  cm. 

7.  Mean  Doppler  and  spectral  variance  at  vertical  in¬ 
cidence 

The  mean  Doppler  vertical  velocity,  I  ^ ,  which  is  a 
parameter  useful  for  the  expression  of  mean  velocity 
components  from  Doppler  data,  is  given  by 
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FtO.  1 1 .  Theoretical  Doppler  spectrum  at  vertical  i  ncidenct  based  on  Mie  scattering  functions, 
a  Gunn-Kinzer  raindrop  terminal  velocity  model,  and  a  M-P  drop-size  distribution  with  A 
“  19.  Frequencies  are  35  and  94  GHz. 


where  is  the  drop-size  distribution  and  V{D) 
the  relationship  between  raindrops  diameter  and  fall 
speed.  It  is  clear  from  { 15)  that  is  indejendent  of 
( (5 )  iVo  so  that  it  is  only  a  function  of  A.  Here  Vj  was  cal¬ 
culated  as  a  function  of  A  using  Eq.  (5)  with  an  ex- 


(A 

X 

U 


[N(D)ff(D)V(D)dD]  ^ 

[N{D)a{D)dDl 


Fig.  12.  As  in  Fig.  10,  but  for  140  and  240  GHz. 
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VERTICAL  VELOCITY  ml* 

Pig.  13,  Doppler  ^Kclrum  actually  otnerved  wiih  a  94  GHz 
Doppler  radar  at  vertical  incidence  in  a  8  mm  h  *'  rainfall.  The  Iheo* 
relical  sp^rum  calculated  using  the  Mie  function  and  a  M-P  drop- 
size  disiribution  is  also  shown  by  the  dashed  line. 


ponential  drop  size  distribution,  Eq.  (6)  for  V{D), 
and  Mie  functions  for  For  each  value  of  A,  the 
rain  intensity  normalized  to  N„,  RINo,  and  the  radar 
reflectivity  normalized  to  No,  jj/A'o,  were  also  com¬ 
puted.  The  computations  were  repeated  at  35, 94, 140, 
and  240  GHz,  The  results  are  presented  in  Fig.  14  as 
log(  versus  \og{ZelNo),  with  Z,  given  by  Eq.  (9).. 
The  use  of  Z,  provides  a  normalization  of  ri  so  that  the 
results  are  independent  of  the  selected  frequency  in  the 

Rayleigh  region.  _ 

For  low  radar  reflectivity,  =  2.3 1 Z  ,  which  in¬ 
dicates  that  Rayleigh  approximation  is  applicable  (A 
is  proportional  Z’^’)  and  also  that  for  D  between  10 
fi  and  I  mm,  the  terminal  velocity  of  a  small  size  rain¬ 
drop  is  roughly  proportional  to  its  diameter,  D,  as 
shown  by  the  Gunn  and  Kinzer  ( 1949)  data.  Still  as¬ 
suming  ^yleigh  scattering,  but  for  greater  i^r  which 
V is  not  proportional  to  D,  the  mean  log(  Vd)  ver^ 
logZ,/A'o  is  no  longer  a  linear  function  but  Vd 
ss  2Z  ’ '*■*  is  a  reasonable  This  expression  is  signif¬ 
icantly  different  from  the  Vd  =  2.88Z  obtained  by 
Rogers  (1964)  using  the  Spilhaus  relationship  men¬ 
tioned  before  for  V{D).  The  Spilhaus  function  exhibits 
laige  deviations  from  the  G-K.data  for  both  very  small 
and  very  large  raindrops  and  that  can  explain  the  dis¬ 
crepancy.  On  the  contrary,  the  V{D)  function  (6)  is 
in  much  better  agreement  with  the  G-K  data,  with 
significant  differences  only  localized  in  the  Stokes'  re¬ 
gime  region  (below  Z)  =  10  ;tm )  for  which  the  terminal 
velocity  is  proportional  to  The  log(F,^)  versus 
log(R/Afo)  function  is  shown  in  Fig.  15  together  with 


fig.  14.  Mean  Doppler  velocily  at  vertical  incktence,  (T^),  as  a 
function  of  radar  reflectivity  normalized  to  wavelength  and  the  I 
term;  Z,/A„  al  the  indicate^  frequencies.  Drop-size  distribution  is 
the  M-P  function  and  the  drop  terminal  velocity  is  derived  from  the 
relationship  (6)  in  the  text. 


the  Rayleigh  scattering  solution  applicable  to  centi¬ 
meter  wavelen^hs. 

The  emulations  above  show  that,  for  high  rain  in¬ 
tensity,  Vd  is  significantly  smaller  at  millimeter  waves 
than  at  centimeter  waves,  as  a  consequence  of  the  Mie 
backscattering  cross  section  being  significantly  smaller 
than  Rayleigh  for  large  raindrops.  For  a  rainfall  inten¬ 
sity  above  10  mm  h"',  the  mean  Doppler  value  at  35 
GHz  is  approximately  1  m  s~'  smaller  than  that  cal¬ 
culated  at  centimeter  waves  for  which  Rayleigh  scat- 


RAIN  INTEVSm  {mm/ht) 


FIG .  1 3.  Same  as  Fig.  1 4  bui  expressed  as  a  function  of  ram  imensily 
with  A'o  =  0.08.  The  curve  labeled  Rayleigh  is  Ihe  solution  for  Rayleigh 
scattering  applicable  to  centimeter  wave  radars 
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tering  appli^  At  94, 140,  240  GHz  and  for  high  rain 
intensity,  Vj  is  approximately  3  m  s~‘  smaller  than 
the  Rayleigh  scattering  value  (S  m  s'*  instead  of  8  m 
s'*  for  a  100  mm  h“'  rainfall).- 

_ We  also  computed  the  spectral  variance  given  by 

-  {Vd)^,  Vd^  is  obtained  by  i^lacing  K(i))  by 
F^(D)  in  1^.  ( 1 5 ) .  The  spectral  variance  values,  which 
are  again  independent  of  No,  are  plotted  against  tain 
intensity  assuming  a  M-P  drop-size  distribution  with 
No  =  0.08  (see  Fig.  16).  At  35  GHz,  the  spectral  vari¬ 
ance  is  nearly  the  same  as  that  obtained  using  Rayleigh 
scattering.  At  higher  frequencies,  the  spectral  variance 
is  less  than  the  Raylei^  value  for  rain  intensity  lower 
than  1  mm  h"‘,  but  increases  markedly  for  a  higher 
rain  intensity.  TTiis  can  be  explained  by  the  bimodal 
or  trimodal  nature  of  the  spectra  seen  in  Figs.  1 1 
and  12. 

8.  Millimeter  wave  attenuation  by  atmospheric  gases 
and  hydrometeors 

Atmospheric  absorption  in  the  windows  is  due  pri¬ 
marily  to  the  skirts  of  water  vapor  lines.  For  a  specific 
humidity  of  0.25  g  m  ■^,  the  absorption  coefficient  is 
0.04,  0.042,  0.05,  and  0.09  dB  km”',  at  35.  94,  140, 
and  240  GHz,  respectively..  For  a  specific  humidity  of 
25  g  m“’,  absorption  values  for  the  same  frequencies 
arc  approximately  0.35, 2.1, 5.0,  and  1 3  dB  km“'.  Hu¬ 
midity  is  concentrated  at  low  altitude  in  the  earth  at¬ 
mosphere  so  that  atmospheric  gases  absorption  is  im¬ 
portant  only  for  gound  based  equipment  acquiring  data 
with  low  elevation  angles. 

The  main  consideration  in  the  evaluation  of  the  per¬ 
formance  of  a  millimeter  wave  meteorological  radar  is 
the  attenuation  of  millimeter  wave  radiation  by  clouds 
or  precipitation. 


UAIN  IM  ENSITY  (mrn/l,r) 


Fic.  16.  As  in  Fig.  IS  but  for  the  square  root  of  spectral  variance 


Table  2.  Absorption  cocRkicnts  oohi  3S,  94, 140, 240  GHz  in 
dB  km''  per  g  m~^  of  liquid  water  at  the  indicated  temperatures. 


“OHj 


Temperature 

Am 

O'C 

0.99 

4.9 

7.66 

9.68 

lO'C 

o.go 

4.2 

7.46 

11.05 

20"C 

0.6J 

3.75 

7.26 

14.06 

Even  at  millimeter  wavelengths,  Rayleigh  scattering 
or  absorption  generally  applies  to  cloud  droplets.  This 
implies  that  contribution  from  scattering  to  attenuation 
is  negliyble  and  that  the  attenuation  coefficient,  «« 
=  /  N{D)Q{D)dD,  in  decibels  per  kilometer  per  gram 
per  cubic  meter  (one-way),  is  given  by 

aa  =  8.l8  lm(-/r)/X.  (16) 

The  results  of  these  computations  for  liquid  water  are 
shown  in  Table  2. 

If  the  hydrometeors  are  not  very  small  compared  to 
the  radar  wavelength,  expression  ( 16)  is  not  valid  and 
a  must  be  calculated  using  Qa  from  Mie  functions. 
Also  the  scattering  cross  section,  Qs,  becomes  signifi¬ 
cant  and  signal  attenuation  calculations  must  be  based 
on  the  extinction  cross  section  Q,  Qa  +  Qi-  Figures 
17  and  18  show  the  absorption  and  scattering  cross 
sections  Q,  and  Q,,  respectively,  as  a  function  of  Z>  at 
the  selected  frequencies.  In  the  Rayleigh  region  (D 
<  O.l  mm)  Qo  =  \r^D^  lm(— /ir)]/\  and  thus  is  pro¬ 
portional  to  raindrop  mass.  For  larger  Z)/X  >  2.5,  Qa 
exceeds  the  Rayleigh  absorption  value  and  exhibits  a 
maximum  deviation  from  Rayleigh  for  D/X  s®  2.8  (see 


Fic .  17.  Absorption  cross  section  computed  from  the  Mie  equations 
at  the  indicated  frequencies  as  a  function  of  diameter  for  liquid  water 
raindrops  (temperature  10°C).  The  da^ed  line  shows  the  94  GH? 
solution  in  the  hypothetical  case  of  Rayleigh  scattering. 
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Flo.  18.  As  in  Fig.  17  bul  for  the  scattering  cross  sections. 


dashed  line  indicating  Rayleigh  absorption  at  94 GHz).. 
Above  I>/X  SB  5,  Ca  becomes  markedly  smaller  than 
its  Rayleigh  value.  Indeed,  above  D  =  I  mm  at  94, 
140,  240  GHz,  and  above  D  =  2  mm  at  35  GHz,  Qa 
is  nearly  independent  of  the  radiation  frequency. 

The  raindrop  size  contribution  to  attenuation  is  il> 
lustrated  in  the  data  presented  in  Fig.  19,  which  show 
the  attenuation  coefficient  per  gram  per  cubic  meter 
of  liquid  water  content  as  a  Ainction  of  the  raindrop 
diameter,  calculated  at  35, 94. 140,  and  240  GHz  (0®C 
water  temperature).  For  small  drop  diameter  (less  than 
100  /irn  for  example),  the  attenuation  coefficient  is 
equal  to  the  Raylet^  absorption  values  shown  in  Table 
2.  However,  for  a  greater  raindrop  diameter,  the  atten¬ 
uation  coefficient  increases,  reaches  a  maximum  for 
D/X  as  3,  greatly  exceeding  Rayleigh  absorption,  and 
then  decreases  slowly  to  a  value  smaller  (except  for  35 
GHz)  than  that  calculation  for  Rayleigh  abwrpUon. 
The  attenuation  peak  occurs  at  approximately:  D 
=  2,5  mm  (35  GHz),  D  =  0.9  mm  (94  GHz),  D 
=  0.6  mm  ( 140  GHz),  and  D  =  0.4  mm  (240  GHz). 

The  attenuation  coefficient,  a  given  by 


a=  [Q,{D)N{D)dD]  (17) 

Jo 


was  computed  using  a  M-P  distribution  and  Q,  eval¬ 
uated  from  the  Mie  functions.  The  results,  which  are 
expressed  as  a  function  of  rain  intensity  calculated  using 
the  same  drop-size  distribution  and  the  raindrop’s  ter¬ 
minal  velocity  relationship  (6),  are  presented  in  Fig. 
20  for  a  20”C  water  temperature.  In  the  Rayleigh  region 
a  is  proportional  to  /?“*'*,  which  agrees  with  Eq.  (8). 
At  35  GHz,  the  expression  a  0.2R(Rin  millimeters 
per  hour  and  a  in  decibels  per  kilometer)  applies  to  a 
wide  rain  intensity  domain  and  is  in  good  agreement 
with  observations.  However,  at  94  GHz  and  above,  the 


Flo,  19.  Altcnuation  coefficient  in  dB  km"'  for  I  g  of  liquid 

water  composed  of  the  same  diameter  as  a  function  of  that  diameter, 
at  the  indicated  frequencies.  The  starting  value  (small  diameter)  is 
the  Rayleigh  absorption  value. 


logo  versus  log/?  relationship  becomes  markedly  non¬ 
linear  for  high  rain  intensity,  which  suggests  that  the 
<t  -  aR'’  form  proposed  by  Olsen  et  al.  ( 1978)  is  not 
an  acceptable  model  at  these  very  short  wavelengths. 
Figure  20  also  shows  that,  for  very  high  rain  intensity, 
the  attenuation  coefficient  tends  to  be  independent  of 
frequency.  For  a  10  mm  h'‘  rainfall,  the  theoretical 
value  of  a  varies  from  2  dB  km"'  at  35  GHz  to  ap¬ 
proximately  7  dB  km"'  at  94  GHz  and  8-9  dB  km"' 


T-T  rsmrt  *  ■»  « r»fn»-  -  rtssi*#  i  »  rtMiN  »  »  1 1 Mfi#  t  •»  *  rvim  i  t  rrmm 

!(■<  lO'*  i(|-»  lO"'  U‘  10’  10*  )0‘ 


IMtNSlT^  imm/trl 

Fto.  20  Aiienualion  coefficient  as  a  function  of  ram  intensity  at 
the  35,  94,  140.  and  240  GHz  frequencies.  The  dashed  line  indicates 
the  relationship  proposed  by  Wallace  at  94  GHz  on  the  basis  of  a 
least  squares  fit  to  evpenmenta!  data.  The  arcled  cioss  shows  the  94 
GHz  measurement  by  Lhermitte  (see  text) 
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at  140  and  240  GHz.  These  theoretical  values  of  signal 
attenuation  at  millimeter  waves  are  derived  from  sin^e 
particle  Mie  cross  sections.  Forward  scattering  and 
possible  multiple  scattering  arc  not  included  in  the  G, 
cross  section. 

To  explore  the  validity  of  the  above  treatment,  the 
theoretic^  results  have  been  compared  to  actual  mea¬ 
surements  of  attenuation.  Two  reliable  measurements 
of  signal  attenuation  at  94  GHz  were  considered.  One 
was  a  set  of  attenuation  measurements  reported  by 
Wallace  (1988)  on  the  basis  of  which  a  linear  least 
squares  fit  between  loga  and  log/?  was  proposed,  The 
other  data  was  a  carefully  done  measurement  of  atten¬ 
uation  in  a  stratiform  rain  in  which  the  drop-size  dis¬ 
tribution  and  the  rain  intensity  were  also  measured 
(Lhermitte  1988b).  Both  data  are  presented  in  Fig.  20 
(see  caption).: The  theoretical  calculations  which,  in¬ 
cidentally,  do  not  seem  to  be  significantly  influenced 
by  temperature,  agree  well  with  the  functions  proposed 
by  Wallace  for  very  low  rain  intensity  but  overestimate 
the  Wallace  function  for/?  >  0.1  mm  h“’  (7  dB  km“' 
calculated  instead  of  4  dB  km~'  measured  at  10  mm 
h  “' ),  Wallace’s  data  are  also  in  close  agreement  with 
Lhermitte’s  data  obtained  in  a  8  mm  h^'  steady  rain 
for  which  both  rain  intensity  and  drop-size  distribution 
were  measured  (6  dB  km  '  theoretical  instead  of  3.5 
dB  km"'  measured).  These  results  point  out  the  need 
for  more  experimental  data  on  signal  attenuation  by 
precipitation  at  94  GHz  and  above. 

If  a  can  be  expressed  as  a  function  of  jj,  this  offers 
a  method  for  signal  attenuation  correction.  The  a  ver¬ 
sus  r]  relationship  was  investigated,  assuming  an  ex¬ 
ponential  drop-size  distribution  and  computing  both 
a  and  >)  as  a  function  of  A  and  for  different  values  of 
No.  The  a  versus  ij  functions  obtained  using  this 
method  are  shown  in  Fig.  21  for  a  94  GHz  frequency. 
In  the  Rayleigh  region,  a  is  proportional  to  (A'o)^^’ 
X  which  indicates  a  significant  dependency  on 
No.  However,  Fig,  21  shows  that  the  No  control  on  ij 
is  significantly  reduced  for  higher  tj  values,  and  that 
for  t)  SB  10"*  cm (a  SB  1  dB  km “' )  the  relationship 
between  a  and  »j  is  virtually  independent  of  No.  The 
same  behavior  is  found  at  35  GHz  (not  shown)  but 
the  independence  from  No  appears  for  much  larger 
and  a  value. 

9.  Conclusion 

Because  of  the  typical  low  power  and  antenna  size 
of  a  millimeter  wave  radar,  the  choice  of  a  millimeter 
wavelength  is  very  attractive  for  space  platform  appli¬ 
cations.  A  0.2®  beamwidth  and  a  -20  dBZ  sensitivity 
for  a  spacebome  radar  at  an  altitude  of  300  km  would 
require  a  6  m  antenna  and  100  watts  average  power  at 
15  GHz,  compared  to  1  m  and  10  watts  at  95  GHz. 

Signal  attenuation  by  clouds  and  precipitation  is  a 
serious  problem  associated  with  airborne  or  spacebome 
millimeter  wave  operation  (Meneghini  et  al.  1986). 


Fig,  21,  Allcnualion  coefficient  as  a  function  of  radar  reflectivity 
using  exponential  drop-size  distributions  with  different  Ng.  Frequency 
is  94  GHz. 


Measuring  precipitation  intensity  directly  (based  on 
radar  reflectivity-rain  intensity  relationship)  may  be 
unachievable  at  very  short  wavelengths,  but  attenuation 
is  already  a  problem  at  1 5  GHz  for  instance,  for  which 
the  attenuation  coefficient  is  approximately  0. 1 2  dB 
per  mm  h“‘  of  rain  intensity,  thus  preventing  reliable 
measurements  of  intense  rainfall. 

The  answer  to  the  question  of  feasibility  of  ground 
rainfall  measurements  from  space  may  be  to  derive 
rainfall  intensity  empirically  from  other  storm  param¬ 
eters  such  as  the  altitude  of— and  radar  reflectivity  gra¬ 
dients  at — cloud  tops  and  also  precipitation  penetration 
depth,  using  precipitation  structure  models.  This  ap¬ 
proach  is  somewhat  analogous  to  the  empirical  meth¬ 
ods  used  in  the  treatment  of  satellite  IR  data  to  infer 
earth’s  surface  precipitation.  However,  the  radar  ap¬ 
proach  relies  on  much  more  elaborate  information  on 
backscattering  signal  intensity  gradients  in  the  upper 
storm  region,  rather  than  only  the  upwelling  radiation 
signal  at  the  storm  top. 

Recently,  a  94  GHz  spacebome  radar  was  proposed 
(Lhermitte  1989).  Because  of  its  low  power  require¬ 
ment  and  small  antenna  size,  the  instrument  does  not 
rely  on  complex  technology.  It  is  nevertheless  capable 
of  providing  information  on  motion  fields  inside  pre¬ 
cipitation  systems  that  can  be  penetrated  by  the  mil¬ 
limeter  wave  radiation. 

Providing  some  of  the  background  material  required 
for  evaluation  of  ground-based  and  spacebome  milli¬ 
meter  wave  radar  performance  was  the  primary  incen¬ 
tive  for  this  work.  The  results  yield  a  basic  understand¬ 
ing  of  the  problem  and  are  directly  applicable  to  clouds 
and  relatively  weak  precipitation.  Considering  ice  par¬ 
ticles  with  complex  structure  and  even  large  raindrops 
that  are  basically  oblate  when  falling  will  require  much 
more  experimental  measurements  of  radar  reflectivity 
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and  agnal  attenuation  using  ground-based  and  airborne 
radars. 
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OBSERVATION  OF  RAIN  AT  VERTICAL  INCIDENCE  WITH  A  94  GHZ  DOPPLER  RADAR; 
AN  INSIGHT  ON  MIE  SCATTERING 

Roger  M.  Lhermitte 

Rosenstid  School  of  Marine  and  Atmospheric  Sciences,  University  of  Miami,  Fla 


Abstract,  A  new  method  aimed  at  differentiating  air  verti¬ 
cal  motion  and  raindrop  terminal  velocity  in  a  velocity  spec¬ 
trum  observed  at  vertical  incidence  identification  based  on  Mie 
backscattering  oscillations  occurring  within  the  raindrop  size 
range  at  the  3.2  mm  wavelength.  The  method  is  applicable 
to  measurement  of  downdraft  or  microburst  in  rainshafts  and 
dropsize  distribution  in  all  types  of  rain. 

INTRODUCTION 

Recent  observations  of  stratiform  precipitation  at  vertical  in¬ 
cidence  with  a  Doppler  radar  operating  at  a  94  gHz  frequency 
(As3.2  mm)  have  revealed  a  unique  capability  of  the  system 
based  on  the  fact  that,  at  this  very  short  wavelength,  the 
backscatter'ng  cross  section  of  a  raindrop  as  a  function  of  its  di¬ 
ameter  oscillates  markedly  (Mie  scattering).  Since  they  occur 
within  the  raindrop  size  range,  the  Mie  hackscattering  oscil¬ 
lations  create  a  pattern  of  maxima  and  minima  in  a  Doppler 
spectrum  observed  at  vertical  incidence,  and  thus  can  serve  as 
a  means  to  Identify  raindrop  size.  The  occurrence  of  air  ver¬ 
tical  velocity,  ui,  merely  shifts  the  Mie  oscillation  pattern  and 
does  not  alter  Its  structure.  The  w  shift  can  he  determined  hy 
comparing  the  observed  spectrum  with  a  predicted  spectrum 
evaluated  assuming  w  =  f>  This  note  presents  tbe  new  method 
and  some  of  the  results. 

The  94  GHz  Doppler  radar  hu  the  following  main  charac¬ 
teristics:  peak  power  1  kw;  pulse  width  200  ns;  antenna  size 
3. foot;  beamwidth  0.3“;  minimum  detectable  signal  with  3  s 
integration  time  -115  dBm.  The  radar  is  equipped  with  a  pro¬ 
cessor  providing  vertical  profiles  of  mean  Doppler  velocity  and 
radar  reflectivity  [Lhermitte,  1987].  In  addition,  the  radar  sig¬ 
nal  containing  backscattering  phase  and  amplitude  information 
(coherent'  video),  can  he  sampled  at  a  selected  range  gate  and 
continuously  recorded.  The  recorded  signal  b  processed  hy  a 
micro  computer  programmed  with  FFT  algorithms  to  generate 
power  density  spectra  (Doppler  spectra)., 

DOPPLER  SPECTRUM  AT  VERTICAL  INCIDENCE  IN 
RAIN 

If  a  Doppler  radar  observes  falliag  raindrops  at  vertical  in¬ 
cidence,  tbe  Doppler  frequency  shift  arises  from  the  raindrops’ 
vertical  velocity,  V  .•  V  is  expressed  by  V  =  +  w  ,  where 

Vi  is  the  raindrop  terminal  velocity  in  still  air  and  w  is  the  air 
motion.  The  fall  velocity  of  raindrops  has  been  accurately  mea¬ 
sured  [Gunn  and  Kinzer,  1948;  hereafter  referred  to  as  G-K], 
Although  complicated  functions  bave  been  proposed  as  a  fit  to 
the  G-K  experimental  data  [i.e.  Beard,  1976],  we  have  found 
that  they  can  be  represented  to  a  sufficient  degree  of  accuracy 
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hy:.  V|  =  Vo[l  -  exp  -(oD*  +  6D)],  where  VJ  is  in  cm/s,  D  In 
cm.  With  Vg  =  925  cm/s,  a  =c  6.8,  and  b  s  4.88,  the  standard 
deviation  of  the  above  expression  from  the  G-K  data  Is  less  than 
3  cm/s  in  the  0.5  to  6  mm  raindrop  size  range  and  may  he  con¬ 
sidered  as  an  acceptable  fit,  possibly  reducing  the  experimental 
data  noise.  The  G-K  measurements  were  made  at  sea  level  and 
these  Vi  values  must  he  corrected  for  observations  made  aloft 
hy  {Po/p)°'*  [Foote  and  duToit,  1969]  where  pt  is  the  air  den¬ 
sity  at  sea  level  and  p  is  the  air  density  at  the  height  at  which 
the  observations  are  made.  The  Doppler  spectrum  observed  at 
vertical  incidence,  S{V)  =  n{D).dD/dV.a{V)  can  he  defined  as 
a  density  distribution,  dq/dV  of  the  radar  reflectivity,  t),  with 
the  vertical  velocity,  V,  of  tbe  scatterers.  The  first  term  con¬ 
tributing  to  5(V)  is  the  dropsize  distribution,  n(D)  =  dNjdD, 
a  density  distribution  of  the  number  of  raindrops  per  unit  vol¬ 
ume  as  a  function  of  equal  diameter  D  intervals.  n{D).dD/dV 
is  the  dropsize  distribution  expressed  as  a  function  of  equal 
velocity  intervals.  n(D)  can  he  represented  by  a  simple  expo¬ 
nential  function  n(D)  =  lV,exp(-AD)  proposed  by  Marshall 
and  Palmer,  [1948]  (M-P),  where  Ko  i*  equal  to  0.08  cm~^  and 
A  is  a  parameter  depending  on  rain  Intensity.  Actual  dropsize 
distributions  do  not  agree  fully  with  M-P  as  there  are  signlfl- 
cani,  departures  at  the  low  and  high  ends  of  the  size  spectrum. 
However,  an  exponential  form  is  a  reasonable  fit  to  any  actual 
size  spectrum  within  Dsl  mm  to  3  '  4  mm  (depending  on  rain 
intensity).  (r(v),  is  the  radar  cross  section  of  raindrops  u  a 
function  of  their  vertical  velocity  using:  i.  o  versus  D  from 
Mie  scattering  tables,  and  li.  the  terminal  velocity  V|  versus  D 
relationship. 

We  can  thus  derive  Jb.n(D)  (Jb  Is  an  unknown  constant)  from 
the  Doppler  spectrum  5(V)  alone.  The  evaluation  of  k  re¬ 
quires  measurement  of  the  radar  reflectivity  which  can  he  seen 
as  the  integrzd  of  S(V)  over  all  possible  vertical  velocities.  This 
method,  which  was  proposed  in  the  esirly  days  of  meteorological 
Doppler  radar  methodology,  is  severely  hindered  hy  the  possible 
presence  of  air  vertical  velocity,  w,  which  shifts  the  V{D)  re¬ 
lationship  and  corrupts  particle  diameter  measurements  based 
on  the  Vi  versus  D  function.  For  Urge  raindrops,  even  a  small 
uncertainty  in  w  produces  a  very  large  error  in  D  measure¬ 
ments  (at  D=4.5  mm  the  error  on  D  u  ±0.4  mm  for  a  ±0.2 
ms~*  uncertainty  on  V)).  Two  methods,  originally  conceived 
for  the  purpose  of  measuring  downdraft. updraft  from  Doppler 
spectra,  were  proposed  to  separate  particle  terminal  velocity 
and  air  velocity.  The  first  method  is  based  on  estimating  the 
velocity  spectrum  boundaries  [Battan,  1970];  tbe  second  relies 
on  a  known  relationship  between  the  mean  Doppler  (first  mo¬ 
ment  of  the  Doppler  spectrum)  and  the  radar  reflectivity  based 
on  an  Af-P  assumption  for  n(D)  (Rogers,  1964].  However,  none 
of  these  methods  can  yield  an  acceptable  estimate  of  w  if  the 
actual  dropsize  distribution  is  unknown. 

MIE  SCATTERING  AS  MEANS  TO  MEASURE  W 

Using  a  radar  operating  at  a  verv  short  wavelength  offers 
a  completely  new  approach  to  the  problem  of  separating  V) 
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Fig.  1.  Backicattering  cross  section  of  spherical  liquid  water 
drops  as  a  function  of  their  diameter  at  a  3.2  mm  wavelength 
and  for  0, 10,  20"  temperatures. 


and  w,  based  on  the  presence  of  maxima  and  minima  of  the 
backscattering  cross  section,  a,  occurring  at  well  defined  rain¬ 
drop  diameters.  As  seen  in  Fig.  1,  at  a  3.2  mm  wavelength 
(droplet  temperature  IS’C),  o  exhibits  maxima  at  D  =  1.0, 
2.26,  3.44,  4.62  mm,  and  minima  at  =  1.67,  2.86,  4.04,  and 
6.22  mm.  These  oscillations  should  modulate  the  Doppler  spec¬ 
trum  observed  at  vertical  incidence  and  thus  can  serve  as  a 
means  for  particle  size  idenfilication.  The  occurrence  of  air 
velocity,  w,  will  merely  shift  these  oscillations  in  the  Doppler 
frequency  scale  and  should  not  destroy  their  integrity.  There¬ 
fore,  comparing  an  observed  Doppler  spectrum  with  a  predicted 
spectrum  evaluated  from  the  Mie  function  should  lead  to  a  mea¬ 
surement  of  ui,  and  then,  correcting  for  w,  a  measurement  of 

Recent  rain  observations  with  our  millimeter  wave  radar 
indeed  confirmed  that  some  of  these  Mie  oscillations  appear 
clearly  in  a  Doppler  spectrum  observed  at  vertical  incidence 
in  moderate  to  heavy  rain  conditions.  Fig.  2.  shows  an  ex¬ 
ample  of  the  observed  spectra,  processed  from  data  recorded 
on  November  30  1987,  at  Sudbury,  Massachusetts.  The  ob¬ 
servations  were  made  at  an  altitude  of  1.3  km  in  steady  rain 
conditions  dominated  by  a  radar  "melting  hand”  (see  Fig.4  cap¬ 
tions).  Four  individual  spectra  were  added  to  yield  the  smooth 
spectrum  shown  in  Fig.  2.  Each  spectrum  was  calculated  using 
4096  complex  samples  wbicb  were  recorded  at  1  s  time  inter¬ 
val.  The  dominant  feature  of  the  spectrum  is  the  well  defined 
dip  at  6.2  m/s  (first  Mie  minimum  for  a  diameter  of  1,67  mm), 
surrounded  by  the  two  maxima  at  approximately  4  and  i  m/s. 
Also  shown  in  Fig.  2  is  a  predicted  spectrum  based  on  the 
Mie  function  evaluated  at  the  3.2  mm  radar  wavelength  as- 
suming  an  exponential  drop  size  distribution  with  a  16  cm~' 
slope  and  a  D=0  intercept  adjusted  for  best  match  with  the 
observed  spectrum  in  the  2  to  8  m/s  range  (1  to  3  mm  in  rain¬ 
drop  size).  Tbe  agreement  between  the  two  spectra  in  the  2  to 
8  m/s  velocity  range  is  remarkable.  Tbe  volume  sampled  by 
tbe  radar  (approximately  10^  m^)  is  quite  sufficient  to  yield  a 
representative  sample  of  tbe  dropsize  -listribution. 

Tbe  contribution  from  air  velocity,  w,  can  be  objectively  de¬ 
termined  by  calculating  tbe  cross  correlation,  p(AV)  between 
predicted  and  observed  spectra  as  a  function  of  a  velocity  lag 
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Fig.  2.  Radial  velocity  (Doppler)  spectrum  observed  in  a 
15  mm/hr  rain  at  vertical  incidence  together  with  a  predicted 
spectrum  (dashed  line)  based  on  Mie  backscattering  cross  sec¬ 
tion  and  an  exponential  dropsize  distribution  with  A  adjusted 
for  best  fit.  The  shaded  area  shows  the  difference  between  tbe 
observed  and  predicted  spectrum.  The  most  noticeable  feature 
is  the  spectral  dip  at  6.2  m/s.  The  function  at  the  top  of  the 
hgure  is  the  cross  correlation  between  the  two  spectra  as  a  func¬ 
tion  of  a  velocity  lag.  The  correlation  peak  (0.966)  is  obtained 
for  a  very  small  velocity  lag  (less  than  5  cm/s).. 


AVt  The  result  of  such  calculations  is  shown  in  the  upper  part 
of  Fig.  2.  The  cross  correlation  is  maximum  (0.966)  for  a  ve¬ 
locity  lag  slightly  less  than  5  cm/s  (downward).  Applying  the 
cross  correlation  method  to  a  few  hundred  individual  spectra 
sampled  during  a  10  minute  time  interval,  indicated  that  the 
peak  to  peak  AV  variation  from  spectrum  to  spectrum  was  ap¬ 
proximately  ±30  cm/s  with  most  of  the  observations  showing 
less  than  a  ±10cm/s  value.  These  results  are  consistent  with 
the  observed  variations  of  mean  Doppler  illustrated  by  Fig.4. 
At  the  altitude  of  the  observation,  the  mean  Doppler  velocv 
ity  varies  by  ±20cm/8  although  the  radar  reflectivity  varies  by 
less  than  0.  i  dB.  Such  a  steadiness  of  the  observed  radar  re 
flectivlty  values  indicates  that  tbe  primary  cause  for  variation 
of  the  mean  Doppler  is  variability  in  air  motion.  Since  the 
observations  presented  here  were  made  in  stratiform  rain  con¬ 
ditions,  AV  (attributed  to  au  updraft  or  a  downdraft)  is  small. 
However,  if  thae  measurements  are  made  in  convective  rain 
conditions  (i.e,  in  situations  where  a  rainshaft  is  associated 
with  a  stror.g  downdraft  or  a  microburst)  much  greater  values 
of  w  are  expected. 

In  Fig.  2  the  spectral  dip  in  the  observed  spectrum  does 
not  quite  reach  the  minimum  shown  in  the  predicted  spectrum. 
Spectral  averaging  required  to  yield  a  smooth  spectrum  could 
produce  this  effect,  but  in  our  case  it  is  done  with  a  8  cm/s 
bandwidth,  too  small  to  produce  significant  spectrum  smear¬ 
ing;  also  amplitude  white  noise  is  removed  in  spectrum  pro¬ 
cessing.  The  possibility  of  a  significant  contribution  from  radar 
frequency  noise  was  rejected  because  backscattering  from  fixed 
targets  was  always  observed  to  be  associated  with  very  nar¬ 
row  spectra.  Due  to  the  radar  narrow  beam  (6  10~^  radian) 


Lbermitte:  Rain  at  Vertical  Incidence  with  a  94  GHz  Doppler  Radar 


1127 


7 


VERTICAL  VELOCITY  ml* 

Fig.  3.  Same  as  Fig.  2  but  including  a  5  cm/s  velocity 
shift  and  convolving  the  predicted  spectrum  with  a  20  cm/s 
standard  deviation  noise  spectrum. 

the  combined  horizontal  wind-beamwidth  effect  on  spectrum 
smearing  was  also  rejected  as  well  as  a  possible  contribution 
due  to  multiple  scattering.  It  was  felt  that  the  primary  contri¬ 
bution  to  the  excess  spectral  density  observed  in  the  spectral 
dip  was  the  pretence  of  a  random  component  attached  to  either 
the  velocity  or  the  diameter  variables,  essentially  smearing  the 
V  versus  D  relationship.  Departures  from  sphericity  typically 
observed  for  large  raindrops  during  their  fall  could  produce 
a  random  modulation  of  the  Mie  backscatterlng  cross  section. 
However,  this  effect  Is  not  likely  to  be  important  for  a  drop  with 
a  diameter  less  than  2  or  3  mm,  so  that  random  fluctuations 
of  raindrop  velocity  due  to  air  turbulence,  or  any  variations  of 
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Fig.  5,  Spectrum  as  a  function  of  equal  raindrop  diame¬ 
ter,  D,  intervals.  Both  the  observed  spectrum  and  a  predicted 
spectrum  (dashed  line)  based  on  the  Mie  function  and  an  ex¬ 
ponential  dropsize  distribution  are  shown.  The  shaded  area 
highlights  the  small  difference  between  the  two  functions.  The 
agreement  is  remarkable  between  D  =  0.5  to  D  =  3mm,  but  for 
diameters  larger  than  approximately  3mm  the  observed  spec¬ 
trum  falls  off  more  rapidly. 

w  during  the  signal  dwell  time  may  be  a  more  general  cause 
for  spectrum  smearing.  Such  an  effect  can  be  reproduced  by 
convolving  the  original  spectrum  with  a  spectrum  simulating 
the  random  velocity  fluctuation  and  its  variance.  This  method 
was  applied  to  the  predicted  spectrum  shown  in  Fig.  2.  The 
results  are  presented  in  Fig.  3  and  indicate  that  a  convolution 
with  a  spectrum  having  a  30  cm/s  bandwidth  provides  a  good 


Fig.  4.  Vertical  profiles  of  mean  Doppler  velocity  and  radar 
reflectivity  in  dBZ  acquired  on  30  November  1987,  The  sudden 
velocity  increase  at  an  altitude  of  2500  meters  is  due  to  particle 
melting.  Four  profiles,  4  s  apart,  are  shown.  The  Doppler 
spectrum  measurements  were  made  slightly  above  the  bottom 
dash  line.  Note  the  variability  of  mean  velocity  compared  to 
the  steadiness  of  radar  reflectivity. 


Fig.  6.  a)  Observed  spectrum,  b)  <t(D)  function,  and  c) 
dropsize  distribution  n[D)  all  shown  in  log  scale  (decibels).  The 
approximate  slope  (in  cm"')  of  n(D)  is  indicated  in  several 
parts  of  the  function. 
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match  between  experimental  and  theoretical  results.  Inspec¬ 
tion  of  Fig.  3  also  Indicates  that  velocity  spectrum  smearing 
produces  a  merging  of  the  closely  spaced  third  and  fourth  Mie 
maxima  {D  —  3.6  and  4.6  mm)  and  yields  a  much  closer  similar¬ 
ity  between  observed  and  predicted  spectra.  This  is  due  to  the 
very  small  value  of  dV(/dD  for  large  raindrops.  Application  of 
the  method  to  several  spectra  indicated  that  the  standard  devi¬ 
ation  needed  to  match  theoretical  and  experimental  data  varies 
from  10  cm/s  to  40  cm/s  and  seems  to  increase  with  the  num¬ 
ber  of  individual  spectra  used  in  the  evaluation  of  an  average 
spectrum.  Increasing  air  motion  turbulence  would  ultimately 
result  in  filling  the  spectral  dip  at  1.67  mm  completely  thereby 
preventing  air  velocity  measurement.  This  would  happen  if  the 
random  velocity  fluctuation  reaches  a  peak-to-peak  value  of  ±5 
m/s  (standard  deviation  greater  than  approximately  2  m/s). 

After  in  is  estimated  using  the  above  correlation  method 
and  the  air  turbulence  correction  introduced,  we  can  calculate 
the  spectrum,  5{D)  =  dr)fdD,  now  expressed  as  a  function  of 
equal  diameter  intervals.  An  example  of  such  results  is  shown 
in  Fig.  5  where  the  observed  spectrum  is  again  shown  with  a 
predicted  spectrum  based  on  the  Mie  hackscattering  function 
and  a  linear  exponential  drop  size  distribution  with  the  slope 
adjusted  for  best  fit.  S(D)  =  n(Z?)<T(jD)  so  that  n(D)  can  he 
readily  computed.  An  example  of  the  results  is  shown  in  Fig.  6 
where  S{D),  o(D),  and  the  calculated  n(D)  function  are  shown 
using  a  logarithmic  scale  for  better  perception.  Incidentally,  we 
found  that  the  smoothness  of  the  calculated  n(D)  function  in 
the  vicinity  of  the  spectral  dip  is  very  sensitive  to  a  correct 
choice  of  both  AV  and  the  spectral  variance  added  to  the  pre¬ 
dicted  spectrum.  Fig  6  shows  that  there  is  a  large  deficiency 
of  raindrops  helow  0.5mm.  Above  D=  0.5  mm  the  n(D) 
function  exhibits  a  near  exponential  decay  hut  the  local  expo¬ 
nential  slope  (in  cm”' )  varies  markedly  from  20  (Ds  0.5  to  1.5 
mm),  to  12  (D=  1.5  to  2.5  mm),  to  21  {D=  2  to  4  mm),  and 
increases  more  rapidly  afterwards. 

CONCLUSION 

We  intend  to  apply  the  method  to  measurements  of  air  ver¬ 
tical  velocity,  and  possibly  air  turbulence  and  dropsize  distrihu- 
tion  in  rainshafts  prouuced  hy  convective  storms.  The  method 
would  also  he  very  useful  in  providing  vertical  profiles  of  drop- 
size  distribution  and  air  vertical  velocity  in  precipitation  re¬ 
gions  observed  hy  an  S-hand  radar  in  an  effort  to  improve  our 
knowledge  of  the  radar  reflectivity-rain  intensity  relationship. 
We  also  intend  to  extend  our  measurements  to  the  spectrum 
of  the  cross  polarized  component,  a  method  which  would  bring 
an  insight  of  large  particle  shape  distortions.  So  far  we  have 


limited  our  observations  to  liquid  raindrops  and  it  is  unknown 
at  this  pmnt  if  tbe  method  u  applicable  to  identification  of 
graupel  or  hailstones  and  their  sizes.  However,  we  will  attempt 
to  obtain  Doppler  spectra  in  these  conditions. 
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paptr  U  coneerntd  with  th*  mh  of  short  milllmeitr- 
w«vi  Doppltr  radars  for  tht  obsarvetka  of  clouds  and  pracipllallon. 
Alianuailon  and  scaiiaHni  llneludlng  Mia  batkscoitarini  by  rain* 
drops)  of  ihasc  short  wavalangihs  radiation  by  hydromalars  la  dis* 
oissad  as  wall  as  tha  saasltlvliy  of  such  radars  fbr  tha  obsarvatlon  of 
clouds. 

Rasults  of  obMrvsiloAS  Indkala  that  tha  latarnal  circulation  of  Ihlr 
waaihar  cumuli  b  alTKtlvaly  obMrrad. 


I.  Backoxound 

I^ETEOROLOGICAL  nidan  operate  typically  at  cea* 
IVltimeter  wavelengths  primarily  to  avoid  problems  In 
the  Interpretation  of  echo  signal  Intensity  In  terms  of  radar 
reflectivity,  arising  from  the  slgntflcant  attenuation  of 
shoner  wavelength  radiation  propagating  through  precip¬ 
itation  regions.  However,  when  It  became  evident  that 
clouds  were  not  observed  by  centimeter-wave  radars,  the 
search  for  an  increased  radar  .'sensitivity  was  directed  to¬ 
ward  the  use  of  shorter  wavelengths  for  which  an  increase 
of  cloud  radar  reflectivity  due  to  Rayleigh  scattering  (see 
(3))  is  expected,  and  the  use  of  mllllmeter-wave  radars 
was  recommended. 

Wavelength  selection  for  a  radar  operated  In  the  earth’s 
atmosphere  is  restricted  to  spectral  regions  where  absorp¬ 
tion  by  atmospheric  gases  Is  low.  The  nillimeier-wave 
absorption  spectrum  is  shown  in  Fig.  I.  Outside  of  the 
water  vapor  line  at  33  GHz.  the  spectrum  is  dominated  by 
three  strong  absorption  regions;  (a)  the  60*GHz  complex 
of  rotational  0:  lines,  (b)  anotherO;  line  at  I  IS  GHz,  and 
(c)  a  strong  water  vapor  line  at  IS3  GHz.  The  spectral 
windows  where  local  minima  of  atmospheric  absorption 
occur  are  centered  on  33.  100,  150,  and  310-290  GHz. 
The  residual  absorption  in  all  the  windows  is  primarily 
due  10  the  skirts  of  water  vapor  lines  and  increases  sys- 
tetitatically  with  the  window  frequency.  At  low  humidity 
(less  than  1  g/m\  for  instance),  the  one-way  absorption 
in  all  the  windows  is  below  0.1  dB/km,  but  for  strong 
humidity  (i.e.,  35  g/m’)  it  increases  to  0.7,  2.5,  6,  and 
15  dB/km,  at  35,  lOO,  150,  and  210  GHz.  respectively. 

Although  there  were  attempts  to  observe  clouds  at  a  1 .2- 
cm  wavelength  in  the  early  days  of  radar  meteorology,  the 
first  systematic  millimeter-wave  radar  observations  of 
clouds  were  done  in  the  35-GHz  window  (i.e. ,  ( 10]).  The 
radar  was  designed  as  a  cloud  monitoring  instrument  and 
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was  used  for  several  years  primarily  as  a  qualitative  in¬ 
dicator  of  cloud  cover  (ceilometer),  There  was  no  attempt 
to  measure  cloud  radar  reflectivity  at  these  times  and  it  is 
regrettable  that  the  radars  were  not  properly  calibrated. 
Measuring  cloud  echo  Intensity  would  have  produced  an 
extensive  documentation  on  cloud  radar  reflectivity  at  35 
gHz,  especially  for  cirrus  clouds,  which  were  well  de¬ 
tected  by  the  radar. 

More  recently,  3S-GHz  radars  were  again  considered 
for  cloud  observations  and  some  of  the  available  equip¬ 
ments  converted  Into  Doppler  radars  [12],  (9],  (I).  These 
radars  are  equipped  with  the  same  magnetron  that  was  de¬ 
veloped  eariler,  now  used  in  conjunction  with  a  Stalo- 
Coho  technique  for  measurement  of  the  echo  phase  con¬ 
taining  the  Doppler  information. 

The  expectation  that  a  shorter  wavelength  would  in¬ 
crease  cloud  radar  reflectivity  fbii.h'r  was  an  Incentive  for 
proposing  the  use  of  frequencies  beyond  35  GHz.  The  re¬ 
cent  development,  at  94  GHz.  of  a  receiver  technology 
comparable  to  that  available  for  centimeter-wave  radar 
design  and  the  availability  of  long-life  high-power  trans¬ 
mitter  tubes  (extended  interaction  oscillator,  ElO)  pre¬ 
sented  the  opportunity  to  design  and  assemble  a  cloud  ob¬ 
servation  radar  working  at  a  wavelength  shoner  than  the 
previous  35-GHz  band  limit  (6).  It  was  decided  to  equip 
the  radar  with  a  modem  signal  processing  system  so  that 
accurate  quantitative  measurements  ot‘  radar  reflectivity, 
mean  Doppler  velocity,  and  velocity  spectrum  variance 
could  be  obtained.  Frtmurily  due  to  the  use  of  the  ElO 
and  a  well-designed  low-noise  local  oscillator  (phase  lock 
oscillator),  the  radar  demonstrated  e.xcellent  Doppler  per¬ 
formance  (3].  Some  of  the  results  ot  cloud  observations 
performed  with  the  radar  were  published  previously  (3). 
[4].  This  papei  summarizes  the  radar  characteristics  and 
performance  as  a  cloud  observing  instrument,  discusses 
signal  attenuation  by  water  vapor  and  hydrometers  at  94 
GHz.  and  presents  some  of  the  results  on  cloud  and  pre¬ 
cipitation  observations. 

II.  Radar  Wavelength  ano  Cloud  OesERVATtoN 
Capability 

Assuming  that  the  radar  beam  is  completely  filled  by 
the  distributed  target  (cloud),  the  performance  of  a  me¬ 
teorological  radar  cun  be  derived  from  ihe  following  radar 
equation: 

10  log  >)  =  10  log  F,  -  10  log  F, 

-  10  log  |/i)  -  10  log  (/),) 

-F  10  log  (37/?’)  +  )J  ( 1 ) 
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where  n  is  the  rador  reflectivity  in  nei{oiive  reciprocal  cen* 
limeters,  P,  and  P„  respectively,  are  the  received  power 
and  the  transmitter  peak  power  In  watts.  A  Is  the  radial 
dimension  of  the  pulse  volume.  A,  is  the  antenna  effective 
area,  and  R  is  the  ranj^e  of  the  iar|et.  0  Is  a  corrective 
renn  that  includes  transmitter  antenna  efficiency  ( ~2  dB ) 
and  two'way  signal  attenuation  between  radar  and  target. 

If  the  cloud  panicles  are  so  small  that  even  the  shortest 
radar  wavelength  Is  associated  with  Rayleigh  scattering 
(cloud  panicle  diameter  smaller  than  approximately,  one- 
fourth  of  the  radar  wavelength),  the  radar  reflectivity  fac¬ 
tor.  Z  *  f  N(D)EPdD  (N  in  panicles  per  cubic  meter 
and  D  In  millimeters);  Is  given  by  Z  «  (qX^)/(r’{/r  |‘) 
where  {/r|'  is  a  scattering  parameter  discussed  below. 
With  the  more  familiar  JBZ  *  10  log  Z  units 

dBZ  =  10  log  q  +  95.14  +  40  log  X  -.10  log  |Ar|* 

(2) 

or 

dBZ  =  10  log  P,  +  40  log  X  -  10  log  P, 

-  10  log  A  -  10  log  Af  +  10  log  {2vR') 

+  10  log  t/r|' +  95.14.  (3) 

The  term  containing  X  show.s  the  drastic  inlliiencc  of 
the  radar  wavelength  on  r.idar  sensitivity,  which  is  a  strong 
incentive  for  the  choice  of  u  very  short  wavelength  for 
cloud  detection.  When  esamining  the  contribution  to  ra¬ 


dar  sensitivity  from  the  terms  In  (3).  It  should  be  remem¬ 
bered  that  echo  intensity  measurements  are  usually  per- 
fonned  using  a  signul  integrator  so  that  the  smallest  signol 
the  receiver  can  detect  Is  a  function  of  the  number  of  sam¬ 
ples  used  in  perfortuing  the  integration.  More  specifi¬ 
cally.  integration  of  a  radar  signal  P^,  which  is  performed 
on  Kceiver  signul  sumpics  digitized  at  the  radar  repetition 
rate  at  adjacent  range  gates,  has  two  effects.  Firstly,  it 
reduces  th£  standard  deviation  of  the  fluctuating  mean  es¬ 
timate  to  ?,/(A(,)''\  where  N,  is  the  equivalent  number 
of  independent  samples  (smaller  than  the  actual  number 
iVof.’iamples  depending  on  sample  to  sample  correlation). 
Secondly,  it  reduces  the  effective  receiver  noise  level  F, 
to  F,/(.V)''\  where  N  is  now  the  actual  number  of  sam¬ 
ples  integrated  us  they  are  statistically  independent  from 
pulse  to  pulse.  Therefore,  for  a  given  signal  dwell  time, 
the  radar  sensitivity  increases  not  only  with  the  pulse 
width  (or  pulse  volume  radial  dimension,  the  term  con¬ 
taining  A  in  the  equation  above),  but  also  with  tite  radar 
pulse  repetition  rate.  Radar  sensitivity  for  distributed  tar¬ 
gets  li.e..  clouds)  is  thus  a  function  of  the  transmitter  av¬ 
erage  power,  not  peak  power. 

We  are  now  equipped  to  e.xamine  the  relative  merit 
(compared  to  a  typical  5-bund  radar)  for  cloud  observa¬ 
tion  of  typical  radars  operating  at  wavelengths  ranging 
from  5-b:ind  (3  GHz)  to  the  present  technological  limit 
of  .tppro.iimatcly  220  GHz.  Table  I  is  based  on  average 
ch.iractcriNtics  of  meteorological  radars  acioally  in  use  and 
show,  ilic  wa\clcngih  depcmtcncc  of  the  iransndtter  av- 
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TABLI:  I 


Have 

len^Ui 

cm 

Transmitter 

Average 

power 

dfl(w) 

Antenna 

size 

dB(ai) 

Rec.  NT 
figure 
(dB) 

X-* 

gain 

dB 

Beam 

wldtli 

degree 

far 

field 

)(M 

Total 

dS 

10.0 

0(500) 

0(10,0) 

0(1) 

0 

0.7 

2,00 

0 

5.6 

-2(300) 

-10(3.6) 

-1(2) 

i-lO 

1.0 

0.23 

-3 

3.2 

-3(250) 

-10(3.0) 

-2(3) 

♦20 

0.7 

0.26 

♦5 

0.15 

-13(25) 

-14(2.0) 

-4(5) 

♦42 

0.3 

0.47 

♦11 

0.32 

-20(5) 

-20(1.0) 

-4(5) 

♦60 

0.2 

0.31 

♦16 

>0.14 

-30(0.5) 

-34(0.2) 

-7(6) 

•73 

0.3 

0.03 

♦2 

erage  power,  antenna  size,  receiver  noise  ievel,  1/X'* 
Rayieigh  scotiering  gain,  and  the  oveniii  merit  figure  of 
the  radar  for  cioud  detection.  The  radar  antennu  beam* 
width  and  the  far  field  factor  X  -  D‘/X  are  also  indi¬ 
cated. 

Table  1.  covering  radar  wovelengths  from  lO  to  0.14 
cm.  shows  that,  despite  large  variations  of  some  of  the 
individual  terms,  there  Is  only  a  relatively  smull  change 
of  the  overall  radar  sensitivity  for  cloud  observation.  Note 
that  the  obove  data  relates  to  typical  meteorological  radar 
characteristics  that  do  not  necessarily  represent  state-of- 
the-art  radar  development. 

Besides  the  far  field  condition,  which,  for  antennae  pro¬ 
ducing  the  same  beamwidth.  Is  more  favorable  at  shorter 
wavelengths,  the  most  Important  consideration  for  close 
range  detection  of  weak  atmospheric  targets  is  not  only 
the  receiver  noise  level,  but  also  ground  echoes  leaking 
through  the  antenna  sidelobes.  This  is  especially  true  for 
high-power  tong-wavelength  radars.  Assuming  the  same 
sidelobe  structure  for  all  antennae  regardless  of  the  wave¬ 
length.  the  effect  of  ground  echo  Interference  in  the  mea¬ 
surement  of  the  weak  backscattering  produced  by  a  small 
cloud  at  close  range  can  be  characterized  by  the  contrast 
between  returns  from  the  ground  and  from  clouds.  Ground 
echoes  are  caused  by  a  variety  of  large  size  targets  such 
as  buildings,  trees,  terrain,  etc.  The  variation  us  a  func¬ 
tion  of  radar  wavelength  of  ground  clutter  radar  cross  sec¬ 
tion  per  unit  area  o*  is  not  well  known  and  depends  on 
the  targets,  but  it  seems  to  increase  only  very  slowly  with 
a  decrease  of  the  wavelength.  On  the  other  hand,  cloud 
echo  intensity  increases  proportionally  to  I /\*.  or  a  dra¬ 
matic  Increase  of  60  dB  from  10-  to  0.32-cm  wavelength.- 
Therefore,  even  with  a  possible  lO-dB  increase  of  a*  with 
the  shoner  wavelength,  the  ground  cluitcr-cloud  return 
contrast  will  improve  by  at  least  50  dB  from  a  10-  to  0.3- 
cm  wavelength.  It  is  our  experience  that,  even  surrounded 
with  buildings,  a  94-GHz  radar,  providing  a  theoretical 
15-dB  (no-clulicr)  cloud  detection  improvement  over  an 
5-hand  radar,  does  not  show  any  measurahic  (le.vs  than  21) 
dB  hclow  receiver  noise)  ground  clutter  .signal  at  short 


ranges  (200  m)  for  elevation  angles  above  approximately 
5*.  Therefore,  a  mlHimeier-wave  radar  operated  in  a  ver¬ 
tically  pointing  mode  (or  slightly  off  venicol  If  three-di- 
menstonol  scanning  Is  required)  Is  an  attractive  solution 
for  the  observation  of  low-altitude  low-refiectivUy  clouds 
such  as  fair  weather  cumuli. 

The  94-GHz  radar  discussed  here  provides  a  Rayleigh 
scattering  gain  of  17  dB  over  a  35*GHz  rador  and  has  the 
following  main  characteristics;  peak  power  1.2  kW  (lO 
kW).  maximum  overage  power  5  W  (25  W),  antenna  size 
90  cm  (180  cm),  receiver  double  side-band  noise  figure 
6  dB  (S  dB).  The  numbers  quoted  within  parentheses  are 
feasible  Improvements  of  the  radar  eharacteristies  that  are 
not  implemented  oi  the  present  time.  The  radar  described 
in  (1)  operates  at  35  GHz  and  has  the  following  cliarae* 
teristies:  average  power  25  W.  antennu  size  1 80  cm.  A 
receiver  noise  power  of  -85  dBm  is  quoted  by  Hobbs  but 
with  the  indicated  receiver  bandwidth  and  an  up-to-date 
mixer,  the  receiver  noise  should  be  down  to  -95  dBm 
(6*dB  DSB  noise  figure). 

Assuming  that  such  a  receiver  noise  Improvement  is  re¬ 
alized  in  the  Hobbs'  radar,  the  94-GHz  radar  exhibits  a 
dB  sensitivity  gain  on  cloud  detection  with  a  3-fi  dish, 
and  1 1  dB  with  a  b-ft  dish.  With  an  average  power  in¬ 
crease  to  25  W.  the  sensitivity  gain  is  17  dB^  Small  fair 
weather  cumuli  are  at  the  threshold  of  detection  (dBZ  - 
-50)  and  even  a  relatively  small  improvement  of  the  ra¬ 
dar  sensitivity  may  make  a  signiHcant  difference  in  their 
detectability. 

The  94-GHz  radar  has  a  much  smaller  size,  is  operated 
at  low  power  (a  few  hundred  wait.s  of  power  .supply)  and 
doe.s  not  require  waveguide  pressurization.  The  94-GH2 
radar  also  use.s  a  gridded  extended  interaction  oscillator 
(ElO).  which  e.\hibit.s  much  less  inirapiilse  FM  or  phase 
noi>e  j6)  than  the  magnetron  u.sed  in  the  35-GHz  radar. 
Thi>  i>  confirmed  hy  oh.\crvaiions  of  fixed  target  signal 
pha>e  coherence  and  also  tests  of  mean  Doppler  stability 
and  -rcctral  variance  on  cloud  echoes  [3)  that  indicate 
lh.ii  the  radar  pli.i.NC  nni.se  eomrihiuion  to  the  measurement 
of  Dcri'lcr  .-ipeciriim  nmmciiis  is  negligible. 
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n|. AlmoiptMric  tncmniion  M  94  OHs  u  i  ftofietioa  of  tpoelflc  hu* 
mtdiiy  iR  gnmi  diMt  cubic  mcitii,  u  iroond  prction  ind  for  two  itm* 
pcnnict,  310  tnd  3i0*K. 

in.  AuoRrnoN  of  94<OHz  Radution  by 

ATMOSFHBRtC  GaSES  AND  HyDXOMBTEORS 
The  eietr  air  retidual  ionospheric  absorption  in  the  94* 
GHz  window  Is  due  to  the  skirts  of  the  Oj  and  HjO  lines 
and  depends  on  altitude  (pressure  broadening)  and  spe* 
clfie  humidity.  There  is  a  small  variation  with  tempera* 
ture.  The  absorption  varies  from  approximately  0.4 
dB/km  for  a  low  0.25  g/m*  humidiQr  to  approximately  2 
dB/km  for  a  high  25  g/m’  humidity  (see  Fig.  2). 

Absorption  is  a  line  Integral  of  the  absorption  eoefll* 
cients  evaluated  from  the  radar  to  the  target.  Since.,  in 
clear  air,  water  vapor  is  the  dominant  absorption  agent, 
the  maximum  zenith  absorption  (l.e.,  to  the  top  of  the 
atmospdtere)  for  a  vertically  pointing  beam  Is  a  function 
of  humidity  conditions  and  can  be  derived  from  the  ver¬ 
tically  integrated  water  vapor  that  Is  often  chancterized 
as  the  precipitable  water.  The  94-CHz  absorption  coeffi¬ 
cients  based  on  data  by  Llebe  (8]  and  Shimabukuro  and 
Epstein  (15)  are  shown  in  Fig.  3.  The  one-way  zenith  ab¬ 
sorption  varies  from  approximately  0.5  dB  for  a  dry  atmo¬ 
sphere  (0.6*cm  precipitable  water)  to  more  than  3  dB  for 
a  very  humid  one  (5-cm  precipitable  water).  It  increases 
by  3  percent  for  a  radar  beam  IS*  off  vertical.  15  percent 
for  30*.  and  doubles  for  60*.  Although  the  vertically 
pointing  beam  mode  yields  the  least  absorption,  off  ver¬ 
tical  scanning  up  to  30*  offering  a  beam  scanning  capa* 
bility  can  be  tolerated  even  in  the  case  of  a  humid  atmo¬ 
sphere. 

Attenuation  of  94-GHz  radiation  by  clouds  and  precip¬ 
itation  is  caused  by:  1)  absorption  of  the  radiative  energy 
being  transferred  into  heating  of  the  panicles:  21  scatter¬ 
ing  of  radiative  energy  by  the  hydro  meteors  <  The  absorp¬ 
tion  cross  section  Q„  and  the  scattering  cross  section  Qi 
of  a  sphere  (cloud  droplet)  much  smaller  than  the  radia¬ 
tion  wavelength.  (Rayleigh  absoqition)  is  given  by 

G,  =  tr'DMm(-A)A  (3b) 

G  -  3r'D'’|/:|7.-tV  (4) 


FfteemTABU  WATER  (CM) 

Ftf.  ).  Ont* way  nnhh  tbMqKien  u  a  (tiiKthm  of ptteipittblt  wiitr.  TIm 
leltd  iiM  ii  Iht  modtl  propoMd  by  Sklimbiilivre  u  t  katt  iqgim  III  of 
liij  diu.  TIm  eireki  in  tht  valun  propowd  by  Lkb*  e«  Uit  bull  ofkU 
oxptdiMRU. 


where  K  *  (m*  —  I)/(m*  +  2),  and  m  ■  /i'  +  in"  is 
the  complex  index  of  refraction  of  water  (or  Ice).  Im 
( -/T)  is  the  imaginary  part  of  K.  The  attenuation  cross 
section  it  Q,  •  Qt  +  Q,.  For  a  cloud  of  particles,  the 
absorption,  scattering,  and  extinction  coefficients  (in  dec¬ 
ibels  per  kilometer),  a,,  er,,  ot„  respectively,  are  ex¬ 
pressed  by 

a,  -  0.434  £  a  (d) 

(similariy  for  a„  a,  with  the  summation  extending  to  all 
particles  present  in  a  cubic  meter  of  air).  In  the  ease  of 
Rayleigh  absorption  (droplets  much  smaller  than  the  radar 
wavelength)  in  decibels  per  kilometer  Is  given  by:  oi« 
*  8.lS<Vf  Im  ( -K)/\  where  M  is  the  cloud  liquid  water 
in  grams  per  cubic  meter. 

At  94  GHz,  the  complex  index  of  refraction  of  liquid 
water  must  be  intetpolated  between  actual  measuremenu 
at  X  *  1  mm  and  X  «■  5  mm,  using  the  Debye  equations 
(1 1]  and  is  seen  to  vary  with  temperature.  Table  II  shows, 
at  94  GHz,  the  real  and  Imaginary  parts  of  m,  the  absorp¬ 
tion  index  Imag  ( -K ),  the  absorption  coefficienu.  and 
the  scattering  index  |  ^ 

The  o«  values  above  are  presented  graphically  in  Fig. 
4  where  a«  is  seen  to  vary  almost  lineady  with  tempera¬ 
ture.  The  scattering  cross  section  Q,  is  negligible  for  Ray¬ 
leigh  absorbers. 

In  the  case  of  larger  size  particles  for  which  Raylei^ 
scattering  is  not  applicable,  the  computation  of  a  for  a 
given  cloud  or  precipitation  must  be  based  on  a  knowl¬ 
edge  of  panicle  size  distribution,.  To  overcome  this  re¬ 
quirement  and  still  provide  an  insight  of  the  imponance 
of  94-GHz  signal  absorption  by  precipitation,  the  cuntri- 
bution  to  a  due  to  1  g/m^  of  water  divided  into  drops  of 
the  same  diameter  D  is  calculated  as  a  function  of  that 
diameter. 

The  ctj  and  a,  coefficients  computed  in  these  conditions 
using  .\Jie  tables  arc  shown  in  Fig,  5  as  a  function  of  D. 
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TABLE  II 

iNon  or  RinACTtON  or  Uouto  Water  AasoaTJO-N  axo  SCArrcaixa  Tkims  at  94  GHi 


TMporatura 

n' 

n* 

X»(-X) 

|K|‘ 

0‘C 

2.84 

1.48 

0.1938 

4.80 

0.711 

IS'C 

3.21 

1.79 

0.1S14 

4.0 

0.787 

20‘C 

3.41 

2.02 

0.1470 

3.64 

0.828 

23 ’C 

3. OS 

2. OS 

0.1402 

3.47 

0.832 

for  a  liquid  water  content  of  I  g/m\  a„.  which  exhibits  leigh  region  for  u  panicle  size  greater  than  3  mm.  The 

a  constant  valoe  of  3.47  dB/km  (temperatore  23*C)  in  contribution  to  a  a„  which  is  negligible  in  the  Rayleigh 

the  Rayleigh  region ^  increases  to  a  nta.ximom  of  9  dB/km  region,  increases  to  a  maximum  of  approximately  8 

for  a  panicle  diameter  of  approximaiely  I  mm.  and  de-  dB/km  for  D  =  I  mm.  and  decreases  to  approximately 

crea.se.s  to  a  value  smaller  than  that  observed  in  the  Ray-  3.5  JB.'km  for  D  =  3.5  mm. 
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In  the  case  of  ice  particles  we  have:  m  «  1.878  +  <4.76 
10-*  (0*0  $0  that  Imf-AT)  ~  2  lO'^'ando,  ~  6.10"* 
dB/km.  Although  it  may  be  concluded  that,  even  at  94 
GHz,  signal  attenuation  by  an  ice  cloud  is  negligible,  this 
is  only  true  for  Rayleigh  scattering,  and  perhaps  for  Mie 
scattering  if  we  consider  only  attenuation  caused  by  ab> 
sorption.  At  94  GHz,  60-70  Ice  spheres  produce  a 
scattering  contribution  to  attenuation  equal  to  that  due  to 
absorption.  For  larger  sizes,  the  scattering  contribution 
becomes  rapidly  overwhelming,  reaching  mote  than  10 
dB/km  for  particle  sizes  of  D  ■■  1-2  mm.  Note  that,  ex* 
cept  for  solid  hailstones,  the  assumption  of  spherical  par* 
tides  is  not  realistic  and  should  be  replaced  by  the  com* 
plex  structure  of  sn  ice  crystal  or  snowflake  for  which  the 
andytica!  Mle  solution  Is  inapplicable, 

IV,  Cloud  and  PMawTATiON  Radar  REaEcnvmf 
AT  94  GHz 

In  the  Rayleigh  region,  the  backscattering  cross  section 
of  s  sphere  Is  c  irD^[K  IV^**  w  thst  cloud  radar  re¬ 
flectivity  1}  Is 

i,-(ir‘|A:lVx')  jN(D)dD.  (6) 

q  Is  In  recinrocal  cemlmeten  and  Z>  and  X  are  In  cen¬ 
timeters.  ( iC  p  is  shown  In  Table  11  and  Is  seen  to  Increase 
only  slightly  with  temperature. 

Table  III  shows,  at  94  OHz  and  23*  temperature,  the 
radar  reflectivity  of  I  g/m*  liquid  wster,  divided  In  drop* 
lets  of  the  same  diameter  D,  as  a  function  of  the  droplet 
diameter  D. 

For  non*RayleI|h  scatterers  (drops  having  a  diameter 
larger  than  0.5  mm  at  94  GHz),  the  computation  of  radar 
reflectivity  must  be  performed  using  the  Mle  backscatter¬ 
ing  cross  section  that  is  shown  In  Fig.  6.  The  radar  re¬ 
flectivity  is  given  by  an  Integral  over  all  particle  sizes, 
thereby  Involving  the  drop  size  distribution.  It  was  shown 
(3]  that,  usumlag  a  Marshall-Palmer  11948]  drop  size 
distribution,  the  maximum  radar  reflectivity  encountered 
In  heavy  rain  will  not  exceed  30  dBZ.  However,  in  order 
to  provide  an  Insight  Into  the  Influence  of  particle  size  on 
ra^r  reflectivity  at  94  GHz,  the  contribution  to  radar  re¬ 
flectivity  arising  from  1  g/m*  of  water  divided  into  drop¬ 
lets  of  a  given  diameter  D  is  shown  In  Fig.  7  as  a  function 
of  D  for  both  water  and  ice.  A  iiuxlmum  reflectivity  Is 
reached  if  the  same  liquid  water  is  divided  into  1-mm 
drops.  Larger  sizes  are  associated  with  less  reflectivity; 
however,  there  are  deep  oscillations  of  the  function  in  the 
precipitation  particle  size  range.  Comparing  the  scattering 
functions  for  ice  and  water  indicates  that,  for  Rayleigh 
scattering,  there  is  an  increase  of  reflectivity  of  approxi- 
maicly  6  dB  if  a  panicle  changes  from  ice  to  water  during 
melting.  However,  this  elTect  is  reversed  for  panicle.s  of 
larger  size  for  which  Mie  scattering  applies.  Indeed,  at  94 
GHz,  backscaiiering  by  ice  spheres  having  a  diameter 
greater  than  2  mm  is  much  stronger  than  thai  produced  by 
water  spheres  of  the  same  size. 

The  presence  of  deep  Mie  scattering  o.scillaiions  within 


TABLE  Itt 


D 

dBZ 

1  U 

4.6  10’^* 

-57 

10  p 

4.6  10'“ 

-27 

100  p 

4.6  10'* 

+03 

the  raindrop  size  range  may  provide  an  opportunity  to  re¬ 
motely  identify  raindrop  size..  This  is  realized,  for  in¬ 
stance.  if  observation  of  the  Doppler  spectrum  (density 
distribution  of  reflectivity  driJdV  as  a  function  of  radial 
speed  k")  is  made  in  rain  conditions  at  vertical  incidence. 
To  illustrate  this  method.  Fig.  8  shows  a  theoretical  eval¬ 
uation  of  a  Doppler  spectrum  observed  at  vertical  inci¬ 
dence  in  u  I  mm/h  rainfall  a.ssuming  a  Mar.shall-Palmer 
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Fli.  7.  Ridsr  nJiMlIvliy  MmritaihM  bu«d  on  Ml<  miMrins  fuiKiluu 
by  yinkiti  Af  liUmcnr  0.  p«r  fnmt  llmt  cubic  mcicn  uf  wticr  con> 
wni.  u  I  fbiKihM  a(  0.  The  ««m  ctlcul«ii«M  for  kc  tpberex  era  cImi 
InOicued. 


FIs.  1  Doppler  spcciniffl  cl  v*nlc»l  UcUcnct  for  i  (•inin/h  nlflfcll  i<- 
Mmlns  «  MtrthalUFalmcr  drop  tin  dlHribuilon.  cakulalcd  w  94>  end 
9.}*OHi  radtr  fit^Muty. 

drop  size  distribution.  The  oscillations  of  the  Mte  scat* 
tering  functions  are  cieariy  visible.  The  Doppler  spectrum 
observed  in  the  same  rainfall  but  at  9  CHz  is  also  shown 
for  comparison.  Fig.  9  shows  the  ratio  r  between  the  two 
(94  and  9  GHz)  spectra.  This  ratio  r  is  independent  of 
the  drop  size  distribuiion  and  only  a  function  of  drop  size 
and.  if  air  vertical  velocity  is  neglected,  a  function  of  the 
raindrop  vertical  velocities  or  Doppler  frequency  shifts  at 
vertical  incidence.  One  sees  that  in  the  Rayleigh  region. 
T  is  40  dB.  which  reflects  the  X'**  relationship,  decreases 
to  IS  dB  for  0  »  I  mm.  and  then  o.scillates  with  well 
delined  ma.xima  at  D  =  2.2.  3.5.  and  4.75  miii<  r  is  also 
shown  for  ice  and  indicate.'!  that  the  rolloff  of  the  function 


:t.» 


Fl|.  9.  lUtIo  between  Doppler  ipeetn  it  9t*  end  9.]-GH>  wevtti«stl>  ** 
e  SjACtioA  of  pnitlele  dlemeier.  Tbie  utio  it  effecilvely  repretemed  by 
dN  laiio  between  ndsr  erect  teettone  n  the  two  wivtlensihi.  The  nculu 
ere  tito  shown  for  Ice  tphctti. 

with  increasing  panicie  diameter  is  iimited  to  lets  than  15 
dB  instead  of  40  dB  in  the  water  case.  This  may  offer  a 
means  to  differentiate  between  raindrops  and  hailstones. 

V.  StONAL  PaOCESSOK  ANO  RaDAR  PlRFOR.VtANCB 

The  radar  is  equipped  with  anauiocovariance  processor 
(pulse  puir  processor)  and  a  signal  integrator  thut  were 
both  described  earlier  (2).  The  two  terms  of  the  autoco* 
variunce  function  and  the  signal  complex  square  (0  lug 
correlation  term)  as  weii  as  the  integration  of  a  iog  re* 
ceiver  signal  are  evaluated  in  real  time  with  a  fast  digital 
hardware  processing  unit  and  recorded  on  magnetic  tupe. 
The  computation  of  mean  Doppier  and  spectrai  width, 
which  is  performed  by  a  microcomputer  from  the  re* 
corded  duia.  was  discussed  previously  (2).  'Hits  method 
allows  fle.xibility  for  the  computation  of  spectni  nin* 
meats,  particularly  useful  for  the  calibration  of  the  spec* 
trum  width  for  which  receiver  noise  has  to  be  estimated 
|7].  The  processor  produces  pulse  pair  data  continuously 
evaluated  at  4(X)*ns  range  intervals.  Most  of  the  data  re* 
corded  so  far  were  obtained  with  a  lO-kHz  PRF  ( ig  nt/s 
unambiguous  Doppler  velocity  range),  a  0.4*ps  pulse 
width  providing  a  range  resolution  of  60  m.  and  32  000 
samples  (3.2'S  dwell  time)  used  in  the  evuluation  of  the 
autocovariance  and  signal  intensity  terms.  In  these  enn* 
ditions.  mean  Doppler,  spectral  width,  and  mean  signal 
intensity  can  be  extracted  from  signals  22  dB  below  the 
receiver  noise  ievci  l3J. 

The  radar  has  been  operated  with  a  transmitter  peak 
power  of  1.2  kw  (pulse  width  0.4  ^s),  and  a  90-ciu  di* 
ameter  antenna  with  a  one-way  efliciency  factor  estimated 
to  be  2  dB.  In  the.se  conditions,  the  minimum  radar  re¬ 
flectivity  that  can  be  detected  by  the  radar  at  a  I -km  ranee 
is  given  by  inserting  these  values  in  ( l)^  We  have 

10  log  1}  ^  P,  (dBm)  -18  +  20  log  R  17) 

or  using  the  more  familiar  dBZ  with  lO  log  tf  =  dBZ  - 
76.2  ai'94  GHz.  we  have 

dBZ  =  P,  (dBm)  +  .57.8  +  20  log  R.  1.81 
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The  above  values  only  apply  to  Rayleigh  scattering.  For 
larger  panicles,  we  must  define  an  equivalent  Z,.  i.e..  the 
radar  reflectiviiy  factor  of  Rayleigh-scattcring  raindrops 
required  to  produce  an  echo  intensity  equal  to  that  ob¬ 
served. 

In  the  radar  receiver,  the  signal  et  the  output  of  a  log¬ 
arithmic  amplifier  is  integrated  using  32  000  samples  so 
that  its  level  Is  determined  with  epproximately  O.Ol-dB 
precision.  The  method  for  extracting  the  signal  power  Pg 
is  based  on  measuring  the  ratio,  a  •  (Pt  +  Ph)/Pm  {P» 
is  the  noise  power),  and  then  computing  *■  10*  - 
1.  Assuming  that  the  noise  at  the  output  of  the  receiver  Is 
entirely  due  to  the  94-GHz  receiver  Input  noise  and  esti¬ 
mating  this  noise  from  receiver  eallbradon  using  solar  ra¬ 
diation  received  by  the  antenna  [3]  provides  a  technique 
for  the  measurement  of  signal  power.  Using  32  000  sam¬ 
ples  In  the  integration  yields  a  theoretical  noise  reduction 
of  22  dB,  which  is  confirmed  by  our  testt.  The  receiver 
noise  power  is  estimated  to  be  *’93  dSm  ( 15-MHz  re¬ 
ceiver  bandwidth,  O-dB  DSB  noise  figure)  so  that  the 
minimum  detectable  slgnel  In  these  conditions  is  about 
-US  dBm.  If  we  report  these  data  in  (8),  the  minimum 
dBZ  velue  that  the  rader  cen  detect  at  a  2-km  range  (typ- 
Icel  boundeiy  layer  cloud  top)  Is  -  -SO  dBZ.  Table  Ill 
indicates  that  -SO  dBZ  is  the  radar  reflectivity  of  a  cloud 
with  10  mg  of  liquid  weter  and  a  particle  size  of  10  pm, 
which  Is  less  than  the  typical  conditions  encountered  In  a 
felr  weather  cumulus. 

VI.  Moot  or  OniiATtoN  and  Risults  or  Ctx>uD 
OiseavATtON 

One  of  the  project  objectives  was  to  observe  the  vertleel 
air  velocity  iv.  which.  In  the  cese  of  vei>  small  droplets 
having  negligible  terminal  velocity.  Is  directly  measured 
using  a  fixed  vertically  pointing  radar  beam.  This  mode 
of  operation  yields  e  venical  cross  section  of  the  observed 
clouds  where  time  can  be  replaced  by  a  space  coordinate 
if  It  Is  essumed  thet  the  clouds  are  at  a  nearly  steady  state 
end  translating  at  a  fixed  speed.  This  Is  ecceptable  for 
stratiform  clouds  or  precipitation  conditions.  In  the  case 
of  an  isolated  cloud  such  as  e  fair  weather  cumulus,  the 
method  can  still  yield  very  useful  data  but  cen  be  advan¬ 
tageously  replaced  by  a  two-dimensional  cloud  scenning 
methodology.  The  radar  beam  is  now  scanned  slightly  off 
venical  (i.e.,  en  elevation  angle  of  90  ±  13*  that  does 
not  appreciably  change  the  signal  attenuation  conditions) 
and  that,  combined  with  the  cloud  motion,  yields  en  ef¬ 
fective  three-dimensional  probing  of  the  cloud.  The  only 
problem  is  the  contamination  of  vertical  velocity  by  hor¬ 
izontal  velocity,  which  occurs  for  nonvenical  beam  an¬ 
gles.  However,  if  the  technique  is  applied  to  the  obser¬ 
vation  of  a  fair  weather  cumulus  moving  in  a  well-defined 
direction,  the  .scanning  can  be  done  in  a  direction  perpen¬ 
dicular  to  the  cloud  motion  .so  that  any  contribution  to 
radial  velocity  hy  the  cloud  horizontal  velocity  is  mini¬ 
mized. 

.Such  three-diiiicn.sional  .scanning  i.s  now  planned  but  all 
the  ‘)4  GH/  rad.ir  ohsorsations  of  clouds  in  the  past  were 


performed  using  a  fixed  vertically  pointing  radar.  This 
methodology  is  appropriate  for  the  observation  of  cloud 
layers  such  es  stratus,  stratocumulus,  cirrus,  and  also 
stratiform  precipitation.  Incidentally,  the  experience  ec- 
quired  using  the  radar  indicates  thet  data  in  stratified, 
slightly  precipitating,  nimbostratus  clouds  are  obtained  up 
to  the  cloud  top  at  8  km.  High-eltitude  ( 12-14  km)  cirrus 
are  always  detected  even  if  barely  visible  with  the  naked 
eye. 

In  the  case  of  Isolated  clouds  such  as  fair  weather  cu¬ 
muli.  there  la  only  a  small  probability  for  the  cloud  to  pcss 
overhead  with  the  radar  Intersecting  the  cloud  center  so 
that  1  representative  vertical  cross  section  of  both  radar 
reflectivity  and  vertical  velocity  could  be  obtained.  Dur¬ 
ing  a  two-month  cloud  observation  campaign  In  the  1986- 
1987  winter  In  Florida,  we  were  fortunate  to  And  that,  at 
eeittin  occasions,  the  radar  happened  to  be  underneath  a 
line  of  cumuli  that  were  systematically  moving  with  the 
environmental  wind  and  that  some  of  these  clouds  were 
passing  directly  overhead  with  the  radar  beam  going 
through  the  cloud  center.  These  were  conditions  favorable 
for  acquisition  of  vertical  velocity  and  radar  reflectivity 
data  In  a  cloud  region  representative  of  the  cloud  circu¬ 
lation.  To  confirm  this  favorable  aspect  of  the  observa¬ 
tions  the  visual  Imege  of  the  clouds  passing  directly  over¬ 
head  was  recorded  by  e  video  camera. 

An  example  of  the  results  Is  shown  in  Fig.  lOfal  end 
(b).  The  date  presented  are  radar  reflectivity  end  vertical 
velocity  In  a  venlcel  cross  section  defined  by  the  trans¬ 
lation  of  the  cloud  through  the  radar  beam.  These  data 
were  presented  before  [4],  but  they  have  been  reaneiyzed 
to  show  more  detell  In  both  the  reflectivity  and  the  radial 
velocity  fields.  The  radar  reflectivity  field  Is  shown  In  Fig. 
lOfa).  The  meximum  radar  reflectivity  Is  -31  dBZ  end 
the  minimum  detectable  reflectivity  is  -50  dBZ.  There  is 
e  sharply  defined  cloud  top  due  to  a  temperature  inversion 
clearly  Indicated  by  the  radiosonde  data  collected  thet  day 
[4].  The  cloud  base  is  not  es  well  defined  and  exhibits  a 
noticeable  variability.  The  venical  velocity  Is  sho«n  in 
Fig.  10(b)  and  reveals  the  presence  of  a  strong  updraft 
region  (updraft  maximum  velocity  3  m/s)  In  the  cloud 
center  and  two  significant  downdrafts  below  and  toward 
the  cloud's  leeding  edge.  Comparing  the  velocity  and  the 
radar  reflectivity  fields  indicates  that  the  bulk  of  the  up¬ 
draft  is  just  below  (and  slightly  lagging  in  lime)  a  region 
of  high  radar  reflectivity,  although  the  downdrafts  are  al¬ 
ways  in  a  region  of  very  low  radar  reflectivity.  The  pres¬ 
ence  of  the  updrah  is  associated  with  local  rising  of  (he 
cloud  top  (still  strongly  capped  by  the  temperature  inver¬ 
sion  lid),  and  the  downdrafts  are  a.ssociaied  with  a  local 
lowering  of  the  cloud  bottom  (hat  can  considerably  disiun 
(he  cloud  base,. 

Fig.  1 1  shows  another  typo  of  data  collected  by  the  mil- 
liraeier-wavc  radar.  These  oh.sers’aiion.s  were  made 
through  the  melting  region  of  .a  stratiform  cloud  previ- 
ousls  named  ’’bright  bund”  because  of  a  systematic  mas- 
inuim  of  echo  intensity  obsened  just  helow  the  0®  Iso^ 
ihemi.  The  previous  bright  hand  obvervaiioiis  with  cen- 
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Fig.  10.  (<)  Januiry  34.  1917  iimf«)Mi|hi  preianuiwn  o(  rtdtr  rtiltcitvUy  Ia  dBZ  oAitrvtd  Ia  a  small  aumutui  ptsslAf  evtrhiad.  Tha  alcud  Wf)  Is 
siroAfly  llmiiad  by  a  cappiAi  lAvaraioA.  Tima  is  EOT.  Ths  cloud  is  movlAp  at  !0  mr'  so  (hut  I  min  is  squivaisAi  to  <00  m.  lb)  Sama  as  (a)  but  for 
iba  vanical  air  valocliy.  Tha  dashad  araa  iAdieaias  updraft  rt|iOAS.  Updrafts  tard  m  pus.h  iha  cloud  top  aad  dowAdrahs  to  lower  the  cloud  base. 
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(a)  <b)  (cl 

Fif.  II,  Venteal  profiles  of  radar  rtfieciivily  aAd  mean  Doppler  observed  at  94  GH:  ihrouyh  the  mcliiAg  baad  Ia  stniiform  laiA  at  difiercAt  localioAS 
aAd  lime  |a)  July  2.  1933.  Miami.  FL:  lb)  February  6.  19S7.  Miami.  FL;  (el  .May  t3.  (937,  Boston.  ,MA. 

limeier-wave  radars  were  characterized  by  a  well  deSiied  94  GH:.  of  vertical  profiles  of  both  radar  reflectivity  and 
refieciivity  maximum  and  laier,  when  Doppler  radars  be-  mean  Doppler  velocity  are  presented  in  Fig.  1 1 ,  Although 
came  available,  by  a  sudden  increase  of  ihc  particle  ver-  these  obsen'aiions  were  acquired  at  differeni  limes  and 
tical  velocity  due  lo  meliing  (Sj.  locations,  they  exhibit  the  same  general  characteristics. 

Several  case.s  of  .such  siraliform  cloud  observations,  at  Firvth.  there  is  the  same  sudden  change  of  the  panicle 
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venicyl  velocity— attributed  to  the  oivict  of  panicle  melt¬ 
ing  wheo  they  fall  through  the  melting  zooc— which  was 
previous  observed  by  centimeter-wave  radars.  There  is 
also,  on  the  top  of  the  melting  zone,  the  same  6-7  dB 
increuse  of  radar  rcHectivity  that  was  observed  by  centi¬ 
meter-wave  radars  and  is  attributed  to  the  change  of  index 
of  refraction  from  ice  to  water.  However,  the  94-GHz  ob¬ 
servations  show  that,  although  the  particles'  acceleration 
in  the  melting  zone  Is  comparable  to  that  observed  with 
centimeter-wave  radan.  it  is  no  longer  associated  with  a 
decrease  of  the  reflectivity  observed  when  descending 
through  the  melting  zone.  The  decrease  of  radar  reflectiv¬ 
ity  below  the  maximum  clearly  observed  with  centimete^ 
wave  radars  was  primarily  attributed  to  the  decrease  of 
particle  concentration— and  thus  radar  reflectivity— due  to 
the  systematic  acceleration  of  the  particles  through  the 
melting  zone  end.  to  a  lesser  degree,  to  the  shrinking  of 
panicle  size  due  to  melting.  Since  that  reflectivity  de¬ 
crease  is  not  observed  by  the  94-GHz  radar  although  the 
punicie  acceleration  through  the  melting  zone  Is  siiil  ob¬ 
served.  we  must  conclude  that,  when  panicles  shrink  dur¬ 
ing  their  melting— and  thus  accelerate— their  radar  cross 
section  (and  therefore  the  radar  reflectivity)  increases.  The 
answer  to  this  apparent  anomaly  may  be  found  in  Fig.  6, 
which  shows  that,  due  to  the  peculiarities  of  Mie  Katte^ 
Ing.  the  radar  cross  section  of  a  1 .5-mm  panicle  witt  in¬ 
crease  by  more  than  10  dB  if  it  shrinks  to  a  I -mm  rain¬ 
drop.  The  minimum  of  radar  reflectivity  above  the  melting 
band  observed  with  the  mmimete^wave  radar,  which  is 
someiimes  associated  with  an  increase  of  the  panicles* 
vcnical  velocity,  may  also  be  attributed  to  the  same  .Mie 
scattering  behavior.  The  duta  presented  in  Fig,  1 1  were 
all  acquired  in  very  light  precipitation  conditions  ( truce 
to  less  than  I  mm/h)  with  a  mean  Doppler  veltK'ity  below 
the  melting  region  reaching  approximately  4  m/s.  which 
is  the  terminal  velocity  of  a  I -mm  raindrop. 

VII.  Conclusion 

The  94-GHz  Doppler  described  here  is  about  the  shon- 
est  wavelength  radar  that  can  be  assembled  for  meteoro¬ 
logical  use.  Radar  components  are  available  off-the-shelf. 
The  radar  is  very  small,  can  be  easily  tran.sponed.  and  has 
demonstrated  e.xcellent  Doppler  perhrrmance.  !t.s  capabil¬ 
ities  for  cloud  detection  and  its  ponabiliiy  make  it  a  good 
cloud  physics  tool,  useful  for  the  study  of  boundary  layer 
cumuli  and  clouds  in  the  initial  stages  of  their  develop¬ 
ment. 

Because  of  the  deep  Mie  backscuitering  oscillations  oc¬ 
curring  in  the  raindrop  panicle  size  range,  the  94-GHz 
radar  is  also  an  attractive  choice  for  drop  size  di.siribuiion 
mea.surements  based  on  a  muliiwavclengih  methodology, 
with  the  94-GHz  radar  observations  done  jointly  with  un 
S-baod  or  an  X-band  radar. 
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ABSTRACT 

A  Doppler  radar  operating  at  3.2  mm  wavelength  was  designed  and  assembled  primarily  for  observation  of 
clouds  and  precipitation.  Phase  detection  of  the  radar  signals  which  is  required  for  Doppler  operation  is  imple¬ 
mented  through  the  use  of  a  coherent  oicillator  phase  locked  on  the  transmitter  pulse  and  used  as  a  reference 
in  the  phase  detector.  The  ndar  and  associated  signal  processing  techniques  tueh  as  signal  Integrator  and  signal 
aulocovatiavxe  estimator  are  discussed  along  with  the  Doppler  performance  of  the  radar.  Also  presented  are 
the  results  of  observation  of  ice  and  water  clouds  and  also  precipitation,  which  show  the  excellent  Doppler 
capabilities  of  the  radar  in  terms  of  accuracy  of  the  mean  Doppler  and  Dopi^er  spectrum  width. 


1.  Introdoctloa 

Observation  of  motion  and  growth  of  small  hydro* 
meteon,  such  as  cloud  droplets,  is  a  powerful  method 
for  the  study  of  the  physics  and  dynamics  of  clouds  in 
eaily  stages  of  their  development.  A  Doppler  radar  can 
be  used  for  this  purpose  as  it  can  monitor  particle 
growth  by  measuring  radar  reflectivity,  as  well  as  si* 
multaneously  observing  updraft  and  entrainment  aris* 
ing  from  latent  heat  relea^  by  condensation  and/or 
firing  processes. 

Centimeter  wavelength  .odan  may  have  the  sensi¬ 
tivity  required  to  detect  small  cloud  droplets  if  these 
radars  are  equipped  with  a  high  power  transmitter  and 
a  large  iqperture  antenna.  However,  this  approach  in¬ 
volves  bulky  and  heavy  equipi  lent  while  short  range 
performance  of  such  radars  is  a’tually  limited  by  ground 
clutter.  Since  the  backscattering  fix>m  small  tar^  such 
as  cloud  particles  increases  proportionally  to  a  rel¬ 
atively  low  power  radar  operating  at  a  shorter  wave¬ 
length  may  constitute  a  better  alternative.  Signal  ab¬ 
sorption  due  to  atmospheric  gases,  clouds  and  precip¬ 
itation  may  be  prohibitive  at  certain  wavelengths. 
Figure  1,  derived  from  Liebe's  (1985)  data,  shows  at- 
mosf^eric  attenuation  for  centimeter  and  millimeter 
wavdengths.  At  low  humidity,  spectral  windows  suit¬ 
able  ibr  radar  operation  are  cleariy  identified.  The  fre¬ 
quency  bands  are:  K,-band  (35  GHz,  A  *  0.85  mm), 
W-band  (94  GHz,  A  =  3.2  mm),  F-band  (140  GHz,  A 
=  2.14  mm)  and  G-band  (220  GHz,  A  =  1.36  ram). 
Note  that  in  high  humidity  conditions,  these  “win¬ 
dows"  are  almost  filled. 

The  shortest  wavelength  used  in  previous  meteo¬ 
rological  radars  was  A  =  0.85  cm  (K«-band).  W-band 
radar  components  are  now  available  at  an  affordable 
cost  and,  since  the  A"^  improvement  at  94  GHz  is  40 
dB  with  respect  to  X-band  and  17  dB  with  respect  to 


K«-band,  a  W-band  Doppler  radar  can  be  considered 
for  cloud  observations.  The  W-band  radar  has  the  ad¬ 
vantage  of  low  power,  lightness,  smlll  overall  size  and 
mobility.  Thus,  it  can  be  easily  transported  in  a  small 
vehicle  or  installed  aboard  an  airplane  or  an  orbiting 
satellite.  When  compared  to  that  of  a  longer  wavelength 
radar,  the  W-band  radar  performance  at  short  range  is 
less  likely  to  be  degraded  by  ground  clutter,  thereby 
allowing  observation  of  targets  in  close  proximity  of 
the  radar  (less  than  100  meters). 

A  94-GHz  meteondogica]  I^ppler  radar  has  been 
proposed  (Lhermitte,  1981;  Lhermitte  and  Frusb, 
198S).  The  radar  has  now  been  constructed  and  op¬ 
erated  for  ctoud  and  storm  observations;  this  paper 
presents  the  radar  design,  an  evaluation  of  the  radar 
performance,  and  some  preliminary  results. 

2.  Pulse  Doppler  Radar  Design 

a.  Fulse  Doppler  techniques 

Pulse  Doppler  radars  measure  the  pulse-to-pulse 
echo  phase  change  due  to  target  motion  and  convert 
it  into  a  Doppler  frequency  shift.  Echo  phase  mea¬ 
surements  can  be  perform^  using  either  one  of  the 
following  techniques: 

1)  The  transmitter  pulse  is  produced  by  high  power 
coherent  amplification  of  a  pul%  sample  of  a  frequency 
stable  signal  (STALO).  The  STALO  signal  is  also  used 
as  a  phase  reference  for  incoming  echoes; 

2)  The  transmitter  is  a  pulsed  high  power  oscillator 
producing  an  RF  pulse  with  a  random  phase.  Either  a 
ftequency  stable  reference  is  used  to  measure  the  phase 
of  both  the  transmitter  pulse  and  the  echoes  associated 
with  it,  or  the  transmitter  pulse  synchronizes  a  coherent 
oscillator  (COHO)  acting  as  a  phase  reference  for  the 
incoming  echoes. 
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Fto.  t.  Atmospheric  (ases  ibsorplion  spectrum  at  the  ground  for  various  humidity  conditions 
indicated  by  the  specific  humidity  values.  The  absorption  due  to  a  100  mm  rainfall  is  also  shown. 


Phase  coherent  power  amplification  at  94  GHz  can 
be  implemented  with  a  high  power  ( I  k W  peak,  5  watts 
average),  low  phase  noise,  pulse  coherent  amplifier  re¬ 
ferred  to  as  an  Extended  Interaction  Amplifier  (ElA). 
The  EIA  appears  to  be  an  optimum  device  for  a  mil¬ 
limeter-wave  Doppler  radar  because  of  its  low  phase 
noise;  however,  its  cost  is  high.  A  much  less  expensive 
power  oscillator  tube,  the  Extended  Interaction  Oscil¬ 
lator  (ElO)  delivers  the  same  power  as  the  EIA  and, 
used  as  a  pulsed  oscillator  assisted  by  phase  locking 
circuits,  represents  a  more  economical  alternative.  The 
ElO  solution  was  selected  for  our  Doppler  radar  design 
based  on  a  phase-locking  COHO  technique. 

A  block  diagram  of  the  millimeter-wave  Doppler 
radar  is  shown  in  Fig.  2;  Fig.  3  shows  a  j^otograph  of 
the  equipment.  The  radar  operates  at  a  frequency  of 
93.95  GHz  and  is  equipped  with  two  (fw  transmission 
and  reception)  Cassegrain  antennae  (3-foot  diameter 
dishes,  0.27®  3  dB  beamwidth).  The  isolation  between 
the  two  antennae  exceeds  80  dB.  The  Cassegrain  feed 
configuration  allows  the  transmiuer  or  receiver  to  be 
attached  directly  to  the  back  of  its  respective  antenna 
so  that  waveguide  losses  (3  dB  rii~')  are  ntinimized. 
The  entire  radar  is  set  on  a  rigid  frame  that  con  be 
tilted  and  rotated  in  azimuth;  thereby  offering  antenna 
scanning  capabilities. 

b.  Transmitter 

Doppler  radars  with  magnetron  transmitters  have 
been  implemented  at  K,-band  (Pasqualucci  et  al.,  1983; 


Hobbs  et  al.,  1985).  A  magnetron  delivering  10  kW  at 
94  GHz  is  available  off-the-shelf,  but  its  choice  was 
rqected  because  of  short  life  expectancy  and  unknown 
intrapulse  phase  noise.  The  EIO  has  a  demonstrated 
longer  life  and  is  believed  to  have  lower  phase  noise, 
so  that  even  if  it  delivers  less  power  (1  kw  peak),  it  was 
considered  a  better  solution.  Also,  the  ElO  can  be  gated 
by  applying  a  low  power  pulse  to  an  input  grid  which 
controls  the  higL  power  output  pulse.  The  dc  high  volt¬ 
age  main  power  source  is  constantly  applied  to  the  tube, 
and  there  is  no  need  to  switch  on  a  large  amount  of 
power  to  produce  the  transmitter  pulse,  a  basic  re¬ 
quirement  for  pulsed  magnetron  operation.  Therefor^ 
frequency  stability  and  full  power  are  rapidly  estab¬ 
lished  at  ^e  start  of  the  transmitted  pulse  and  intrapulse 
frequency  modulation  is  leduced. 

Cbnsiderable  effort  was  devoted  to  the  modulator 
design.  This  and  the  ElO  performance  resulted  in  a 
transmitter  pulse  with  a  3  ns  rise  time  and  a  5  ns  fall 
time  and  intrapulse  frequency  modulation  and  phase 
noise  email  enough  to  ensure  satisfactory  Doppler  op¬ 
eration.  Figure  4  shows  an  example  of  the  transmitter 
pulse  at  the  output  of  the  mixer  used  in  the  phase-lock 
operation.  Note  the  very  short  rise  and  fall  times.  The 
pulse  width  can  be  set  between  0. 1  and  2  and  the 
pulse  repetition  rate  can  be  adjusted  between  2.5  and 
20  kHz.  Although  not  considered  in  this  early  design, 
a  transmitter  pulse  as  narrow  as  20  ns  (3-meter  range 
resolution)  can  be  produced.  Such  an  improvement  of 
range  resolution  would  require  an  increase  of  the  IF 
frequency  from  its  present  value  of  approximately  60 
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MHz  to  perhaps  160  MHz  or  more  to  accommodate 
the  need  for  a  larger  bandwidth  (50  MHz). 

c.  Receiver 

Phase-coherent  translation  of  rf  echoes  to  the  55 
MHz  IF  frequency  is  performed  by  a  low-noise  mixer- 
preamplifier  and  a  94-GHz  frequency  stable  local  os¬ 
cillator  (STALO).  The  STALO  is  a  Gunn  diode  phase 
locked  oscillator  (PLO),  referenced  to  a  100-MHz 
crystal  oscillator..  The  PLO  signal  has  the  frequency 
stability  and  low  phase  noise  of  the  crystal  oscillator 
reference.  The  mixer-preamplifier  noise  figure  is  6.5 
dB  double-sideband,  and  the  equivalent  receiver  noise 
power  is  P„  =  -93  dBm  (30  MHz  double  sideband 
as  a  preamplifier-filter  front-end  is  not  available  at 
94  GHz). 

A  separate  mixer-preamplifier  provides  a  sample  of 
the  transmitter  pulse  downconverted  to  IF,  needed  for 
COHO  synchronization.  The  94  GHz  rf  signal  input 
to  the  mixer  is  derived  from  a  20  dB  coupler  inserted 
between  the  ElO  output  and  the  transmitter  antenna; 
an  additional  80  dB  attenuation  in  an  attenuator- 
waveguide  assemHy  sets  the  RF  pulse  at  the  input  of 


Fio.  3.  Photograph  of  the  Doppler  radar^ 


the  mixer  to  a  -40  dBm  power  level.  The  IF  phase 
reference  pulse  (3  GHz  bandwidth)  at  the  mixer  output 
shown  in  Fig.  4  synchronizes  the  COHO  oscillator 
when  it  starts  to  oscillate.  The  COHO  continues  to 
oscillate,  retaining  the  phase  of  the  reference  pulse  du^ 
ing  the  interpulse  time,  until  it  is  disabled  a  few  mi¬ 
croseconds  before  the  next  transmitter  pulse. 

A  conventional  dual-mixer  system  acts  as  a  phase 
detector  producing  the  familiar  coherent  video  com¬ 
posed  of  two  in-phase,  I,  and  quadrature,  Q,  signals. 
Figure  5  presents  an  example  of  the  I  or  Q  signal  as¬ 
sociated  with  the  transmitter  pulse  as  a  test  of  radar 
phase  coherence.  Note  that  the  COHO  frequency  is 
intentionally  shifted  (approximately  2  MHz)  with  re¬ 
spect  to  the  IF  frequency  and  that  phase  coherence  is 
still  obtained.  The  low  phase  noise  performance  of  the 
radar  was  further  tested  by  observing  the  coherent  video 
for  fixed  targets  shown  in  Fig.  6. 

3.  Signal  processing,  radar  sensitivity  and  calibration 

a.  Signal  processing 

The  I  and  Q  coherent  video  signals  and  the  signal 
at  the  output  of  a  logarithmic  amplifier  are  sampled 
and  digitized  at  range  gate  positions  by  three  8-lMt  an- 


FiC.  4.  IF  signal  (5S  MHz)  at  the  output  of  the  mixer-preamplifier 
used  for  phase  lock.  Mixer  inputs  are  94  GHz  ElO  transmitter  pulse 
and  PLO  CW  signal.  Pulse  width  is  0.1  tiS. 
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Flo.  S.  1  or  Q  coherent  video  »|nal  at  the  output  of  the  phase 
demoduiator  occurring  during  the  transmitter  pulse.  PRF;  10  KHz, 
2-second  ezpoture  time,  pulse  width  1.5  iit,  COHO/IF  f,equency 
offset  2X)MKz. 


ak}g-t(MUgita]  converters.  The  converters  are  of  the 
flash  type  with  a  conversion  time  less  than  100  ns. 
Video  samples  can  be  recorded  at  100  ns  time  intervals 
if  a  fast  OufTer  memory  is  used,  and  Doppler  spectra 
calculated  off-line.  This  method  is  comtflex  and  leads 
to  a  relatively  slow  data  throughput  rate,  considerably 
smaller  flian  that  achieved  using  a  signal  processing 
autocovariance  technique  originally  proposed  by 
Woodman  and  Hagfon  (1969)  and  named  '*pulse-pair" 
(PP)  later.  In  the  first  phase  of  this  imject,  directed 
primarily  toward  exploring  the  millimeter  wave  Dopp¬ 
ler  radar  capal^ities  for  meteorolo^cal  research,  we 
have  selected  the  PP  method  to  provide  continuous 
pTOOcising  of  mean  Doppler  at  alt  available  range  gates. 
Our  PP  processor  requires  a  minimum  of  3(X)  ns  to 
perform  the  complex  multiplication  and  integration 
needed  for  calculation  of  the  autocovariance  in  real¬ 
time,  therefore  setting  the  limit  to  the  minimum  avail¬ 
able  time  interval  between  range  gates.  The  processor 
can  calculate  signal  autocovariance  at  all  the  range  gates 
available  in  the  pulse  repetition  period  up  to  a  maxi¬ 
mum  of  1024. 

The  PP  method  has  often  been  discussed  (i.e.,  Dov- 
iak  and  Zmic,  1983),  and  the  following  is  a  brief  sum¬ 
mary  of  the  theory  and  operation  of  the  PP  processor 
designed  for  this  project.  The  complex  autocont'ation 
function  of  the  bacl^ttered  signal  is: 

p(r)  =  pXr)e-^''^’  (1) 

where  r  is  the  time  lag  and  p,(r)  is  the  autocorrelation 
function  amplitude  only  dependent  on  signal  band¬ 
width.  The  fo  is  rigorously  equal  to  the  mean  frequency 
(first  specti^  moment)  for  symmetric  spectra  but  only 
very  riightly  dependent  of  spectrum  nonsyrametry 
(Sirmans,  1975). 

We  can  express  pj(t)  by: 

=  (2) 

where  is  the  spectral  variance.  Equation  (2)  applies 
rigorously  to  a  Gaussian  shape  spectrum;  however,  it 
is  a  very  close  approximation  for  most  spectrum  shapes 
(Lhermitte  and  Serafin,  1984). 

In  our  processor,  the  real  and  imaginary  terms  rep¬ 
resenting  the  complex  signal  autocovariance  for  a  time 
lag  equal  to  the  interpulse  time  and  the  +  Q^)  term 


representing  the  zero-lag  reference  are  calculated  in  real 
time  at  each  range  pte  and  recorded  on  magnetic  tape 
( 1 6-bit  words).  An  off-line  computer  performs  the  cal¬ 
culation  of^  and  following  Eqs.  (1 )  and  (2). 

Some  of  the  observations  present^  here  were  made 
in  very  low  signal-to-noise  conditions.  White  noise  does 
not  bias  the  estimate  of  ^  (it  oitly  increases  its  variance). 
However,  the  noise  contribution,  P„,  to  the  (1^  +  (}*) 
term  representing  the  zero  lag  autocovariance  amifli- 
tude  r.fierence  must  be  remo^  in  the  spectral  width 
co'nputations,  as  this  term  must  relate  to  signal  au¬ 
tocorrelation  function  amplitude  only.  The  long  signal 
integration  time  used  in  our  observations  allows  precise 
determination  of  Pn  in  echo  free  regions,  which  can 
thus  be  subtracted  from  ^1^  +  (^)  for  spectral  width 
computation.  We  have  used  this  proc^ure  for  the 
spectrum-width  data  presented  here.  However,  this 
method  only  removes  white  noise  contribution  from 
the  receiver.  All  white  noise  contributions  including 
the  presence  of  noise  in  backscatter  or  radar  phase  noise 
can  be  canceled  by  calculating  p|  and  pj  for  two  pulse- 
repetition  periods  Ti  and  tj,  TTie  signal  spectral  vari¬ 
ance  is  then  expressed  by  (Lhermitte  and  Serafin,  1984) 

<ri'“ln(p,/P2)/(2ir’(r,'-T:')]  (3) 

This  improved  noise  correction  method  will  be  im¬ 
plemented  later  when  the  radar  is  equipped  with  dual 
pulse  repetition  rate  capabilities. 

b.  Atmospheric  and  hydrometeor  absorption 

The  absorption  spectrum  for  atmospheric  gases  is 
shown  in  Fig.  1 .  The  3-mro  wavelength  is  inside  a  win¬ 
dow  that  contains  the  skirts  of  the  60  GHz  and  1 18 
GHz  O3  lines,  and  also  the  skirt  of  the  intense  184 
GHz  H2O  line  which  makes  absorption  in  this  window 
strongly  dependent  on  atmospheric  water  vapor.  Figure 


a 


b 


V 


FtG.  6.  Coherent  video  ji^al  (I  or  Q)  on  6xed  targeu.  a.  targets  at 
7S0  meters  and  1  km;  b.  targets  at  ~2  km  on  an  expanded  scale. 
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7  shows  the  absorption  coefficient  at  94  GHz  derived 
from  Liebe's  ( 1 985)  data,  as  a  function  of  specific  hu-  ^ 
midity.  One  sees  that  with  25  g  ra"^  humidity,  the  ab-  I 
sorption  coefficient  is  about  2  dB  km" '  (one-way).  For  W 
a  vertically  pointing  radar  beam,  zenith  absorption  is  ® 
much  less  agnificant  because  of  the  systematic  decrease  S 

of  both  water  vapor  amount  and  pressure  broadening,  i 
with  altitude.  Shimabukuro  and  Epstein  (1970)  have  | 
collected  a  large  number  of  observations  of  clear  air  ^ 
attenuation  at  94  GHz  using  solar  radiation  as  a  source;  !< 
an  empirical  least  squares  fit  of  their  measurements  of  ; 
zenith  absorption  versus  precipitable  water  is  shown  | 
in  Fig.  8.  The  same  data  computed  by  Liebe  (1985) 
from  laboratory  experiments  are  also  presented  for 
comparison. 

For  small  spheres  of  diameter  A  the  absorption  eross 
section  Qa  =  (ir^D^/\)  Im(-/r)  (Rayleigh  approxima¬ 
tion),  where  K={m^  -  I  )/(m’  +  2)  and  m  is  the  com¬ 
plex  index  of  refraction.  At  94  GHz,  m  =  3.5  +  2.05/ 
for  water  (25®C)  and  m=  1.78  +  4.76  X  10"^/ for  ice 
(Ray,  1972),  so  lra(~K)  »  0.1402  for  water  and  2.5 
10"^  for  ice.  The  d,  is  proportional  to  particle  mass  so 
the  absorption  coefficient  is  proportional  to  cloud  or 
fog  water  content.  The  calculated  one-way  absorption 
per  g  m"^  of  liquid  water  varies  almost  linearly  with 
temperature  from  3.47  dB  km"‘  at  23®C  to  4.8  dB 
km"'  at0*C 

Precipitation  particles  whose  diameter  exceeds  ap¬ 
proximately  0.2  mm  are  not  Rayleigh  absorbers  at  94 
GHz,  and  d,  must  be  calculated  using  the  Mie  equa¬ 
tions.  Furthermore,  the  scattering  cross  section  Q,  ex¬ 
ceeds  di  for  water  particles  with  a  diameter  greater 
than  1  mm.  At  94  GHz,  the  attenuation  coefficient 
(including  scattering)  ^  (in  dB  km"'),  calculated  for 
rainfall  and  assuming  a  Marshall-Palmer  (Marshall 


Fig.  7.  Atmospheric  atienualion  ai  94  GHz  as  a  function  of  specific 
humidity  in  g  ro~’,  at  ground  pressure  and  for  two  lemperaiures.  280 
and  310  K. 


Flo.  8.  One-way  zenith  absorption  as  a  function  of  precipitable 
water.  The  solid  line  is  the  model  proposed  by  Shimabukuro  as  a 
least  squares  fit  of  his  data.  The  heavier  line  indicates  the  range  of 
Ihe  Shimabukuro  data.  The  circles  are  the  values  proposed  by  Liebe 
on  the  basis  of  his  experiments.  NtHe  the  very  good  agreement  between 
two  completely  independent  sets  of  data. 


and  Palmer,  1948)  drop  size  distribution  N(,D) 
=  A^o  exp(-AA  A’o  =  8  •  10'‘  and  A  =  R  rainfall 
in  mm  h"',  is  approximately  /3  =  0.9R^^*-,  this  equation 
agrees  reasonably  well  with  experimental  measure¬ 
ments. 

c.  Rayleigh-Mie  backscattering 

Water  or  ice  spheres  with  a  diameter  smaller  than 
approximately  k/4  are  Rayleigh  scattereis  and  theii  ra¬ 
dar  eross  section,  S„  is 

5,=ir*|/(iW.  (4) 

At  94  GHz,  lATp  =  0,813  for  water  (temperature 
20®C)  and  0. 1 73  for  ice.  Since  l/iTp  docs  not  vary  sig¬ 
nificantly  for  longer  wavelengths,  S,  for  Rayleigh  scat- 
terers  at  W-band  is  1 7  dB.  This  is  40  dB  higher  than 
that  at  Ka  and  X-band,  respectively.  The  Sr  values  cal¬ 
culated  at  94  GHz  using  the  Mie  equations  are  shown 
in  Fig.  9  as  a  function  of  particle  diameter,  D,  for  water 
(0®  and  20®C)  and  ice.  Rayleigh  scattering,  which  is 
acceptable  for  Z)  >  8  mm  particles,  is  also  indicated  in 
Fig.  9..  Above  Z)  ~  1.8  mm,  the  Mie  backscattering 
values  oscillate  and  are  consistently  1 5  to  25  dB  below 
Rayleigh  backscattering.  The  Mie  computations  for  K«- 
band  and  X-band  are  also  shown  in  Fig.  9  for  com¬ 
parison. 

For  Rayleigh  scatterers,  S,  is  6.7  dB  less  for  ice  than 
it  is  for  water.  However,  at  94  GHz,  precipitation  par¬ 
ticles  with  D  >  I  mm  are  more  efficient  scatterers  if 
they  are  made  of  ice  rather  than  water  (see  Fip  This 
fact  was  recognized  tong  ago  for  centimetr  uars  de¬ 
tecting  large  hailstones  (Atlas  et  al.,  1 96U).  However, 
at  3-mm  wavelength  the  increased  backscatter  for  ice 
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Fig.  9.  Backscattering  cross  section  of  water  and  ice  spheres  versus 
their  diameter  at  94  GHa.  The  results  for  35  GHz  and  9.3  GHz  are 
also  indicated  for  reference. 


occurs  for  millimeter  size  particles  commonly  found 
in  clouds.  As  a  comparison,  the  increased  backscatter 
for  ice  at  K«-band  occurs  for  particle  sizes  £>  >  7  to  8 
mm,  definitely  greater  than  the  size  of  ice  particles  in 
clouds  except  for  large  graupel  and  hailstones. 

Radar  reflectivity  is  expressed  by  tj  =  JN(D) 
X  S^D)dD,  where  N{D)  is  the  drop  size  distribution. 
The  reflectivity  factor  is  Z  =  / N(D)D^dD  and  implies 
Rayleigh  backscattering.  At  94  GHz  and  assuming 
Rayleigh  scattering  the  relationship  between  Z in  mm^ 
m“^  and  in  cm“'  is 

Z=4.21X10V  (5) 

Precipitation  particles  at  94  GHz  are  Mie  scatterers 
and  an  equivalent  radar  reflectivity  factor,  Z^,  derived 


from  Eq.  (5)  using  the  actual  t)  calculated  using  Mie 
backscattering  cross  sections,  must  be  expressed. 

The  Zr  computations  at  94  GHz,  bas^  on  a  Mar¬ 
shall-Palmer  dropsize  distribution,  are  presented  as  a 
function  of  rainf^l  rate,  R,  in  Fig.  10  in  mm  h~'.  For 
Rayleigh  scatterers  Z  is  also  shown  and  one  sees  that 
Z,  is  always  smaller  than  Z,  with  Z/Z,  increasing  to 
25  dB  for  .R  =  100  mm  h”'.  Since  the  ratio  S^M,  where 
M  is  the  particle  mass,  decreases  with  increasing  D,  a 
narrow  dropsize  distribution  concentrated  around  the 
virtual  maximum  raindrop  size  (5-6  mm),  produces 
even  slightly  less  Z,  values  for  the  same  precipitation 
intensity.  The  3-mm  wavelength  radar  reflectivity  of 
rainfall  cannot  be  greater  than  approximately  1.5 
X  10"^  cm"’,  or  an  equivalent  Z,  of  30  dBZ,  ■'s  com¬ 
pared  with  55  dBZ  values  or  more  for  centimeter 
wavelength  radars,  .(..arge  ice  particles  (i.e.  hailstones) 
can  increase  Z,  to  note  than  30  dBZ,  still  less  than 
the  radar  reflectivity  observed  using  centimeter-wave 
radars.  The  fact  (hat,  at  94  GHz,  large  raindrops  do 
not  dominate  backscatter  renders  the  relationship  be¬ 
tween  radar  reflectivity  and  rainfall  rate  much  more 
sensitive  to  small  particle  deficiency  (compared  with  a 
M-P  distribution)  often  observed  in  drop  size  distri¬ 
bution  measurements.  This  may  invalidate  the  use  of 
a  M-P  drop  size  distribution  for  the  interpretation  of 
rack;’'  reflectivity  data  at  a  3-mm  wavelength. 

The  expression  for  mean  Doppler  velocity,  V,  at 
vertical  incidence  is 


N{D)SXD)V{D)dD 


// 


NiD)S^D)dD 


(6) 


where  V(D)  is  the  particle  vertical  velocity.  Assuming 
that  precipitation  particles  fall  at  their  terminal  velocity 
and  using  a  M-P  dropsize  distribution  and  the  Mie 
scattering  function  allows  calculation  of  the  V-R  re- 


Fig.  10.  Radar  reflectivity  factor,  Z,  as  a  function  of  rainfall  rate 
for  a  long-wavelength  radar  (Rayleigh  scattering)  and  a  94-GHz  radar. 
Marshall- Palmer  dropsize  and  the  Gunn-Kinzer  vertical  velocity/ 
diameter  relationship  are  used. 
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lationship  at  94  GHz  shown  in  Rg.  1 1  for  rainfall  at 
sea  level,  together  with  the  V-R  relationship  proposed 
byJossand  Waldvogel  (1977)  for  microwave  radar.  At 
94  GHz  P'for  very  heavy  rainfall  (100  mm  h~')is5m 
s"'  instead  of  9-10  m  s"'  for  the  J-W  relationship,  ow 
ing  to  the  relatively  stronger  contribution  to  bacl^t* 
tering  (compared  to  microwave)  by  smaller  particles. 

d.  Radar  equation  and  radar  parameters 

With  the  acceptable  assumption  that  meteorological 
targets  fill  the  radar  beam  completely,  radar  reflectivity 
can  be  expressed  by  the  following  simplified  rads^ 
equation,  presented  in  a  logarithm  form  for  conve> 
nience. 

1 0  logif  =  1 0  logPr  - 1 0  logP,  + 1 0  log4ir/?  ^ 

-lOlog/i-lOlog/i,  (7) 

where  v  is  the  radar  reflectivity  of  the  meteorological 
scatterers,  P,  the  echo  power,  P,  the  transmitter  peak 
power,  A  =  tcfl  the  scattering  volume  radial  dimension. 
At  the  effective  area  of  the  receiving  antenna,  and  R 
the  target  distance. 

The  transmitter  antenna  efficiency  including  the  loss 
of  signal  in  the  antenna  sidelobes  must  be  included  in 
At  as  well  as  the  loss  of  sensitivity  arising  from  the  only 
partial  merging  at  close  range  of  essentially  parallel  an¬ 
tenna  beams.  The  antenna  radiation  pattern  is  not  ac¬ 
curately  known,  but  it  is  unlikely  that  the  antenna  ef¬ 
ficiency  and  beam  merging  correction  exceeds  3.S  dB. 

With  this  correction  included,  and  using  P,  =  1  kW 
peak,  a  3-foot  antenna  diameter,  and  A  -  60  m  (400 
ns  pulse),  the  radar  reflectivity  is 

101og»r  =  P,(dBm)-23.5  +  201ogP  (8) 

where  r;  is  in  cm”',  P,  is  the  received  signal  intensity 
in  dBm,  and  R  the  target  distance  in  km. 

Or,  in  dBZ: 

dBZ  =  P,(dBm)  +  53  +  20logP.  (9) 


FlC.  1 1 .  Mean  Doppler  velociry  at  vertical  incidence  as  a  function 
of  rainfall  lolensity  for  cenlimeter  wave  radar  (Rayleigh  scattering 
accepted),  and  for  940Hz  radar  data.  M-P  dro|»ize  distribution 
and  Ounn-Kinzer  vertical  velocity/size  relationship  are  used.  The 
spectrum  width,  v,,  for  the  94  GHz  data  is  also  presented. 


At  a  3  km  range,  a  cloud  reflectivity  of  approximately 
-30  dBZ  will  produce  a  signal  equal  to  noise  (-93 
dBm);  this  is  the  radar  reflectivity  at  94  GHz  of  a  1  g 
m”^  liquid  water  cloud  composed  of  10  particles. 

e.  SignaUo-noise  improvement  and  radar  receiver  cal¬ 
ibration 

This  section  discusses  the  effective  signal-to-noise 
gain  due  to  signal  intention,  and  signal  intensity  cal¬ 
ibration  methods.  The  results  presented  here  were  ob¬ 
tained  u.sing  a  3.2-sec  signal  dwell  time  and  a  10  kHz 
PRF,  providing  32  000  samples  for  the  compulation 
of  signal  autocovariance  and  signal  intensity  for  each 
range  gate.  This  unusually  large  sample  size  results  in 
very  small  variance  of  the  estimates  and  enables  mea¬ 
surement  of  signal  intensities  which  are  much  less  than 
the  receiver  noise. 

Mean  intensity  which  yields  radar  reflectivity  is  the 
mean  square  of  the  backscatter  amplitude.  However, 
mean  signal  intenaty  is  usually  derived  from  integrat¬ 
ing  the  signal  at  the  output  of  a  radar  receiver  and,  if 
the  radar  receiver  input/output  transfer  function  is  not 
exactly  quadratic,  the  receiver  calibration  performed 
using  a  nonfluctuating  signal  must  be  corrected.  For  a 
logarithmic  amplifier,  this  correction  is  -f'2.5  dB  for 
large  signals  and  O.S  to  1  dB  for  weak  signals  for  which 
the  logarithmic  amplifier  transfer  function  becomes  al¬ 
most  linear  (Lhermitte  and  Kessler,  1 966).  In  our  radar 
processor,  an  arbitrary  correction  of  +0.5  dB,  below 
-80  dBm  and  +2.5  dB  above  -80  dBm,  is  applied  to 
the  logarithmic  signal  integrator  output.  The  I  and  Q 
signals  below  saturation  are  produced  by  an  essentially 
linear  receiver-phase  detector  transfer  function  and  the 
<1^  +  Q^)  integration  data  need  not  be  corrected. 

A  calibrated  signal  source  at  94  GHz  was  not  avail¬ 
able  for  calibration  of  the  radar  receiver  and  we  had 
to  consider  an  alternate  radar  receiver/integrator  cal¬ 
ibration  technique.  Let  P„  be  the  94-GHz  noise  power 
and  Ps  the  received  signal  at  the  input  of  the  millimeter 
wave  receiver.  Assuming  independent  (Gaussian)  pro¬ 
cesses  the  signal  plus  noise  power  is  Pj+it  “  Ft  +  P„. 
The  change  from  P„  to  associated  with  the  oc¬ 
currence  of  a  signal,  can  be  expressed  by  the  ratio 
“  =  Ps^nlPn-  The  signal-to-noise  ratio  is  thus  Ps/P„ 
=  a  -  1 ,  or: 

(dBm)„^  =  (dBm)noi«.+  10log(a-l).  (10) 

The  relationship  between  a  and  PJP„  is  shown 
graphically  in  Fig.  12.  If  the  integrator  output,  F,,  as¬ 
sociated  with  an  echo  exceeds  the  noise  rms  fluctuation, 
the  signal  is  detected.  With  our  dataset  (32  000  samples 
at  each  range  gate,  see  below)  the  detection  threshold 
was  found  to  be  -20  to  -25  dB  below  the  noise.  This 
result  agrees  well  with  the  \/(Nt)'^^  theoretical  signal- 
to-noise  improvement  expected  from  Nt  statistically 
independent  samples. 

Only  the  signal-to-noise  ratio,  PJP„,  is  measured  by 
the  integrator  and  the  noise  power  P„  must  be  assessed 
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can  be  established.  The  accuracy  in  measurement  of 
receiver  noise  power  and  the  knowledge  of  the  receiver 
transfer  function  are  the  main  factors  governing  the 
accuracy  in  measuring  mean  signal  intensity.  For  P, 
greater  than  we  believe  that  the  uncertainty  in  our 
P,  measurement  does  not  exceed  ±2  dB  and  can  de¬ 
grade  to  ±3  dB  for  S/N  *  -20  dB. 

f.  Precision  of  mean  velocity  and  spectral  width  esti¬ 
mates 

The  ability  of  the  radar  to  provide  reliable  measure¬ 
ment  of  mean  Doppler  and  spectrum  width  for  a  signal 


Fio.  12.  Signal'to-noise  ratio  S/N  aa  a  function  of  the 
(S  +  N)/N  ratio  obsen'ed  by  the  radar. 

for  the  received  signal  intensity  to  be  measured. 
Knowing  the  receiver  noise  figure  and  bandwidth,  P„ 
can  be  calculated.  A  laboratory  measurement  of  the 
mixer-preamplifier  noise  figure  can  be  performed,  but 
it  requires  extensive  equipment  and  a  more  convenient 
way  is  to  measure  the  receiver  performance  using  solar 
radiation  as  a  reference  source.  The  quiet  sun  radiation 
temperature  is  8000®K  at  94-GHz.  The  sun  can  be 
considered  as  a  blackbody  emitting  a  signal  power  (ra¬ 
diation  flux)  14.2  dB  above  a  3(X)'’K  reference.  The 
radar  antenna  3  dB  beamwidth  is  approximately  0.27”, 
and  the  solar  disc  diameter  is  (0.53”).  If  the  radar  beam 
axis  is  centered  on  the  solar  disc  and  assuming  a  gaus- 
sian  shape  radar  beam  and  a  uniform  sun  brightness, 
yields  less  than  2%  (0. 1  dB)  signal  loss  with  respect  to 
a  radiation  source  filling  the  beam  completely. 

Using  an  integrator  with  approximately  O.S-sec  in¬ 
tegration  time  (a  few  degrees  Kelvin  precision  with  30 
MHz  receiver  Irandwidth),  it  was  found  that  when  the 
antenna  was  pointed  toward  the  solar  disc,  the  solar 
radiation  power  level  finally  reaching  the  mixer  was 
(fortuitously)  equal  to  the  mixer  noise  within  ±0.2  dB. 
Neglecting  the  lower  temperature  atmospheric  radia¬ 
tion,  we  calculated  the  solar  signal  power  at  the  input 
of  the  radar  receiver  by  considering  only  solar  signal 
attenuation  in  the  sun-radar  receiver  path.  We  esti¬ 
mated;  1 .5  dB  atmospheric  attenuation  from  temper¬ 
ature  and  humidity  soundings,  2  dB  ±  1  dB  loss  for 
antenna  efficiency  including  sidelobe  effects,  1  dB  loss 
in  waveguides,  and  3  dB  for  natural  radiation  received 
by  a  polarized  receiver,  or  a  total  of  7,5  dB.  This  yields 
a  6.6  ±  1.5  dB  noise  figure  which  agrees  with  the  6.5 
dB  quoted  by  the  manufacturer  of  the  mixer.  In  ad¬ 
dition  to  its  simplicity,  this  method  allows  easy  periodic 
checks  of  the  mixer  performance. 

A  6.5  dB  noise  figure  and  a  30  MHz  double  sideband 
bandwidth  is  equivalent  to  a  -93  dBm  receiver  noise 
power,  P„.  Knowing  P„,  F,  can  be  derived  from  the 
PJP„  measurements  and  an  absolute  calibration  of  the 
radar  receiver  from  antenna  input  to  integrator  output 


HEIGHT  (Km) 

Fig.  !3.  Power  terra,  /?„  =  +  Q*),  and  aulocovariance  real, 

A 12,  and  imaginary  Bi2  terms  for  weak  signals  with  signal-lo-noise 
indicated  for  R, , .  32  OCO  samples  contribute  to  the  estimates  evaluated 
at  each  range  gale.  The  mean  Doppler  estimates  are  also  shown. 
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well  below  the  receiver  noise,  is  illustrated  in  Fig.  13 
which  shows  weak  echoes  at  an  altitude  of  ^proxi- 
mately  7  km.  The  real  and  imaginary  autocovaiiance 
terms  are  shown  together  with  the  <1^  +  Q^)  power 
term.  Fig.  13  shows  the  well>defined  riean  Dopier 
velocity  (0.5  to  0.7  m  s"'  downward)  obtained  for  two 
echoes  (-96  dBwj  and  —113  dBm).  T  e  fluctuation 
level  of  the  autocovariance  real  and  imaginary  terms 
indicates  that  the  mean  Doppler  uncertainty  is  about 
±1.5  and  ±16  cm  s~'  for  the  -96  dBm  and  -113  dBm 
echo,  rcH)ectively.  The  theoretical  expression  for  the 
variance,  of  a  PP  mean  Doppler  estimate  for  large 

S/N  is  [Zmic,  1979]: 

2  _  /1 1  \ 

where  is  the  Doppler  spectrum  standard  deviation, 
X  the  radar  wavelength,  and  T  the  signal  dwell-time. 
With  X  =  3  mm,  T  -  3.2  s,  and  ==  1  m  s”‘,  a„ 
-  0.9  cm  s”'. 

This  theoretical  precision  of  the  mean  Dpppler  es¬ 
timate  was  compared  with  the  time  stability  of  actual 
measurements  of  mean  Doppler  obtained  in  stratified 
cloud  conditions  for  which  steadiness  of  the  velocity 
profiles  is  expected.  Two  examples  of  data  obtained  in 
such  stratifi^  conditions  are  shown  in  Fig.  14.  Persis¬ 
tent  fine  structures  in  the  two  profiles  are  observed. 
With  an  echo  translation  speed  of  5  m  s”‘  and  a  radar 
beam  cross  section  of  5  to  50  meters  ( 1  to  1 0  km  range), 
the  profiles  should  be  nearly  statistically  independent. 
Above  the  0®C  level  (4  km  altitude)  where  particles 
are  expected  to  be  ice  crystals,  60%  of  the  profile-to- 
profile  peak  velocity  differences  are  less  than  3  cm  s"*, 
and  85%  are  less  than  6  cm  s“*.  Between  the  4  km 
altitude  and  ground  level,  in  a  region  where  liquid 
droplets  are  expected,  60%  of  the  same  velocity  differ¬ 
ences  are  less  than  2  cm  s"‘  and  90%  less  than  5  cm 
s"\-  There  are  signs  of  convection  above  6  km  and 
some  of  the  velocity  variability  may  be  due  to  actual 
change  of  the  velocity  in  5  seconds,  Further  statistical 
analysis  of  profile  to  profile  velocity  differences  indi¬ 
cates  that  the  standaid  deviation  of  the  5-second  ve¬ 
locity  difference  is  typically  less  than  1.3  cm  s“'  in  S/ 
N  >  10  dB  regions.  This  is  slightly  greater  than  that 
predicted  by  Eq.  (1 1)  and  may  be  explained  by  actual 
time  variability  of  the  mean  vertical  velocity. 

4.  Observations 

In  this  section  we  present  selected  observations  to 
illustrate  the  capabilities  of  the  3  mm  wavelength 
Doppler  radar.  The  radar  was  originally  assembled  and 
tested  in  Miami.  Observations  of  thunderstorms  in  the 
Miami  area  were  made  just  afler  testing  of  the  radar. 
The  thunderstorms  were  at  a  20  to  30  km  range.  With 
the  high  humidity  conditions  (6  cm  of  precipitable  wa¬ 
ter  and  28  g  m“^  specific  humidity  at  ground  level). 


Fig.  14.  Vertical  profiles  of  mean  Doppler  observed  in  a  slraiifted 
cloud  which  may  contain  very  light  precipitation.  The  two  profiles 
are  obtained  from  processing  of  two  independent  sets  of  data  taken 
S  seconds  apart.  Note  the  sudden  acceleration  of  the  mean  Doppler 
velocity  at  the  3.9  km  altitude  associated  with  particle  melting. 


the  one-way  zenith  absorption  was  almost  3  dB  and 
the  one-way  signal  attenuation  at  the  ground  was  2.5 
dB  km“’..:  The  two-way  absorption  between  a  thun¬ 
derstorm  at  20  km  and  the  radar  was  1 2  dB  for  a  30" 
elevation  angle,  increasing  to  50  dB  for  a  T  elevation 
angle.  In  such  high  humidity  conditions,  the  minimum 
elevation  angle  for  which  echoes  were  deteaed  was  20® 
to  30°  so  only  the  top  of  thunderstorms  could  be  ob¬ 
served.  For  example,  a  storm  at  20  km  observed  with 
a  34°  elevation  angle  (altitude  in  the  storm  1 1.2  km), 
produced  a  -71  dBm  echo.  With  an  estimated  two- 
way  clear  air  absorption  of  1 1  dB,  and  an  additional 
estimated  attenuation  of  10  dB  (2  dB  km~'  two-way, 
5  km  depth)  in  an  intervening  cell  revealed  by  the  pres¬ 
ence  of  an  echo,  the  backscattered  signal  intensity  in 
the  storm  core  was  estimated  to  be  -50  dBm.  This  is 
equivalent  to  a  26  dBZ  radar  reflectivity  which  is  rea¬ 
sonable  for  targets  at  a  height  of  1 1  km  in  a  Florida 
thunderstorm,  and  also  considering  the  virtual  20  to 
28  dBZ  maximum  which  can  be  observed  at  94  GHz. 
Scanning  horizontally  and  vertically  in  search  of  the 
strongest  echo  in  the  storms,  indicated  that  the  maxi¬ 
mum  echo  intensity  at  a  mean  range  of  20  km  was 
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about  -70  to  -75  dBm  and  was  observed  at  30  to  40* 
elevation  angles. 

Observations  using  a  vertically  pointing  beam  were 
made  later  at  Miami  and  also  at  Boulder,  Colorado  (in 
association  with  the  National  Center  for  Atmospheric 
Research).  Cirrus  layers  were  always  detected  with  the 
millimeter-wave  radar,  even  if  they  were  not  observed 
visually.  This  must  be  due  to  the  fact  that  most  ice 
particles  are  large  and  likely  to  be  efficient  Mie  scat- 
terers  at  the  short  3-mm  wavelength  (see  Fig.  9).  Ice 
particles  in  cirrus  are  prisms  with  a  typical  len^  of 
0.5  mm  and  a  ratio  of  length  to  breadth  varying  from 
one  to  five  (sec  the  classical  study  by  Weickman,  1947, 
for  instance)  and  clusters  of  crystals  reaching  a  size  of 
2  mm  can  be  found.  Note  that  a  cirrus  cloud  containing 
one  spherical  ice  particle  of  2  mm  size  (Mie  backscat- 
tering  cross  section  2  mm^)  per  m^  has  a  2  X  10'*  cm”  ‘ 
reflectivity  equivalent  to  approximately  0  dBZ  at  X 
*  3  mm.  In  comparison,  a  fair  weather  cumulus  having 
the  same  water  content  and  with  droplet  size  distri¬ 
bution  extending  to  10  to  15  m  has  a  reflectivity  equiv¬ 
alent  to  -40  to  -35  dBZ. 

Observations  of  vertical  velocity  in  cirrus  clouds  at 
altitudes  of  7-10  km  indicated  that  particles  are  typi¬ 
cally  moving  downward  with  a  velocity  on  the  order 
of  0.5  to  1  m  s'*,  sometimes  reaching  2  m  s"'..  Occa¬ 
sionally,  upward  motion  of  approximately  1  m  s~'  was 
observ^  but  only  at  the  cirrus  top  with  downward 
motion  observed  ^low.  Since  cirrus  clouds  are  so  well 
observed  at  X  =  3  mm,  a  Doppler  radar  operating  at 
this  wavelength  may  be  an  excellent  tool  for  observa¬ 
tions  of  cirrus  kinematics  and  could  provide  insight 
into  high  altitude  turbulence. 

Radar  observations  using  the  same  vertically  point¬ 
ing  radar  beam  mode  were  also  made  at  the  Miami 
site  on  a  cloud  apparently  related  to  the  outflow  of  a 
dissipating  thunderstorm.  The  thunderstorm  was  lo¬ 
cated  approximately  30  km  West  of  the  radar  site  and 
the  first  echo  seen  by  the  radar  came  from  a  high  al¬ 
titude  (12  to  14  km)  cirrus  cloud  which  was  probably 
an  extension  of  the  storm  anvil.  Later,  the  echo  de- 
veltved  downward  and  finaUy  reached  the  ground.  The 
downward  development  of  the  cloud  seemed  to  be  re¬ 
lated  to  the  sinking  motion  of  ice  crystals  which  ap¬ 
peared  first  in  the  cirrus  layer  10  to  12  km  over  the 
radar..  There  was  never  any  significant  precipitation  on 
the  ground,  although  a  few  sparse  small  drops  (diameter 
less  than  1  mm)  were  observed  a  few  minutes  before 
the  observations  were  terminated. 

Stratified  precipitation  is  usually  characterized  by 
the  presence  of  a  well  defined  “bright  band”  (Austin 
and  Bemis,  1950),  or  more  properly  “melting  band” 
(MB),  whidi  appears  as  a  maximum  of  radar  reflectivity 
50  to  1 50  m  below  the  0®C  level.  This  enhancement 
of  radar  reflectivity  is  attributed  partly  to  the  increase 
of  the  iXTp  term  in  the  expression  of  Sr  associated  with 
the  change  from  ice  to  water.  The  passage  of  precipi¬ 
tation  particles  through  the  0*C  altitude  level  is  also 


characterized  by  a  sudden  increase  of  their  terminal 
velocity  as  they  evolve  from  a  high  drag,  low  density 
ice  cry^-snowflake  structure  to  a  high  density  droplet 
or  raindrop  (i.c.,  Lhcrmittc  and  Atlas,  1963). 

Although  many  MB  oteervations  with  centimeter- 
wave  radars  have  been  done  in  the  past,  we  were  de¬ 
listed  to  find  that  a  MB  could  be  clearly  identified 
using  a  3-mm  wavelength  radar  for  which  particles 
larger  than  0.8  mm  are  not  Rayleigh  scatterers.  Samples 
of  mean  Doppler  velocity  vertical  profiles  observed  with 
our  millimeter-wave  r^r  in  the  cloud  mentioned 
above  are  presented  in  Fig.  1 4  and  also  in  Fig.  1 5  where 
vertical  profiles  of  mean  signal  intensity  range,  nor¬ 
malized  to  1  km,  and  spectral  variance  are  also  shown. 
The  most  striking  feature  in  the  velocity  profiles  is  in¬ 
deed  the  drastic  change  of  mean  Doppler  from  ap¬ 
proximately  1 .2  m  s“*  to  approximately  4.5  m  s~',  when 
particles  go  through  the  melting  level.  These  observa¬ 
tions  were  made  at  a  time  when  the  parent  storm  had 
almost  completely  dissipated  and  the  vertical  velocity 
profiles  were  fairly  steady.  This  is  confirmed  by  in¬ 
spection  of  Fig.  14  which  exhibits  remarkable  similarity 
for  two  velocity  profiles  taken  five  seconds  apart. 

At  the  altitude  where  particle  melting  is  suspected 
from  the  increase  of  mean  Doppler,  there  is  a  sudden 
6  to  8  dB  reflectivity  increase.  Further  inspection  of 
many  profiles  acquired  successively  at  five-second  time 
interval,  indicates  that  this  very  localized  reflectivity 
increase  attributed  to  melting  was  between  6  and  9.5 
dB  with  a  remarkably  repeatable  sl(^  of  0.6  dB/10 
m.  The  presence  of  such  a  reflectivity  enhancement 
attributed  to  melting  indicates  that  the  particles  are 
Rayleigh  scatterers,  i.e.  they  are  smaller  than  1.5  mm 
(sec  Fig.  9).  Detail  study  of  the  reflectivity  and  velocity 
vertical  profiles  indicates  the  sharp  reflectivity  change 
occurs  at  30  to  50  m  above  the  altitude  at  which  par¬ 
ticles  start  to  accelerate,  and  is  completed  at  an  altitude 
where  a  2.8-2.9  m  s“'  mean  Doppler  (middle  of  the 
accelerating  curve)  is  observed.  There  is  no  significant 
decrease  of  reflectivity  just  below  the  MB  as  seen  in 
classical  observations  of  MB  by  centimetric  wavelength 
radars  (i.e.,  Lhermitte  and  Atlas,  1963). 

Previous  observations  of  the  MB  in  precipitation 
conditions  showed  a  systematic  downward  increase  of 
radar  reflectivity  starting  a  kilometer  or  so  above  the 
0®C  level  MB.  This  is  not  observed  here.  Rather,  a 
continuous  downward  reflectivity  decrease  of  both  sig¬ 
nal  intensity  and  vertical  velocity  is  observed  starting 
a  few  hundr^  meters  above  the  MB.  This  may  indicate 
that  the  cloud  particles  we  observe  just  above  the  MB 
are  single  crystals  (or  more  likely  prisms)  having  a  small 
collision  efficiency  with  the  surrounding  supercooled 
droplets  because  of  their  smooth  and  simple  shape  and 
their  small  terminal  velocity.  On  the  contrary,  the  large 
snowflakes  observed  using  a  microwave  radar  in  pre¬ 
cipitation  conditions  have  a  complex  structure  and  a 
substantial  terminal  velocity.  The  snowflakes  reside  in 
a  high  liquid  water  environment  and  they  are  likely  to 
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SPECTRUM  MEAN  DOPPLER  VELOCITY  m»-'  NORMAL.  SIGNAL 
WIDTH  INTENSITY 


Fig.  is.  Vertical  profiles  of  radar  signal  intensity  corrected  for  range  (right),  mean  Doppler 
velocity  (center),  and  spectrum  width  (left).  Profiles  labeled  2  are  obser/ed  10  minutes  after 
profiles  labeled  I. 


grow  rapidly  through  capture  (riming)  of  the  sur¬ 
rounding  cloud  droplets,  which  results  in  an  increase 
of  both  radar  reflectivity  and  mean  Doppler  velocity 
in  the  layer  just  above  the  MB. 

Our  millimeter- wave  radar  observations  indicate  that 
at  higher  altitudes  above  the  MB,  the  velocity  and  re¬ 
flectivity  profiles  remained  fairly  steady  over  a  long 
period  of  time  (one  hour).  However,  below  the  MB, 
there  are  significant  time  variations  of  the  velocity  and 
reflectivity  profiles  over  5  to  10  minute  periods.  This 
is  illustrated  by  the  profiles  presented  in  Fig.  15  which 
are  observed  10  minutes  apart.  In  the  profiles  labeled 
1,  the  mean  Doppler  velocity,  the  radar  reflectivity, 
and  the  spectrum  width  below  the  MB  do  not  exhibit 
much  vanation  with  altitude.  The  spectral  width  has 
an  average  value  of  (1.2  m  s'*')  which  is  comparable 
with  the  value  predicted  for  a  few  millimeter  per  hour 
M-P  rainfall  not  observed  at  the  ground.  This  relatively 
high  spectrum  width  may  be  due  to  the  presence  of  a 


large  concentration  of  cloud  droplets  coexisting  with 
another  much  smaller  population  of  0.5  mm  to  I  mm 
diameter  raindrops  falling  at  4  to  6  m  s"'. 

The  vertical  profiles  labeled  2  which  are  observed 
10  min  after  profiles  I  are  quite  similar  to  profiles  1 
above  the  MB.  However,  below  the  MB,  velocity  profile 
2  exhibits  greater  velocity  (5-6.5  m  s"')  than  that  ob¬ 
served  in  profile  1 ,  and  also  much  greater  velocity  and 
reflectivity  variations  with  height.  Just  below  the  MB 
layer,  the  range  normalized  signal  intensity  in  profile 
2  is  only  slightly  higher  than  that  observed  in  profile  I 
at  the  same  altitude  (-86.5  dBm  to  -77.5  dBm  instead 
of-87.0  dBm  to  -80.5  dBm).  However,  below  the  1.3 
km  altitude,  there  is  a  significant  decrease  of  radar  re¬ 
flectivity.  Also,  the  increase  of  mean  velocity  across 
the  MB  is  appreciably  greater  (1.3-5  m  s“‘  instead  of 
1 .1-4.2  m  s"').  The  greatest  difference  with  respect  to 
profile  I  is  the  Doppler  velocity  maximum  at  an  alti¬ 
tude  of  approximately  1.2  km. 
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In  the  absence  of  significant  air  vertical  velocity,  a 
mean  velocity  increase  must  relate  to  a  corresponding 
increase  in  particle  size  and,  for  Rayleigh  scatterers,  a 
radar  reflectivity  increase.  However,  inspection  of  ve¬ 
locity  profile  2  indicates  that  the  velocity  maximuni  at 
1 .2  to  is  not  associated  with  a  coire^nding  increase 
of  signal  intensity.  The  only  rational  explanation  for 
this  behavior  is  that  scatterers  produdng  ^e  reflectivity 
maximum  are  indeed  larger,  but  these  particles  have 
a  Mie  backscattering  cross  section  decreasing  with  an 
increase  of  their  size.  The  altitude  at  which  a  weaker 
signal  intensity  is  observed  coincides  with  a  mean 
Dopirfer  of  5-6  m  s"‘.  This  is  in  agreement  with  Fig. 
9  which  shows  a  decrease  of  the  backscattering  cross 
section  Sr  when  particle  size  increases  from  1.3  mm 
(terminal  velocity  5  m  s"‘)  to  1 .8  mm  (terminal  velocity 
6  m  s~').  There  is  also  a  local  qpectrd  variance  maxi¬ 
mum  at  an  altitude  of  about  1 . 1  km  which  may  be  due 
to  the  flattening  of  the  Eloppler  spectrum  associated 
with  a  depression  in  the  velocity  spectrum  center  cre¬ 
ated  by  the  decrease  of  backscattering  cross  section  for 
particle  sizes  between  1  and  2  mm.  These  preliminary 
observations  of  a  MB  observed  with  a  millimeter  wave 
Doppler  illustrate  the  capability  of  the  radar  and  also 
the  need  for  further  studies. 

Observations  were  also  made  inside  precipitating 
thunderstorms  at  Boulder,  Colorado.  A  thin  nylon 
cover  could  effectively  protect  the  radar  against  rain 
without  producing  appreciable  attenuation  of  the  94- 
GHz  signal  During  the  observations  only  light  rainfall 
occurr^,  although  sparse  4  to  5  mm  raindrops  could 
be  occasionally  observed.  Figure  16  shows  an  example 
of  updraft  structures  observed  on  a  thunderstorm 
overhang.  The  vertical  profiles  were  obtained  every  3 
seconds  and  exhibit  notable  variations  of  both  mean 
Doppler  and  signal  intensity  in  the  time  interval  shown. 
Figure  17  shows  the  time-height  presentation  of  the 
updraft  event  which  lasted  2  to  3  minutes  (see  captions 
for  a  discussion  of  the  data). 

5.  Conclusion 

A  3.2-mm  wavelength  (W-band)  is  about  the  shortest 
wavelength  that  can  be  used  for  a  meteorological 
Doppler  radar.  Radar  components  are  available  off- 
the-shelf  and  the  sensitivity  of  the  radar  for  small  size 
hydrometeors  is  comparable  to  that  of  a  K«-band  radar. 
However,  the  W-band  radar  has  the  advantage  of  re¬ 
duced  size,  weight  and  power  requirements.  The  use 
of  a  gridded  Klystron  oscillator  instead  of  a  magnetron 
also  provides  a  more  efficient  Doppler  system,  as  the 
phase  noise  is  considerably  reduced.  The  use  of  a  signal 
integrator  and  a  high  pulse-repetition  rate  (10  kHz) 
provides  an  effective  processing  gain  of  approximately 
20  dB  compared  with  single  pulse  data,  so  that  our 
radar  should  detect  a  -50  dBZ  cloud  reflectivity  at 
three  kilometer  range  and  -40  dBZ  at  10  km.  Radar 


Fig.  16.  Time-height  presentation  of  mean  Doppler  velocity  in  a 
thunderstorm  outflow,  llie  numbers  indicate  downward  motion  in 
m  s*'.  The  letters  indicate  upward  motion  with  -  I  m  s”',  J9  *  2 
ro  s*',  etc.  The  echoes  below  the  cloud  layer  are  clear  air  targets 
showing  a  primarily  upward  motion.  Note  the  aliased  velocity  above 
8  m  $■'  in  the  bottom  right  of  the  figure.  C  =  8  m  s"',  f  ■  9  m  s"‘, 
etc. 


sensitivity  estimates  derived  from  observations  of 
thunderstorm  tops  at  a  range  of  1 5-20  km  in  high  hu¬ 
midity  conditions  indicate  that  the  radar  performance 
agrees  with  the  radar  equation  in  this  paper  and  the 
radar  characteristics.  Also,  all  cirrus  layers  are  detected 
and  nearly  all  clouds  are  above  the  freezing  level.  We 
have  a  very  limited  experience  in  observing  fair  weather 
cumulus  clouds,  vvhich  seem  to  be  detected  only  if  their 
base  shows  some  sign  of  darkness  indicative  of  their 
vertical  development.  This  is,  however,  a  very  subjec¬ 
tive  estimate  and  we  are  planning  more  systematic  ex¬ 
periments  in  conjunction  with  K,-band  and  X-band 
radars  and  also  aircrafts  for  particle  sampling. 

Mie  backscattering  oscillations  and  their  effects  on 
mean  Doppler  and  reflectivity  arc  important  at  such  a 
short  wavelength  because  the  transition  from  Rayleigh 
to  Mie  scatterers  occurs  for  relatively  small  (D  ~  0.8 
mm)  particles.  At  94  GHz,  the  equivalent  radar  reflec¬ 
tivity  factor,  Z„  cannot  exceed  30  dBZ,  even  for  the 
highest  rainfall  rates.  At  3-mm  wavelength,  backscat¬ 
tering  maxima  occur  for  1. 1, 2.2  and  3.3  mm  diameter 
raindrops,  and  backscattering  minima  for  1.8  and  3.9 
mm  diameter.  The  3-mm  wavelength  is  therefore  an 
attractive  choice  for  multiwavdength-derived  drop  size 
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Ra.  17.  Selected  vertical  profiles  of  mean  Doppler  through 
downdraft  and  updraft  regions,  observed  every  5  seconds. 


distribution  measurements  if  the  W-band  observations 
are  done  jointly  with  K»  and  X-band  Doppler  radars. 
The  dropsize  distribution  can  be  measui^  using  the 
complete  spectral  information  available  at  the  three 
wavelengths  with  vertically  pointing  observations.  Such 
a  multiwavelength  experiment  can  produce  much 
needed  information  on  raindrop  size  distribution  even 
in  thunderstorms. 
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